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RESEARCH  IN  ADAPTIVE  BEAMFORMING  FOR  SATELLITE  COMMUNICATIONS 

1.  INTRODUCTION 

Reported  below  are  the  results  of  a  study  carried  out  in  the  in¬ 
terval,  November  1,  1981  -  October  31,  1982  on  adaptive  beamforming  and 
interference  cancelling  in  large  sparse  satellite-borne  arrays.  The 
application  in  mind  is  the  uplink  of  a  communication  satellite  using 
beam-switching  to  ground  users  who  may  be  close  to  potential  interferers. 

In  an  earlier  exploratory  study,  [1],  the  use  of  large  sparse  arrays  was 
examined,  giving  their  anticipated  benefits  and  describing  their  problems. 
In  brief,  the  array  size  is  made  large  to  put  potential  interferer  on  the 
quiescent  pattern  sidelobes  where  they  can  be  suppressed  effectively;  it  is 
made  sparse  to  keep  system  complexity  within  bounds.  Large  size  is  us¬ 
ually  attended  by  element  position  uncertainty  so  that  forming  accurate 
steering  vectors  toward  desired  sources  becomes  a  problem;  both,  beam 
placement  error  and  random  error  in  forming  the  pointing  vector  phasors 
were  anticipated.  A  processing  system  configuration,  including  alternative 
processing  schemes  was,  nevertheless,  found  and  this  became  the  focus 
of  the  study  reported  here. 

In  section  (2)  below  we  present  a  technical  summary  of  the  study  and 
its  results.  Detailed  reports  on  which  this  summary  is  based  are  given  in 
the  Appendices  of  this  document.  Section  3  contains  the  conclusions  and 
recommendations  for  further  study. 

2.  TECHNICAL  SUMMARY 

2.1  Hybrid  Array  Analyses 

The  principal  part  of  our  work  dealt  with  the  properties  of 
an  array  processing  scheme  which  we  term  the  "Hybrid  Array".  It  is  a 
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composite  of  the  Howells-Applebaum  processor  which  utilizes  known  direction 
of  signal  arrival,  and  the  Widrow-Compton  processor  which  uses  a  reference 
signal  extracted  from  the  incoming  field  through  its  known  superimposed 
spread-spectrum  code.  This  concept  was  arrived  at  by  first  supposing  that 
rapid  beam  switching  is  achieved  with  a  processor  that  utilizes  pointing 
vector  injection,  then,  recognizing  that  such  arrays  are  highly  sensitive 
to  pointing  vector  errors,  ways  were  sought  to  overcome  the  effect  of 
such  errors.  The  addition  of  a  self-generating  reference  circuit  which 
makes  use  of  known  signal  structure,  was  conjectured  to  be  one  such  way. 

An  alternative  scheme  which  bootstraps  a  pointing  correction  was  also 
conceived,  but  remains  to  be  analyzed.  Analyses  and  simulations  of  the 
Hybrid  Array  were  carried  to  determine  its  steady  state  and  transient 
properties,  and  the  improvement  it  provides  when  subjected  to  pointing 
error.  The  latter  includes  error  in  orienting  the  main  beam  toward  the 
desired  source,  and  independent  phase  errors  in  the  pointing  vector  tend¬ 
ing  to  defocus  the  main  beam. 

2.1.1  Steady  State  Behavior  of  Hybrid  Array 

Appendix  A  gives  the  results  of  the  steady  state  analysis  with  the 
array  processor  shown  in  Figure  1,  and  the  principal  result  shown  in 
Figure  2.  The  latter  shows  the  effect  on  the  output  signal  to  noise  ratio 
(SNR)  of  adding  the  reference  generating  circuit  to  an  Applebaum  array 
in  which  the  steering  vector  may  be  pointed  off  target,  p  is  a  measure  of 
perfection  of  the  reference  generating  circuit  and  e  is  a  measure  of  the 
pointing  error.  The  loss  in  SNR  with  pointing  error,  and  the  decreasing 
sensitivity  to  pointing  error  as  the  reference  generator  approaches  the 
ideal,  is  evident  from  these  curves.  Withp-1,  meaning  that  the  reference 
generator  perfectly  matches  the  desired  signal  output  so  that  the  residue 
Gj  (see  Figure  1)  contains  no  signal,  the  effect  of  pointing  error  is 
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essentially  the  same  as  It  would  be  in  a  focused  nonadaptive  array.  The 
Applebaum  Array  which  is  a  limiting  case  of  the  hybrid  array  with  the 
lowest  value  of  p  is  therefore  poorer  than  the  Hybrid  Array  for  any 
pointing  error. 

Figures  4  through  8  of  Appendix  A  show  the  results  of  computations 

carried  out  with  a  7  element  array  comprised  of  non-uniformly  spaced 

elements  for  a  wide  variety  of  conditions.  The  parameter  Fg  represents 

the  amplitude  and  phase  of  the  reference  relative  to  the  signal  component 

entering  the  reference  generator;  F  -1  is  ideal.  Typically,  one  might 

s 

expect  the  reference  amplitude  to  be  a  few  percent  away  from  ideal  and 
several  degrees  off  in  phase.  For  such  departures  from  ideal  the  curves 
all  show  a  substantial  reduction  of  sensitivity  to  pointing  error  over 
the  Applebaum  Array  (the  Fg*0  case).  It  should  be  noted  that  cases  of 
multiple  inter ferers  and  interference  inside  the  main  lobe  of  the  quies¬ 
cent  array  are  included  in  these  results. 

The  observed  behavior  can  be  explained  by  examining  the  analytical 
results.  Without  reference,  the  array  processor  treats  the  desired  signal 
as  interference  when  pointing  is  inaccurate.  With  perfect  reference,  the 
mechanism  which  marks  the  incoming  signals  for  suppression  (the  co- 
variance  matrix  of  the  array  element  outputs,  less  a  constant  times 
the  covariance  matrix  of  the  desired  signal  alone)  has  the  signal  component 
removed  by  the  action  of  the  reference.  As  a  consequence  there  is  no 
tendency  to  suppress  it  when  there  is  pointing  error,  though  there  is  the 
normal  loss  that  one  would  encounter  in  a  miss-aimed  conventional  beam- 
forming  array.  With  imperfect  reference,  some  vestige  of  signal  remains 
in  the  suppression  mechanism.  The  tendency  to  suppression  is  blunted  but 
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not  totally  eliminated.  The  conclusion  to  be  drawn  form  these  results  is 
that  the  Hybrid  Array  is  a  significant  improvement  over  the  directionally 
constrained  array  though  it  may  fall  short  of  ultimate  perfection  as  it 
stands. 

2.1.2  Transient  Behavior  of  Hybrid  Array 

The  results  of  analysis  and  simulation  of  the  transient  behavior  of 
the  Hybrid  Array  are  given  in  Appendix  B.  Because  there  are  two  filters  in 
the  Hybrid  Array,  one  in  the  sidelobe  canceller  weight  setting  loop  and 
one  in  the  reference  generator,  two  coupled  set  of  differential  equations 
result.  Both  filters  are  assumed  to  be  of  first  order  in  this  analysis. 

By  assuming  that  the  wei^It  setting  loop  is  much  slower  than  the  re¬ 
ference  generator,  one  set  of  equations  can  be  solved  without  involving 
the  other.  The  result  is  that  a  single  set  of  first  order  linear  differ¬ 
ential  equations  with  time  variable  coefficients  is  obtained  for  the 
time  variable  weights  in  the  sidelobe  canceller.  Explicit  solution  of 
the  output  transient  SINR  for  an  H-elemeant  array  remains  to  be 
carried  out.  We  have  however  dealt  more  completely  with  the  two-element 
array.  The  instantaneous  values  of  the  single  complex  weight  was  obtained 
in  this  case.  From  this,  one  can  with  some  calculation  determine  the 
response  as  a  function  of  time. 

Simulations  were  carried  out  on  the  two-element  array  with  A/2  spacing 
assuming  a  signal  and  one  interferer  present.  Results  were  obtained  with 
and  without  noise.  The  signal  was  taken  to  be  modulated  with  a  spread 
spectrum  code  and  the  interference  was  taken  to  be  a  pure  sinusoid.  The 
array  processor  was  assumed  perfectly  synchronized  to  the  spectrum  spreading 
code  so  that  the  transient  behavior  observed  concerns  only  the  response  of 
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Che  weight  setting  loop  and  the  reference  generator.  Runs  were  made 
with  various  sets  of  parameter  values,  input  amplitudes,  and  pointing 
angles.  Parameter  values  could  be  found  with  which  convergence  of  the 
signal  to  interference  (and  noise,  where  applicable)  ratio  occurs  in 
a  fraction  of  the  bit  interval  provided  the  pointing  angle  is  not  greatly 
in  error.  We  found,  for  instance,  that  pointing  errors  up  to  5°  could  be 
tolerated  but  with  a  20°  error  convergence  was  slowed  to  the  point  where 
the  output  signal  to  interference  ratio  never  reaches  an  adequate  level. 

In  our  simulations  idealized  signal  pulses  are  used  with  zero  rise 
time.  Clearly  the  filter  in  the  reference  generator  will  not  follow  such 
an  input  so  that  at  the  beginning  of  each  signal  polarity  reversal  a  sub¬ 
stantial  difference  exists  between  signal  input  and  output  of  the  reference 
generator.  This  reflects  itself  as  a  sudden  drop  in  signal  to  interference 
ratio  of  the  array  output.  The  processor  overcomes  this  effect  in  less 
than  half  the  pulse  interval  in  the  examples  simulated,  but  this  effect 
will  reduce  the  receiver  output  quality.  The  loss  can  be  overcome  by  an 
increase  in  transmitted  power  of  the  order  of  3dB.  Or,  the  weight  obtained 
at  the  end  of  the  first  signal  interval  can  be  frozen  for  the  duration  of 
the  transmission.  The  latter  approach  may,  however,  not  be  acceptable  in 
the  case  of  blinking  interferers  where  nulls  may  have  to  be  moved  during 
the  brief  interval  of  a  users  on-time. 

As  a  check  on  the  analytical  work,  comparisons  were  made  between  the 
transient  weight  variation  obtained  by  analysis  and  by  simulation.  The 
results  matched  very  well  suggesting  that  certain  approximations  used  in 
the  analysis  are  valid. 
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2.1.3  Effect  of  Random  Errors  in  Pointing  Vector  on  Steady  State 
Behavior 

In  Appendix  C  an  analysis  is  reported  on  the  effect  of  random  phase 
and  amplitude  errors  on  the  steady  state  signal  to  interference  plus  noise 
ratio  of  the  Hybrid  Array.  The  analysis  reported  previously  [1]  on  the 
effect  of  pointing  error  was  concerned  with  correlated  deviations  in  the 
phases  of  the  pointing  vector  components  which  result  in  a  shift  of  the 
direction  of  focus.  This  sort  of  error  will  occur  where  there  is  un¬ 
certainty  in  the  signal  arrival  direction.  For  the  case  treated  now  the 
direction  of  arrival  is  assumed  known  correctly  but  the  amplitudes  and 
phases  of  the  components  of  the  pointing  vector  generated  by  the  processor 
are  assumed  subject  to  independent  errors.  Such  errors  may  be  a  consequence 
of  circuit  inadequacies,  quantization  errors  when  digital  processing  is  done, 
and  element  position  uncertainties.  In  section  2.2  below  we  will  discuss 
an  analysis  of  a  method  for  estimating  the  pointing  vector  by  calibrating 
to  a  pilot  station  then  steering  to  the  desired  ground  user.  It  is 
inevitable  that  this  estimation  procedure,  like  any  other  such  procedure, 
is  less  than  perfect  and  the  principal  point  of  that  study  is  to  find  the 
error  In  the  pointing  phases  so  determined. 

Results  were  obtained  from  which  the  SINR  as  a  function  of  phase 
and  amplitude  standard  deviations  can  be  obtained;  the  reference  generator 
quality  parameter,  Fg,  and  the  number  of  elements,  K,  are  principal  para¬ 
meters  in  the  results.  The  analysis  was  carried  out  for  a  version  of 
the  Hybrid  Array  derived  from  a  variant  of  the  Applebaum  Array  in  which 
the  signal  output  magnitude  is  not  constrained  (See  Figure  2  of  Appendix 
C  for  this  version;  the  version  based  on  the  constrained  Applebaum 
Array  is  shown  in  Figure  1  of  Appendix  C.).  Calculations  were  carried  out 
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for  che  ratio  of  mean  signal  power  output  to  mean  noise  power  output 
(i.e.,  interference  was  not  included)  as  a  function  of  the  standard 
deviation  of  amplitude  and  phase  for  arrays  of  7,  30,  and  100  elements 
and  for  values  of  Fg  ranging  from  1  (perfect  reference)  to  0  (no  re¬ 
ference,  identical  to  Applebaum  Array).  The  results  shown  in  Figures 
3-8  of  Appendix  C  were  similar  to  those  obtained  for  the  case  of  error 
in  pointing. 

With  no  reference  the  loss  of  SNR  is  substantial  and  with  perfect 
reference  the  loss  is  zero.  With  reference  within  5%  in  amplitude  and 
10°  in  phase  the  loss  is  generally  small  enough  to  be  acceptable  in  the 
case  of  the  7  element  array.  Of  great  interest,  though,  is  the  effect 
of  increasing  the  number  of  array  elements.  The  sensitivity  to  the 
random  errors  increases  with  increases  in  number  of  array  elements.  In 
particular,  in  the  case  of  a  30  element  array,  the  loss  is  substantial 
unless  the  reference  generator  is  close  to  perfect,  or  the  errors  in 
phase  and  amplitude  are  very  small.  We  point  out  though  that  the  im¬ 
provement  over  the  Applebaum  Array  is  great  in  all  cases  of  random  pointing 
error. 

Because  the  analysis  was  carried  out  for  the  unconstrained  version 
of  the  Hybrid  array  (as  shown  in  Figure  2  of  Appendix  C)  we  compared 
calculated  results  for  SNR  of  the  7  element  array  for  unconstrained  and 
constrained  versions  (the  latter  as  in  Figure  1  of  Appendix  C) .  The 
differences  were  found  to  be  only  slight  so  that  we  are  inclined  to  view 
the  other  results  obtained  as  being  approximately  applicable  to  both 
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2.2  Pointing  Vector  Estimation 


We  examined  a  method  for  estimating  the  pointing  vector  to 
a  desired  source  when  interferers  may  be  close  to  the  desired  source  and 
the  array  element  positions  are  not  perfectly  known.  The  procedure  is 
based  on  the  beamforming  and  scanning  processes.  The  analysis  of  the 
method  appears  in  Appendix  D,  below. 

First,  the  array  generates  a  pointing  vector  at  the  reference  beacon 

by  the  use  of  a  reflector  antenna.  The  assumption  is  made  that  interferers 

are  seen  by  this  antenna  only  through  its  sidelobes  and  that  sidelobe 

levels  are  reasonably  low.  The  pointing  vector  is  determined  by  correlating 

array  outputs  with  the  output  of  the  reflector  antenna.  The  vector  obtained 

in  this  way  has  a  bias  caused  by  some  residual  interference  at  the  reflector 

antenna  output  correlating  with  the  interference  at  array  element  outputs. 

There  is  also  a  random  component  which  is  averaged  out  if  the  integration 

time  is  long  enough.  Analysis  shows  that  both  of  these  errors  can  be 

treated  as  random  phase  errors.  For  integration  time  long  enough,  the 

-4 

standard  deviation  of  those  random  phase  errors  is  about  5x10  radian  if 
there  is  one  interferer  with  the  same  power  as  the  beacon  source  and  if 
the  sidelobe  reduction  is  30dB. 

The  second  step  is  to  steer  the  array  from  the  beacon  source  to  the 
desired  ground  source  by  using  the  information  on  angular  displacement 
between  them  and  the  imperfect  information  on  array  element  positions.  The 
total  error  in  the  pointing  vector  aimed  at  a  desired  user  is  the  sum  of 
the  error  in  beamforming  and  that  in  steering.  For  a  seven  element  10X 
Hybrid  Array  with  Fg*0.95  and  SNR  *  10  at  input,  the  tolerance  of  element 
position  is  about  0.44X  for  ldB  output  SINR  loss  if  the  error  in  beamforming 
on  the  beacon  is  5x10  ^  radian.  If  the  error  induced  in  steering  can  be 
eliminated,  the  error  of  5x10  ^  radians  causes  a  negligible  loss  in  SINR 

for  F  *0. 95e~'* 10  . 

s 
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2.3  LMS  Interference  Canceller  with  Reference  Generator 

Appendix  E  contains  an  analysis  of  an  LMS  adaptive  array 
processor  which  utilizes  a  self-generated  reference  only;  that  is,  one 
which  does  not  use  direction  of  signal  arrival  information.  An  array 
processor  to  accomplish  this  is  described  by  Compton  (reference  [3], 
Appendix  E)  and  the  problem  of  weight  cycling  encountered  with  this 
scheme  is  analyzed  by  DiCarlo  and  Compton  (references  [4],  [5],  Appendix 
E) .  We  have  here  analyzed  two  variants  of  the  Compton  scheme  both  of 
which  differ  from  the  original  in  that  one  element  is  left  unweighted  and 
limiters  in  the  reference  generators  are  omitted.  The  difference  between 
the  two  schemes  analyzed  here  resides  in  how  the  reference  generator  is 
driven. 

It  is  shown  in  Appendix  E  that  simple  and  attainable  conditions  exist 
for  the  weights  to  converge  to  a  constant  value.  Also,  the  SNR  for  the 
two  schemes  devised  is  obtained.  In  one  scheme,  where  the  reference  genera¬ 
tor  Is  driven  by  the  output  of  all  array  elements,  the  SNR  turns  out  to  be 
generally  low  -  as  low  In  some  cases  as  the  SNR  of  a  single  element.  The 
other  realization,  in  which  the  reference  generator  is  driven  by  the  out¬ 
put  of  the  weighted  elements  onljv  the  output  SNR  is  the  full  value  ex¬ 
pected  out  of  (N-l)  coherently  combined  element  outputs;  N  is  the  number 
of  array  elements. 


We  point  out  that  the  earlier  analyses  by  DiCarlo  and  Compton  and 
the  analysis  of  Appendix  E  were  carried  out  without  interference  present. 


2.4  Effect  of  Interference  on  the  Behavior  of  the  LMS  Array 
with  Internally  Generated  Reference 

As  pointed  out  in  section  2.3  above  the  published  analyses  on  the 
subject  of  adaptive  arrays  with  internally  generated  reference  were 
carried  out  with  a  single  signal  and  Independent  element  noise  only;  that 
is,  interference  was  not  considered.  In  Appendix  F  the  earlier  analyses 
are  extended  to  include  inter ferers. 

The  array  processing  scheme  of  DiCarlo  (reference  [4]  Appendix  F) 
is  assumed,  along  with  the  phase  shift  compensation  scheme  of  Bar-ness 
(reference  [3]  Appandix  F).  The  analytical  difficulties  which  arise  with 
interference  present  is  associated  with  the  greater  complexity  of  the 
covariance  matrix  with  more  than  one  signal  present.  The  covariance 
matrix  of  the  multiple  signals  is  of  rank  equal  to  the  number  of  signals. 

The  uncoupling  of  the  resulting  simultaneous  differential  equations  presented 
difficulty  though  a  transformation  was  found  which  made  it  possible  to 
separate  the  effects  of  desired  signal  and  interference. 

It  is  here  shown  that  the  phase  shift  compensation  is  effective  in 
the  multiple  interference  case  and  that  weight  cycling  can  be  avoided.  The 
conditions  to  be  satisfied  though  have  not  yet  been  found  for  the  general 
case.  One  expects  the  conditions  and  the  transient  behavior  o{  the  array  to 
be  affected  by  the  relative  levels  of  the  incoming  signals  and  on  their  ar¬ 
rival  directions.  It  was  possible  however  to  get  an  explicit  condition  for 
convergence  in  the  case  of  a  single  interferer  in  addition  to  the  desired 
signal. 

Further  work  is  to  be  done  on  the  conditions  for  convergence  in  multiple 
interference  as  well  as  on  convergence  rates  and  achievable  signal  to 
interference  plus  noise  ratios. 
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3. 


CONCLUSIONS  AND  RECOMMENDATIONS 


We  have  seen  Chat  the  Hybrid  array  processor  is  a  substantial 
improvement  in  the  steady  state  over  an  array  using  only  direction  of 
arrival  information  when  systematic  or  random  errors  are  encountered  in 
the  focusing  phases  and  amplitudes.  The  transient  analysis  also  suggests 
that  if  signal  structure  alone  is  used  to  drive  the  adaption  process 
the  response  time  will  be  greater  than  with  the  Hybrid  Array.  We  point 
out  that  the  Hybrid  Array  as  described  here  does  not  suffer  from  the 
weight  cycling  problem  associated  with  some  processors  which  utilize  known 
signal  structure  to  self-generate  a  reference.  Though  we  have  not  dealt 
with  the  problem  of  code  aquisition  time  in  this  work  it  is  to  be  noted 
that  the  use  of  a  pointing  vector,  even  an  approximate  one,  insures  a 
reasonably  large  signal  component  immediately  from  which  the  code 
aquisition  circuit  can  more  readily  extract  synchronizing  Information, 

We  have  analyzed  a  pointing  vector  estimation  method  based  on  the 
use  of  a  pilot  source  sufficiently  out  of  the  way  of  interferers.  We 
have  derived  the  estimation  errors  and  with  the  pointing  error  sensitivity 
of  the  Hybrid  processor,  we  have  demonstrated  feasibility  of  the  scheme. 

Finally,  we  have  extended  previous  work  on  adaptive  processing 
utilizing  signal  structure  alone  to  self-generate  a  reference.  Schemes 
were  devised  and  analyzed  which  can  be  prevented  from  generating  oscillating 
weights,  and  the  behavior  of  such  arrays  in  the  presence  of  interfering 
sources  has  been  analyzed. 

The  work  carried  out  leads  us  to  make  a  number  of  recommendations 
for  further  study. 
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While  the  Hybrid  Array  promises  significant  improvement  over  earlier 
schemes  it  depends  on  being  able  to  generate  a  good  reference.  What  is 
needed  is  either  a  mechanism  for  insuring  that  the  self-generated  reference 
is  near  perfect  or  a  mechanism  for  correcting  the  steering  vector.  A 
possible  circuit  to  accomplish  this  by  the  addition  of  another  steering 
phase  correction  feedback  loop  for  a  two-element  configuration  had  been 
conceived  before  [1].  Variants  to  accomplish  this  same  end  have  been  devised 
since  then.  These  schemes  should  be  further  pursued.  It  may  also  be 
valuable  to  pursue  methods  of  steering  correction  in  one  of  the  Applebaum 
arrays  -  that  is,  with  processors  which  do  not  generate  a  reference  in¬ 
ternally.  In  addition,  or  as  an  alternative,  ways  should  be  devised  for 
controlling  the  reference  generator  output  through  a  feedback  loop  to 
minimize  the  signal  component  in  the  residue. 

The  improvements  obtained  with  the  Hybrid  Array  can  be  viewed  as  a 
consequence  of  utilizing  more  information  about  the  desired  signal  than 
heretofore:  the  general  problem  of  how  to  make  optimal  use  of  other  bits 
of  prior  information  is  worth  considering.  For  instance,  approximate 
location  of  interferers,  or  their  signal  structure  may  also  be  known. 

There  is  also  additional  prior  information  regarding  the  desired  signal 
which  is  potentially  useful  -  the  waveforms  and  error  control  codes  used, 
for  instance.  This  also  leads  to  the  question  of  how  to  design  signals 
which  will  most  effectively  aid  the  array  focusing  mechanism. 

We  had  originally  adopted  the  view  that  the  array  size  to  be  used  would 
be  sufficiently  large  to  put  potential  interferers  out  of  the  main  beam  of 
the  quiescent  array  pattern. 
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We  have  nevertheless  obtained  some  computational  results  for  the  Hybrid 
Array  assuming  interferers  veil  Inside  the  beamvidth.  Naturally,  the 
SINR  was  found  to  be  adversely  affected.  There  is  reason  to  expect 
that  the  SINS  with  main-beam  interference  is  to  some  extent  dependent 
on  array  element  layout.  The  problem  of  element  layout  was  one  of  the 
items  considered  originally  by  us  but  ruled  out  as  a  significant  factor 
for  large  arrays  vith  no  main-beam  interference.  Since  close  interfere 
cannot  be  ruled  out  totally  the  effect  of  main-beam  interference  and 
its  dependence  on  element  layout  ought  to  be  considered. 

The  schemes  with  which  we  have  dealt  are  all  narrowband;  that  is, 
they  utilize  complex  non-frequency  dependent  weights  prior  to  element  output 
combining.  As  long  as  the  time  displacement  of  signals  seen  by  the  different 
elements  is  small  relative  to  the  inverse  of  the  signal  bandwidth  this 
method  is  satisfactory.  Tor  broadband  signals,  and/or  for  wavefronts 
arriving  at  angles  very  far  off  array  broadside  this  may  not  be  the  case. 

Ways  of  extending  the  methods  here  described  to  wideband  signals  would  be 
useful. 

One  of  the  principal  advantages  we  see  in  arrays  utilizing  directional 
Information  in  addition  to  signal  structure  information  is  the  potential 
for  high  speed  beam  switching.  We  have  verified  fast  convergence  of 
the  Hybrid  Array  with  an  interferer  present  in  a  particular  location. 

However,  the  functioning  of  the  array  in  a  rapidly  changing  dynamic  en¬ 
vironment  of  multiple  desired  sources  and  dispersed  blinking  interferers 
requires  further  attention.  One  can  imagine  a  simulation  experiment  with 
such  a  scenario.  Also  deserving  attention  for  multiple  access  applications 
is  the  use  of  multiple  directional  constraints  in  which  two  or  more 
simultaneous  beams  are  focused. 
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Appendix  A 


STEADY  STATE  BEHAVIOR  OF  HYBRID  ARRAY 
by 

F.  Haber,  Y.  Bar-ness,  C.  C.  Yeh 

Abstract 

The  steady  state  properties  of  an  adaptive  array  utilizing  prior  know¬ 
ledge  of  both  approximate  signal  arrival  direction  and  signal  characteris¬ 
tics  are  here  presented.  The  method  combines  the  features  of  a  directionally 
constrained  array  and  one  with  a  self-generated  reference  signal.  Explicit 
results  are  obtained  for  output  signal,  interference,  and  noise  powers 
assuming  a  single  lnterferer  is  present.  The  inclusion  of  a  self-generated 
reference  circuit  is  shown  to  reduce  the  sensitivity  to  pointing  error  typical 
of  arrays  utilizing  a  zero  order  directional  constraint,  the  improvement  being 
a  consequence  of  the  reduction  of  the  desired  signal  component  feel  back  to  the 
sidelobe  cancelling  circuit.  A  relationship  between  the  degree  of  sensi¬ 
tivity  reduction  and  the  quality  of  the  reference  signal  is  developed. 
Results  of  computations  of  signal  to  interference  plus  noise  ratios  for  a 
7-element  10  wavelength  non-uniformly  spaced  array  as  a  function  of  pointing 
error  are  presented.  These  results  show  the  behavior  with  one  interferer 
inside  and  outside  the  beamwidth  of  the  quiescent  array,  and  with  multiple 
interferers,  for  various  degrees  of  perfection  of  the  reference  generating 
circuit.  In  all  cases  the  computations  confirm  that  the  otherwise  severe 
effects  of  small  pointing  errors  are  substantially  reduced. 

This  research  is  supported  by  Rome  Air  Development  Center  under  Contract 
No.  F30602-81-K-0211. 
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INTRODUCTION 

Adaptive  arrays  suitable  for  use  in  conmunication  systems  can  be 
categorized  according  to  the  prior  information  utilized  to  distinguish 
desired  from  undesired  signals.  One  scheme  described  by  Applebaum  and 
Chapman  [1]  is  suited  to  point-to-point  communication  where  direction  of 
signal  arrival  is  known;  the  directional  information,  through  an  input 
steering  vector,  constrains  the  array  gain  in  this  direction,  maximally 
rejecting  other,  unwanted  sources.*  Another  scheme,  based  on  a  concept 
by  Widrow,  et.  al.  [2]  recognizes  the  desired  signal  via  an  externally 
supplied  reference.  For  communication  use,  an  extension  of  the  latter 
was  made  by  Compton,  et.  al.  [3,  4]  wherein  the  reference  is  internally 
generated  utilizing  some  prior  known  signal  structure  information 
(e.g.,  the  spread  spectrum  code  superimposed  on  the  signal). 

One  of  the  principal  shortcomings  of  the  first  method,  the  directionally 
constrained  method,  is  that  small  pointing  errors  result  in  large  losses 
in  the  output  signal  to  interference  plus  noise  ratio  (SINR).  An  analysis 
of  this  effect  was  made  by  Compton  [5]  and  was  reinforced  by  results  obtained 
by  us  [6].  Unfortunately,  pointing  error  is  endemic  in  many  applications 
where  direction  of  arrival  is  ostensibly  known.  The  second  method  is 
however  not  free  of  difficulties  among  which  are  problems  relating 
to  imperfect  reference  generation  [7,  8],  Furthermore,  one  expects  better 
transient  response  using  the  directionally  constrained  scheme  because 
a  substantial  signal  power  becomes  immediately  available  at  the  array  output 
on  injection  of  an  approximately  correct  steering  vector. 

*  With  a  constraint  on  the  gain  in  a  chosen  direction  the  array  is  said  to 
be  zero  order  directionally  constrained  [1]. 
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A  comparison  study  of  these  two  methods  [9]  suggested  that  the  two 
can  be  made  to  complement  one  another  and  has  led  us  to  look  at  hybrid 
schemes  which  take  advantage  of  both  kinds  of  prior  information — 
directionality  and  signal  structure.  We  present  below  analyses  and 
computational  results  of  the  steady  state  properties  of  one  such  scheme 
as  defined  in  Figure  1.  The  array  processor  is  represented  in  a  form 
similar  to  that  in  Applebaum  and  Chapman  [1,  Figures  2  and  A]  with  the 
addition  of  a  reference  generating  loop.  We  treat  the  case  of  narrow- 
band  information  bearing  signals  for  both  desired  and  undeslred  arrivals 
so  that  steering  ana  nulling  operations  are  carried  out  by  control  of 
amplitudes  and  phases  of  s^  and  £  respectively,  at  band  center.  Further¬ 
more,  we  assume  the  reference  generating  circuit  is  operating  in  the 
synchronous  mode  with  respect  to  the  desired  signal;  the  acquisition  circuit 
for  code  timing  is  therefore  not  shown. 


2. 


ANALYSIS 


The  array  system  under  consideration  is  shown  in  Figure  1.  Output 
signals  from  the  array  elements  are  represented  by  the  complex  vector 


v(t) 


J 


*  l 

J-l 


8j(t)a  +  n(t) 


(1) 


where  s  ,  is  the  arrival  phase  vector  of  the  desired  signal  comprised  of  N  unit 
— d 

amplitude  components,  a (t)  is  its  complex  envelope;  s^  is  the  arrival  phase 


vector  of  the  j — -  interference  signal,  6j  (t)  is  its  complex  envelope; 

and  n(t)  is  the  complex  noise  envelope  vector,  the  components  of  which 

are  assumed  independent.  If,  for  instance,  the  elements  are  arranged 

along  the  x-axis  at  positions  x^,  i  =  1,2,...N,  a  desired  signal  arriving 

as  a  plane  wave  with  its  plane  of  arrival  at  an  angle  9  relative  to  the 

a 

x-axis  induces  signal  components 


a(t)sdi  *  a(t)exp(jkxicos9d),  i  -  1,2,...N  (2) 

in  the  array  elements.  s^  =  expCjkx^ose^  is  the  i—  unit  amplitude  component 

of  s,.  The  rf  wave  at  the  i—  element  is  given  by  v.  (t)exp(jcot  +  <|>).  The  de- 
U  til  ^ 

sired  signal  component  at  the  i  element,  as  represented  by  (2),  has  the 

carrier  exponential  suppressed. 

The  beamformer  defined  by  the  steering  vector,  £*  in  Figure  1  generates  an 
output  given  by  the  inner  product  of  and  v  given  by  (1);  that  is. 


e 

m 


T* 

s  v 


(3) 


(the  time  variable  in  e  and  v  is  omitted  but  implied).  The  output  of  the 

m  — 


A- 5 


system  is  given  by 


e 

o 


T 

e  -  y  u 
m  ^ — 


where  is  the  weight  vector  on  the  sidelobe  canceller  and 


(4) 


u  =  Av 


(5) 


A  is  an  (N  -  1)  x  N  matrix  of  rank  (N  -  1)  chosen  such  that  Ai3  =  0^;  and,  as  a 
consequence,  such  that  signals  arriving  with  direction  vector  s^  do  not  con¬ 
tribute  to  u.  (3),  (4),  and  (5)  together  give 


X*  T  T  T 

e„  =  s  v  -  Av  =  (s*  -  A  y)  v 


(6) 


Following  Applebaum  and  Chapman  [1]  the  equivalent  weight  vector  w  is  defined  by 


T 

e  =  w  v 
o - 


(7) 


so  that  from  (6) 


w  =  s*  - 


(8) 


Because  of  the  way  A  is  chosen  this  weight  vector  w  insures  the  response  in 

T 

the  direction  of  That  is,  premultiplying  (8)  by  a  we  have 


T  T  *  T.T 
s  w  -  s  s  -  £  A  £  ”  N 


(9) 
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t  T  T  t  * 

In  (9)  we  have  used  a  A  £  “  £  As  »  0,  and  a_  _s  •  N,  the  number  of  elements. 

The  weight  vector  of  the  sidelobe  canceller  is  formed  via  a  circuit 
which  typically  approximates  the  LMS  algorithm  using  steepest  descent 
search.  In  the  steady  state  that  weight  vector  is  given  by 


y  -  g  u*(eQ  -  e  )  (10) 

where  er  is  the  reference  signal  as  identified  in  Figure  1  and  g  is  the 
gain  in  the  weight  setting  loop.  The  overbar  stands  for  expected  value. 

The  reference  signal  is  assumed  formed  from  the  array  output  eQ  by 
despreading  using  the  desired  signal' s  spread  spectrum  code,  filtering, 
and  respreading  using  the  same  code.  In  the  process,  inter ferers  in 
whatever  form  they  initially  appear  (i.e.,  wideband  or  narrowband)  are 
reduced  to  a  broadband  low  power  spectral  density  noise.  e^_  is  therefore 
viewed  as  being  comprised  of  a  component  closely  similar  to  the  desired 
signal  and  a  noise  component,  denoted  n^  ;  that  is. 


e 

r 


F  e  +  n 
s  od  r 


(ID 


where  e^  is  the  signal  component  of  eQ  and  Fg  is  an  operator  representing 
the  effect  of  the  bandpass  filter  on  the  desired  output  signal  component. 
We  will  specialize  it  in  our  work  to  a  complex  constant  implying  a  phase 
and  amplitude  shift  only. 
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Utilizing  these  preliminaries,  the  equivalent  weight  vector  from  which 


the  output  response  is  determined  by  (7),  is  shown  in  Appendix  A  to  be  given 


by  (see  A-12,  A-13), 


where 


-1  * 

w  ■  s 


M  -  a  F  M, 
1  s  d 


M  T  -1  *  ^ 

£  £ 

M  and  are  defined  by  (15)  and  (16)  below.  When  F  =0  the  hybrid  array  re¬ 
duces  to  the  conventional  directionally  constrained  array  and  the  weight 

-1  * 

vector  as  given  by  (12)  reduces  to  w  =  uM  £  . 

We  now  examine  the  implications  of  the  change  brought  about  by  inclusion 
of  the  reference  generating  circuit.  From  (1)  we  write  for  the  covariance 
matrix  of  the  inpucs 


M  *  v  vT  -  | a  |Md  +  £  I B j I MIJ  +  Mn 


where 


*  T 

Md  = 


and  MIj  "  4Ai  (1J 

M  =  n*nT  =  a2 I  (V 

n -  n 

(15)  is  obtained  using  the  assumed  independence  of  all  separately  arriving 
signals  (desired,  undesired  and  noise)  and  the  independence  among  the  noise 
components.  With  as  given  by  (13)  we  see  that  the  effect  of  the  self¬ 
generated  reference  signal  is  to  alter  the  contribution  of  the  desired  signal 
component  to  the  covariance  matrix  by  a  factor  (1-Fg).  When  Fg  =  1,  has, 
in  fact,  no  component  due  to  signal.  We  point  out  that  the  field 


environment  is  communicated  to  the  array  processor  via  the  matrix 
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M  in  the  conventional  directionally  constrained  array.  The  tendency  of 
the  array  processor  is  to  suppress  all  signals  recorded  in  the  matrix  M — 
except  for  the  signal  coming  from  the  constraint  direction  defined  by  the 
steering  vector  £.  If  £  is  different  from  the  direction  of  the  desired 
signal  represented  by  s^  the  processor  deals  with  it  as  if  it  were  an  unwanted 
signal,  suppressing  it  substantially  even  for  small  differences  of  angle. 

The  effect  of  removing  the  desired  signal  contribution  to  M,  or  of  reducing 
its  influence  on  M  by  making  Fg  close  to  unity,  is  to  reduce  the  tendency 
to  suppression  of  the  desired  signal  when  £  and  s^  are  not  coincident.  This 


may  be  seen  most  easily  by  imagining  only  a  single  signal  present  with 

arrival  direction  inherent  in  s.  and  with  F  -  1.  Then 

— d  s 


M  -  M  -  o2I 
Inn 


and  the  weight  vector  is  from  (22) 


w  » 


* 


s 


a?) 


The  weight  vector  is  identically  that  of  a  conventional  bearafortner  and  the 
desired  signal  will  produce  an  output 


6o(sig) 


T  * 
s  ,s  a 


(20) 


The  pattern  of  e0(g^gj  as  a  function  of  angular  difference  between  pointing 
vectors  s^  and  £  will  have  the  usual  mainbeam  beamwidth  of  approximately  \/l 
where  l  is  the  linear  dimension  of  the  array.  We  later  show  computational 
results  bearing  out  this  effect. 
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We  now  return  to  (12)  for  the  purpose  of  putting  it  into  a  more  explicit 
form  in  terms  of  the  signal,  interference,  and  noise  components.  The  result 
will  allow  us  to  draw  additional  conclusions  about  the  array  behavior  including 
the  effect  of  imperfect  elimination  of  the  signal  component  from  M.  Using  (13) 
and  (15)  we  have,  assuming  a  single  interferer  denoted  B(t)s^, 


(1  -  Fg)|a‘ 


|Md  + 


I  s 2 ! mt  +  o2i 

1  I  n 


-  “iw  -  Wd  +  yiMi  +  11  <21> 


where 
measures 
Append ix 


(i-gV^^)Md-(1-VviWV^[H>^f1)1Mi 

[ (1-Fs) YdU+ll (YIN+l)-(l-Fs)ydY j | S^S*| 2 


ty  ry  .  2  |  2 

=  |a  |  /a  measures  the  SNR  of  the  desired  signal  and  Yj  •  !  B  I  /<*n 
the  SNR  of  the  undesired  signal  at  each  element.  It  is  shown  in 
B  that  the  inverse  of  (21)  is  given  by 


To  see  the  effect  of  imperfect  reference  generation  we  specialize  (22)  to  the 
case  of  desired  signal  plus  noise  alone  so  that  Yj  *  0.  Then  (22)  becomes 


1 


-  Fs)YdMd 

(1  -  F  )Y.N  +  1 
s  a 


The  weight  vector  is  then,  using  (12), 


(23) 


w 


i  l£- 

a 

n 


w*] 


(24' 
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where 


(r  -  Fs)yd 

(1  -  Fs)YdN  +  1 


(25) 


The  output  signal  to  noise  ratio  is 


SNR 


E[|asJ(s*  -  kjMjS* ) | 2 ] 
E[|nT(s*  -  kjMjS* ) | 2 ] 


l*dV  -  klMds*)|2 

'd  111*  -  k^/l2 


(26) 


II  ‘  II  represents  the  norm  of  the  vector  inside  the  bars. 

T  * 

Making  use  of  (16)  and  of  the  fact  that  s^j^  =  N  the  numerator  of  the 
ratio  in  (26)  becomes 


U  - 


V»|2| 


T  *,2 

V-  1 


and  the  denominator  of  that  ratio  becomes 


N[1  +  (| 1  -  kjNl2-  1) I sjS*  j  2 / N2 ] 


so  that  (26)  can  be  written 


SNR 


2  i  T  *  i  2 


I1  "  kjN|  1 5^5.' 


N[1  +  ( J 1  -  kLNr  “  Dl^js  |Z/NZ] 


(27) 


This  result  embodies  both  the  degree  of  imperfect  pointing  and  imperfect 
reference.  For  the  case  where  the  reference  generator  is  perfect  k^  ■  0 
so  that 


A-ll 


(28) 


SNR 


T d i  T  *.2 


I'**!- 


The  maximum  of  this  quantity  is  obtained  when  pointing  at  s^  and  in  this 

*r  ft  j  t  *•  2  2 

case  1 s^s  |  «  I 3j8jl  "  N  so  that 

SNR  -  Nyd  (29) 


an  expected  result. 

For  the  case  where  s^  *  £,  but  with  imperfect  reference  (27)  also  gives 

SNR  -  Ny  . 

a 

just  as  one  would  expect  from  perfect  steering. 

f  *  2 

The  pointing  error  is  embodied  in  the  quantity  |  s^s  [  which,  as  stated 
above,  has  a  maximum  value  given  by  N2  when  s^  -  A  convenient  indicator 
of  pointing  error  is  obtained  by  writing 

l4/l 2  »  (1  -  ON2  O0) 

with  0  e  <  1  being  the  measure  of  pointing  error.  Using  (30)  and,  for 
brevity,  denoting 


_ 1 _ 

1(1  -  F8)YdN+  ll2 


(31) 


where  is  given  by  (25),  (27)  becomes 


SNR  -  y.N[ 
a 


(1  -  e)p  1 
c  +  (1  -  e)p 


(32) 
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The  right  factor  in  brackets  which  we  might  term  the  normalized  SNR  shows 
how  the  SNR  varies  with  pointing  error  as  measured  by  e,  for  degrees 
of  imperfection  in  the  reference  signal  as  measured  by  P.  Me  show  in 
Figure  2  the  behavior  of  that  factor  denoted  by 


(1  ~  e)p 
e  +  (1  -  e)p 


(33) 


for  p  in  the  range  (0.01,  1).  The  larger  values  of  p  represent  closer 

tracking  by  the  reference  generator  and  vice-versa.  We  see  here  clearly 

the  increasing  sensitivity  to  pointing  error  as  p  gets  smaller.  If,  for 

_2 

instance,  F  *■  0.9,  N  *  10,  «  10  we  get  p  =  10  .  It  should  be  noted 

that  with  Fg  -  0,  i.e.  when  no  reference  is  used,  the  corresponding  value  of 
-4 

p  is  about  10  implying  greater  pointing  error  sensitivity  than  for  any  of 
the  cases  shown  in  Fig.  2.  These  results  are  further  confirmed  in  Section 
3  where  computations  are  carried  out  in  the  presence  of  both  noise  and  in¬ 
terference. 

For  further  insight  into  the  behavior  of  the  hybrid  processor  we  ex¬ 
amine  the  output  signal  to  interference  plus  noise  ratio  (SINR)  with  an 
inter ferer  present.  We  will,  however,  now  assume  a  perfect  reference;  i.e. 
F  “1.  The  inverse  covariance  matrix  given  by  (22)  now  becomes 


y  M 
'II 

TjN  +  1 


) 


(34) 


and  the  weight  vector  is  now 


-1  *  u  ,  *  * 

w  »  uM^  £  =  -j  (s  -  k2Ml—  ) 

a 

n 
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where 


k 


2 


YtN  +  1 


05) 


Output  signal  (Sfl) ,  interference  (IQ) ,  and  noise  (Nq)  are  obtained  by 
premultiplying  (34)  by  asj,  Bs*,  and  nT,  respectively.  The  output  SINR 
is  defined  by 


SINR 


E(lS^l) _ 

E( j Iq 1 )  +  E(|No|2) 


06) 


The  signal  output  power  is 


n 


T  *  T  *i 2 
k2  1 


(37) 


If  pointing  is  accurate  so  that  s  -  ^  (37)  becomes 


E(  |S2  I  ) 
o 


siL 


j  *T  1 2 

(NZ  -  2^13^  Sjl 


+  k; 


2 1  *T  | 

1 2d  ^I1 


(38) 


The  noise  output  power  is 


E ( | N  | 2 )  -  ^  E|nT(s* 


,  *  T  *v  |  2 

-  k2-i-i- 


-  [N  -  (2k2  -  k2N)is^s*|2] 

a 

n 


09) 
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and  the  interference  output  power  is 


T  *  T  *,2 


E<II0I  >  *  \  |5  ll^s  -  kjSjSjSjS  | 


-  *4  |82 |(1  -  k2N/ls^s*|2 
a 

n 

Note  from  (35)  if  the  interference  to  noise  ratio  Yj  and/or  the  number  of 
elements  N  is  large  enough  to  make  >>  1,  k2  approaches  1/N.  The  output 
interference  power  as  given  by  (40)  is,  in  this  case,  strongly  suppressed 
through  the  factor  (1  -  k2N). 

With  accurate  pointing  (39)  and  (40)  as  well  as  (38)  involve  the  quantity 
|s^  s^|  ,  the  inner  product  of  the  desired  and  undesired  signal  pointing 
vectors.  If,  for  instance,  the  array  elements  were  strung  out  along  the  x-axis 
with  positions  x^. 


I  T* 


«  N  jkx. (cos0T  -  cos8,) 
l  r  1  I  a 

-  I  e 

i-l 


where  0^.  and  9^  are  the  undesired  and  desired  signal  arrival  angles  relative 
to  array  broadside  and  k  is  the  wavenumber  ("2u/X).  For  a  thinned  array  with 
the  element  positions  x^  far  apart  in  units  of  wavelength  and  with  0^ 
sufficiently  different  from  0^,  the  sum  in  (41)  is  the  vector  sum  of  N  arbitrarily 
oriented  unit  length  phasors.  With  N  large  enough  the  squared  length  tends 
to  be  close  to  N. 

T*  2 

If  we  take  ]s^  s_j|  *  N  and  also  take  k2  to  be  N  we  get  from  (38) 


E(|S2 |) 


(N^  -  2N  +  1) 


2 1  2 1 

4 

a 

n 


(N  -  1)' 


If  kj  had  been  zero,  meaning  therq  were  no  interference  to  suppress, 
2  2 

the  factor  (N  -  1)  in  (42)  would  be  replaced  by  N  .  As  typically 
found,  the  suppression  process  causes  a  reduction  in  signal  output 
power  equal  to  that  of  a  loss  of  one  element. 

Under  the  same  assumptions,  the  noise  output  power  is  from  (39) 


e(!No!2)  -  ^2  (N  -  1)  (43) 

a 

n 


Taking  the  view  that  the  interference  is  totally  suppressed,  the  output 


SINR  is  the  same  as  the  output  SNR  which  is,  using  (42)  and  (43) 


SNR  *  SINR 


(44) 


More  generally,  the  output  SINR  is  given  by  substituting  (37),  (39),  and  (40) 
into  (36). 
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3.  COMPUTATIONAL  RESULTS 

Based  on  Che  results  of  Section  2,  above,  calculations  of  the  steady- 
state  output  slgnal-to-interference  plus  noise  ratio  (SINR)  were  carried 
out  assuming  a  linear  array  exposed  to  a  desired  signal  and  one  or  more 
undeslred  signals.  The  equivalent  weight  vector  defined  by  (7)  is  given  by 

w  “  U[(l  -  Fg)la7TMd  +  +  Mn]_1s*  (45) 

(45)  which  is  a  composite  of  (12),  (13),  and  (15),  was  used  to 
calculate  the  equivalent  weights  and  the  corresponding  mean  square  signal, 
noise,  and  interference  outputs. 

Figure  3(a)  shows  the  the  linear  array  element  layout  used  for  most 

of  the  computations.  Seven  elements  were  placed  non-uniformlv  over  an  interval 

of  10  wavelenghts  and  a  desired  signal  was  assumed  to  be  arriving  broadside 

to  the  array.  Computations  of  SINR  were  made  for  various  deployments  of 

interferers  at  various  power  levels  relative  to  noise  and  signal.  Since  the 

behavior  of  a  sparse  array  with  arbitrarily  selected  element  positions  can  be 

expected  to  depend  on  the  particular  realization  of  those  positions,  a  second 

element  layout  as  shown  in  Figure  3(b)  was  used  to  repeat  a  number  of  these 

computations.  A  sampling  of  Ae  results  obtained  is  presented  below. 

Of  particular  interest  are  the  results  shown  in  Figure  4  for  the  case  of 

a  desired  signal  arriving  broadside  to  the  array  (6d  «  90°)  and  a  single 

interferer  arriving  5°  off  broadside.  Signal,  interference,  and  noise  power 

are  respectively  10  dB,  lOdB,  and  0  dB.  The  SINR  versus  pointing  angle 

is  shown  for  the  case  of  an  array  with  directional  constraint  only  (the 

Applebaum-Chapman  case)  and  for  the  hybrid  array  with  varying  degrees  of 

perfection  in  the  self-generated  reference,  i.e.,  with  various  F  .  With  an 

s 

ideal  reference,  i.e.,  Fg  *  1,  we  see  very  little  sensitivity  to  pointing 
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error  (pointing  error  is  pointing  angle-90o),  while  with  a  directional 
constraint  only,  the  pointing  error  sensitivity  is  extremely  high.  As 
the  reference  generating  loop  tends  toward  the  ideal,  the  sensitivity 
decreases  with  reasonably  tolerable  levels  achieved  (3  dB  loss  in  SINR  for 
0.5°  pointing  error)  with  moderate  residual  reference  loop  attenuation  (5X) 
and  phase  shift  (10°).  Results  not  substantially  dissimilar  were  obtained 
with  the  interference  level  raised  to  20  dB  with  the  alternative  deployment 
of  the  array  elements. 

The  interferer  used  to  generate  Figure  4  being  5°  off  broadside  is 
outside  the  beamwidth  of  the  quiescent  array.  By  quiescent  array  we  mean 
the  array  with  adaption  circuits  inactive.  In  Figure  1  the  quiescent  output 

is  e  .  The  terms  "beamwidth"  and  "mainbeam"  will  be  used  only 

m 

in  connection  with  the  quiescent  array.  The  behavior  of  the  hybrid  array 
with  mainbeam  interference  is  of  interest  and  results  of  computations  illustrat¬ 
ing  its  properties  under  such  a  conditionare  shown  in  Figure  5.  A  single 
interferer  is  assumed  1°  off  broadside  with  all  other  conditions  identical 
to  those  used  to  generate  Figure  4.  Interestingly,  the  sensitivity  to  pointing 
error  of  the  array  with  directional  constraint  only.  Figure  5(c),  is  not  as 
severe  as  for  the  corresponding  case  of  Figure  4  with  the  interferer  off  the 
mainbeam.  The  SINR  in  Figure  5  with  correct  pointing  to  the  broadside  signal 
is  however  about  8  dB  less  than  for  the  case  of  interference  off  the  mainbeam. 
This  confirms  a  result  obtained  by  Bar-Ness  [9]  in  a  comparison  study  of  the 
two  array  processing  methods  which  have  been  fused  here  into  the  Hybrid  Array; 
the  directionally  constrained  array  is  strongly  affected  by  mainbeam  inter¬ 
ference.  Figures  5(a)  -  5(c)  show  the  behavior  of  the  Hybrid  Array  with  ideal 
and  non-ideal  reference  loop.  The  progressive  increase  in  pointing  error 
sensitivity  is  evident  as  the  reference  loop  moves  away  from  the  ideal.  The 
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SINR  with  correct  pointing  is  seen  however  to  be  unchanged  from  that  of  the 
pure  directionally  constrained  case.  Sensitivity  to  interference  level  in 
the  case  of  mainbeam  interference  is  however  more  evident  than  in  the  case 
of  off-mainbeam  interference.  In  Figures  6  and  7  are  shown  results  of  com¬ 
putations  with  interference  level  of  20  dB  and  0  dB  relative  to  the  noise 
level  and  with  all  other  conditions  identical  to  those  of  Figure  5.  With 
20  dB  interference  to  noise  ratio.  Figure  6,  the  SINR  with  perfect  pointing 
is  slightly  below  (about  0.3  dB)  that  obtained  with  a  10  dB  ratio.  With 
0  dB  interference  to  noise  ratio.  Figure  7,  the  SINR  with  perfect  pointing 
is  above  (about  2  dB)  that  obtained  with  a  10  dB  ratio. 

It  should  be  noted  that  the  dips  in  the  SINR  curves,  Figures  5-7,  do 

not  occur  at  values  of  pointing  angle  equal  to  the  respective  interference 

arrival  angles.  Rather,  as  the  level  of  interference  is  raised  relative  to 

that  of  the  desired  signal  level  the  dip  moves  from  a  value  somewhat  above 

88°  to  a  value  close  to  89s,  the  actual  interference  arrival  angle.  This 

behavior  appears  to  be  a  consequence  of  desired  signal  reduction  induced  by 

the  interference  and  the  effect  of  element  noise.  It  should  be  noted  that  the 

T 

multiple  sidelobe  canceller  output,  u  £  in  Figure  1,  contains  both  signal  and 

interference  when  pointing  at  an  angle  below  89°.  Because  signal  and  inter- 

T 

ference  are  angularly  close  to  one  another  the  output  u  j  which  acts  to  cancel 

interference  when  pointing  below  89°  tends  also  to  cancel  signal.  The  signal 

and  interference  powers  in  the  numerator  and  denominator,  respectively,  of  the 

SINR  both  decrease.  Because,  for  Figure  7,  input  interference  and  noise  are 

taken  to  be  of  equal  power,  the  output  interference,  being  partly  suppressed 

when  pointing  away  from  89°,  loses  its  effect  on  the  SINR.  The  net  SINR 

will  therefore  continue  to  decrease  as  the  pointing  angle  moves  further 

below  89°.  The  minimum  actually  observed  occurs  when  the  pointing  angle 

T 

becomes  such  that  the  cancelling  effect  of  u  jr  on  the  desired  signal  goes 
through  a  maximum.  When  the  undesired  signal  is  large  as  in  Figure  6  the 
SINR  immediately  below  89°  is  largely  determined  by  signal  and  interference  - 
not  noise.  Here,  as  the  pointing  angle  moves  below  89°,  though  both 
signal  and  interference  may  be  decreasing,  their  ratio  is  increasing. 


A-19/20 


4. 


CONCLUSION 


The  Hybrid  array  processor  analyzed  here  utilizes  prior  information 
on  expected  direction  of  signal  arrival  and  on  signal  structure.  It  is  a 
fusion  of  techniques  using  a  directional  constraint  and  a  self-generated 
reference  signal.  It  was  here  shown  that  the  sensitivity  to  pointing 
error  typical  of  (zero  order)  directionally  constrained  adaptive  arrays 
is  substantially  reduced,  the  degree  of  reduction  depending  on  the  self¬ 
generating  reference  circuit.  When  the  latter  has  unit  gain  and  zero 
phase-shift  the  sensitivity  is  all  but  eliminated.  With  modest  departures 
from  these  conditions  there  is  still  a  significant  improvement.  In  the 
particular  case  of  a  10  wavelength  linear  array  comprised  of  7  non-uniformly 
spaced  elements  with  a  desired  signal  arriving  broadside  to  the  array  and  an 
interferer  5“  off  broadside,  the  output  SINR  falls  off  about  20  dB  for  a 
pointing  error  of  1/2°  when  no  reference  is  used.  On  the  other  hand,  with 
a  reference  generating  circuit  which  tracks  the  amplitude  of  the  desired 
signal  within  5%  and  without  phase  shift  the  SINR  is  off  by  about  3  dB;  with 
amplitude  5%  off  and  a  phase  shift  of  10°  the  SINR  is  off  by  about  8  dB.  Im¬ 
provements  in  sensitivity  to  pointing  error  are  also  observed  when  the  inter¬ 
ference  is  within  the  beamwidth  of  the  quiescent  array.  However,  such  inter¬ 
ference  is  observed  to  blunt  the  effect  of  pointing  error  even  when  no  re¬ 
ference  is  used.  The  improvement  with  reference  is  therefore  less  dramatic 
in  this  case. 


We  have  alluded  to  the  advantage  expected  in  transient  response  of  this 
scheme  over  one  which  does  not  utilize  the  directional  information.  This 
is  subject  of  another  analysis  and  will  be  reported  separately. 
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APPENDIX  A1  -  EQUIVALENT  WEIGHT  VECTOR 


We  here  extend  the  result  obtained  by  Applebaum  and  Chapman  [1]  for  a 
directionally  constrained  array,  to  the  case  of  the  Hybrid  Array  represented 
in  Figure  1. 

* 

Using  (3),  (4)  and  (10)  in  the  main  text  and  writing  u  ■  Av  we  have 


•ff  £ 

£  =  g  A  v  (£  v  -  £  Av  -  e^) 


(Al-1) 


*  * 

=  g[A  Ms 


★  *  * 

A  MA  y  -  A  v  e 
^  —  r 


] 


where  M  is  the  covariance  matrix  of  inputs;  that  is 


and  it  is  positive  definite  (independent  noise  voltages  are  always  assumed  at  the 
array  elements). 

•k 

To  evaluate  v  e  (1)  is  rewritten 
—  r 

v(t)  =  a(t)sd  +  £(t) 


with  £(t)  representing  all  interference  and  noise  components.  Then  using 
(6),  replacing  v  with  the  signal  part  above  we  get 


e 


od 


A)^ 


(Al-2 ) 


From  (11) 


e 

r 


F  a  (s 
s  — 


*T 


A)  s^ 
— d 


+  n 

r 


so  that 


(Al-3) 


*Equaticn,i  in  the  main  text  wi  t  (  hence  faith  be  t  ecatted  without  {wither 
comment. 


(Al-4) 


£  ft  ft  ★  T 

V  er  -  (a  s^j  +  £  )  [Fgct  (s  -  £  A)^  +  nf] 

-  |a|2Fg  Md(s*  -  ATZ) 

where 

Md  ■  44  (A1-5> 

is  the  covariance  matrix  corresponding  to  the  desired  signal  direction  vector  and 
where  we  have  assumed  that  £,  n^  and  the  desired  signal  are  all  uncorrelated. 
(Since  nf  arises  from  £  some  dependence  is  apt  to  exist;  it  is  here  assumed 
negligible).  Substituting  (Al-4)  into  (Al-1)  we  get 

£  -  g(A*Ms*  -  A*MAT£  -  |a2jFgA*Md(s*  -  a\)  ]  (Al-6) 

-  gA*(M  -  To^ly^Ks*  -  AT2_) 

T 

Premultiplying  by  A  and  using  (8)  we  get 

s  *  -  w  =  gATA*[M  -  |a2|FgMd]w  (Al-7) 

This  result,  from  which  the  weight  vector  w  can  be  obtained  is  in  a  form 

similar  to  that  in  Applebaum  and  Chapman  [1,  eq.  (4)],  the  latter  being 

obtained  from  (Al-7)  with  F  =  0. 

s 

For  the  case  g  »  1  we  can  also  get  from  (Al-6) 
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A  [M  -  |a2  |  F^M^lw  () 


(Al-8) 


Since  A^  =  0,  meaning  that  the  (N  -  1)  linearly  independent  N-dimensional 
row  vectors  of  A  are  orthogonal  to  £,  the  solution  to  (Al-8)  must  satisfy 


[M 


* 

us 


or 


w  =  u[M  -  |a2|FgMd]_1s* 


where  u  is  an  appropriate  constant.  From  (9)  we  have 


T  T  i  2 i  —1  * 

sAw  =  N  =  us  [M  -  jot  | FgM  ]  s 


(Al-9) 


(Al-10) 


so  that 


u 


sT[M 


-1  * 
W  ± 


(Al-11) 


and 


-1  * 

-1  * 

-  =  uMi  -  =  TTT 

£  s_ 


where 


M.  ■  M  -  a  F  M. 
1  s  d 


can  be  shown  to  be  non-singular. 


(Al-12) 


(Al-13) 
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APPENDIX  A2  -  INVERSION  OF  THE  COVARIANCE  MATRIX 


We  obtain  the  inverse  of  the  covariance  matrix,  M,  of  the  array 
element  output  vector  as  given  by  (1)  when  a  single  interferer  denoted 
g(t)^j  is  present.  M  may  be  written  (see  (15)- (18)) 

M  =  o*[I  +  y.M  +  YtMt]  =  .  (A2-1) 

nn  da  II  no 

*  T  *  T 

where  =  s^s^,  =  s^s^.,  IR  is  an  n  by  n  identity  matrix,  and  Yd  and  y 

are  SNR's  of  desired  signal  and  interference,  respectively.  M,  and  MT 

d  i 

are  positive  semi-definite  so  that  M  is  positive  definite  and  therefore 
invertible. 

For  any  two  matrices  A  and  B  of  dimension  (n  x  r)  and  (r  x  n) ,  respectively, 
if  the  inverse  of  ([  +  AB) exists,  then  [10] 

(I  +  AB)"1  -  ln  -  A(Ir  +  BA)_1B  (A2-2) 

I  is  an  r  x  r  identity  matrix. 

Eq.  (A2-1)  may  be  rewritten 

r  * 

M  =  I  +  Y  jSj 
o  n  d-d 

L 


•  T1 

i 


(A2-3) 


T 

% 


and  its  inverse  may  be  written 
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where  A  ■ 


(Nyd  +  1)(NYi  +  1)  -  VdVIl«j«1l  •  Thus 


1,  *  *, 

r  ,  t  *  "i 

T 

ln  "  A*  Yd— d  YA] 

Nyd  +  1  ^iV-I 

T  * 

T 

_Yd~ I— d  ^d  +  XJ 

L%J 

+  l)»d  -  Vl<MdMI  +  MtMd)  *  ,I<Nyd  +  DKj) 


Equation  (22)  in  the  text  is  obtained  from  this  result  by  replacing  by 

2 

y^ (1  -  Fg)  and  multiplying  by  l/an  to  give  the  inverse  of  the  effective 


I 


covariance  matrix,  M, . 


DETAIL  OF  MULTIPLE 
SIOELOBE  CANCELLER 


FIGURE  2.  POINTING  ERROR  SENSITIVITY 
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7-element  10X  array;  Signal,  Interference  power  10  dB 
relative  to  Noise  power,  Interferer  at  85° 

(a)  F  =  1,  (b)  F  =  0.95,  (c)  Fc  =  0. 95exp(j 10°) , 


7-element  10X  array;  Signal,  Interference  power  -  10  dB 
relative  to  Noise  power,  Interference  at  89°. 

(a)  F  “  1,  (b)  F  -  0.95exp(-jl0°), 
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7-element  10X  array;  Signal,  Interference  powers  10  dB 
and  0  dB  relative  to  Noise  power.  Interference  at  89°. 

(a)  F  -  1,  (b)  F  =  0.95exp(-j 10°),  (c)  F  -  0. 
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List  of  Symbols  -  Appendix  A 


A 

E[  ] 

e 

ra 

e 

o 

eod 

eG(sig) 


g 

I 

I 

c 

J 

k 

k. 


1 

-1 


‘U 


N 

N 


transformation  matrix 
taking  expectation  of  [  ] 
output  of  the  beamformer 
array  output 

signal  component  of  the  array  output 
output  desired  signal  of  a  beamformer 
reference  signal 

operator  representing  the  effect  of  the  bandpass  filter  in  the 
reference  loop  on  the  desired  output  signal  component 

gain  in  the  weight  setting  loops 

identity  matrix 

output  interference 

number  of  interferers 

wave  number 

constant  related  to  reference  quality 

constant  related  to  interference  to  noise  ratio 

linear  dimension  of  the  array 

covariance  matrix  of  array  element  outputs 

equivalent  covariance  matrix 

inverse  of  M 

covariance  matrix  corresponding  to  the  desired  signal  direction  vector 

covariance  matrix  corresponding  to  the  jth  interferer  direction  vector 

covariance  matrix  corresponding  to  noise 

normalized  covariance  matrix 

number  of  array  elements 

output  noise 
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(Continued) 


List  of  Symbols  -  Appendix  A 
(Continued) 


V  1 


noise  component  of  the  reference  signal 

complex  noise  envelope  vector 

interference  and  noise  components  of  v(t) 

normal  output  signal  to  noise  ratio  -  SNR/ct.N 

a 

real  part  of  [  ] 
steering  vector 
complex  conjugate  of  s^ 

arrival  phase  vector  of  the  desired  signal 

the  iC^  unit  amplitude  component  of  S, 

— d 

arrival  phase  vector  of  the  j1"*1  interferer 
output  signal 
transpose  of  £ 

signal  to  interference  plus  noise  ratio 
signal  to  noise  ratio 

input  vector  of  the  multiple  sidelobe  canceller 


v.(t) 


6  (t) 


output  of  the  i  array  element 


output  vector  of  array  elements 

equivalent  weight  vector 

iC^  array  element  position 

weight  vector  on  the  sidelobe  canceller 

desired  signal  waveform 

power  of  a(t) 

power  of  B^(t) 

jC^  interferer  waveform 

signal  to  noise  ratio  at  array  elements 


A- 38 


(Continued) 


List  of  Symbols  -  Appendix  A 
(Continued) 

interference  to  noise  ratio  at  array  elements 
determinant 

measure  of  pointing  error 
arrival  angle  of  the  desired  signal 
interferer  arrival  angle 
wavelength 
a  constant 

a  constant  related  to  reference  quality 
noise  power  at  array  elements 


carrier  frequency 
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TRANSIENT  RESPONSE  OF  THE  HYBRID  ARRAY 
by 

Chien-Chung  Yeh 

1.  Analysis 

We  examine  here  the  transient  behavior  of  the  Hybrid  array.  The 
circuit  is  that  shown  in  Fig.  1  of  Appendix  A  and  the  notation  is  also 
that  of  Appendix  A.  The  input  vector  is  written 

V(t)  =  a(t)S^  +  BCOSj  +  N(t)  (1) 

where  a(t)  and  B(t)  are  complex  envelopes  of  desired  and  undesired 
signals  respectively,  and  S^.  are  their  direction  of  arrival  vectors, 
and  N(t)  is  the  random  noise  vector. 

The  output  of  the  beamformer  is 

em(t)  =  VT(t)S*  -  S*TV(t> 

with  ;>  being  the  steering  vector.  We  also  have 

U(t)  =  A  V(t)  (2) 

where  A  S  =  0.  The  sidelobe  canceller  generates  y^(t)  »  Y^U(t)  “  Y^A  V(t) 
so  that 


6o(t)  "  6m  '  yl(t) 

-  (S*  -  ATY)TV(t)  (3) 

Define 

eQ'  (t)  -  eQ(t)g(t)  (4) 

where  g(t)  is  the  spreading  code.  eQ'  (t)  is  the  array  output  with  the 
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desired  signal  despread.  The  output:,  eQut(t),  £be  bandpass  filter 
assuming  it  to  be  a  single  pole  with  time  constant  is  given  by 
de  (t) 

t.  — -  +  ert  ft)  -  e  *(t) 

1  dt  out  o 


(5) 


The  reference  signal,  er(t),  is  obtained  from  eQut(t)  by  respreading 
using  the  code  g(t); 

er(t)  =  eout(t)g(t) 

The  residue  feedback  is 


ef (t)  *  eQ(t)-er(t) 

=  [S*-ATY(t)]TV(t)  -  e  .(t)g(t) 

—  —  —  out 

Tor  the  weight  control  loop,  we  have 

dY(t)  * 

t 2  +  Y(t)  -  G  e(t)  0  ft) 

Substituting  (3)  and  (4)  into  (5),  we  get 

de  (t)  .  _  t 

Tl  - -  +  eQut(t)  -  [S  -  A1Y(t)]iV(t)g(t) 


(6) 


(7) 


(8) 


We  examine  the  transient  response  by  assuming  that  the  signal  component  in 
V(t)  is  a  step  function.  Let 

a(t)  =  ag(t),  t  i  0 

0  t  <  0 

The  foregoing  may  be  viewed  as  the  waveform  seen  with  the  first  informa¬ 
tion  digit  received;  the  amplitude  is  a  and  the  signal  is  modulated  by 
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the  spreading  code  g(t).  Taking  expectations  on  both  sides  of  (8) 


de  ..(t) 
out 

dt 


+  e  (t) 
out 


IS  -  ATY(t)]T  V(t)g(t) 


(9) 


It  is  reasonable  to  assume  that  the  interference  and  noise  are 
independent  of  the  spreading  code.  Also  in  an  adaptive  array,  the  weights 
generally  vary  much  more  slowly  than  the  input,  i.e.  the  control  loop 
bandwidth  is  much  smaller  than  the  bandwidth  of  V(t).  See  e.g.  [1],  where 
a  similar  argument  is  used  on  the  Independence  of  weight  fluctuation  and 
input.  We  therefore  assume  that  to  an  adequate  approximation  the  weights 
may  be  treated  as  being  independent  of  V(t).  (9)  then  becomes 


T1 


de  (t) 
out 


dt 


+  eout(t)  =  [S*-ATY(t)JT 
-  [S*-ATY(t)]T  otS^ 


V(t)g(t) 


(10) 


Next,  substituting  (2)  and  (6)  into  (7),  we  have 

dY(t)  * 

t2  +  Y(t)  -  G  e(t)V  (t) 

-  G( [ S *- ATY ( t ) ] TV ( t ) -eQut ( t ) g ( t ) } A *V* ( t ) 

-  G{A*V* (t)VT (t) (S*-ATY(t) l~A*V*(t)g(t)eout (t)}  (11) 


Taking  expectations  on  both  sides  of  (11)  and  using  the  same  assumptions 


on  Independence  we  get 


T 


2 


dY(t)  _ 


G(AV(t)VT(t)[S*-ATY(t)] 


-  A‘(t)g(t)eout(t)} 


(12) 
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Since  e  (t)  is  the  output  of  a  narrowband  L.P.F.  or  B.  P.  F. ,  as  com- 
out 

pared  to  the  spreading  code  chip  rate,  it  should  be  admissible  to  treat 
it  as  approximately  independent  of  the  product  of  V(t)  and  g(t).  That  is, 

it  ft  ft 

we  can  assume  V  (t)g(t)  «  oS^  +■  N^(t)  and  eout(t)  *  oS^  +  Hg(t)  with 
the  fluctuation  components  N^t)  and  N2(t)  essentially  independent.  Thus 

(12)  becomes 

_  *  *  T_ 
t2  — —  +  Y(t)  -  G(A  M[S  -AAY(t)] 

"  A*  I  (t)g(t)  •  eout<t> 

-  G(A*M(S*-ATY(t>]  -  A*aS,*e  ft)}  (13) 

—  —  — a  out 

where 

M  -  V*(t)VT(t) 

■  *2s,V  +  b2st*s_t  +  a  2i  (14) 

— d  —I  —I  n 

with 

a2  -  |«2(t) | 

B2  -  | S 2 ( t )  | 

2  T 

a  X  -  N  (t)N  (t) 
n  —  — 

(13)  can  be  arranged  as 

dY(t)  *  T  _  *  *  *  * _ 

T.  — —  =  -(I  +  GA  MA  )Y(t)  +  G[A  MS  -  A  aS.  e  „(t)]  (15) 

i  at  —  —a  out 

(10)  and  (15)  are  coupled  sets  of  linear  equations  which  are  to  be  solved. 

As  they  stand  they  are  tedious  to  solve.  However  with  the  following 
approximation  a  manageable  result  is  obtained.  Assume  that  is  small. 
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is  large  and  G  is  moderate  so  that  eout(t)  varies  much  faster  than 

Y(t).  Then  in  (10)  we  may  treat  Y(t)  as  a  constant  and  get  the  solution 

of  eout(t>  as  _t 

e  At)  -  a(S*  -  A'TCOlVa-e  Tl)  (16) 

out  —  ~ ■  “U 


Substituting  (16)  into  (15), 


dY(t) 
T2  dt 


t 


*  2  *  T  T  .  * 

+  GA  [M  -  aS,  S,  (1  -  e  1)]S 
-tl  — a  — 


(17) 


In  the  following  we  look  at  the  physical  meaning  and  justification 
of  (17).  First,  if  the  reference  loop  were  infinitely  responsive,  would 

_  *  T_  t 

have  to  be  near  zero  and  (10)  would  reduce  to  e  It)  *  a[S  -A  Y(t)]  S,. 

out  —  —  — d 

This  substituted  into  (15)  would  give 

^7(t)  it  3  It  T  2  T  _ 

T2~dt~  *  ~[I  +  GA  (6  Sj  Sj  +  an  I) A  )Y(t) 

*  2  *  T  2  * 

+GA  (B  ST  ST  +  a  I)S  (18) 

— I  - 1  n  — 

this  result  is  embodied  in  (17);  note  that  if  t  is  large  enough  so  that 
e  t/,Tl  -  0  (17)  becomes  identical  to  (18).  Next,  we  know  that  without 
the  reference  loop,  the  weight  equation  would  be 


dI<c>  *  T  -  *  * 

x„— —  -  -(I  +  GA  MA  )Y(t)  +  GA  MS 

_JL 

With  t  =  0,  e  T1  =  1,  so  that  (17)  turns  out  to  be 


(19) 


dI<c)  *  T  -  *  * 

x„— - —  =  -(I  +  GA  MA  )Y(t)  +  GA  MS 
l  dt  —  — 

It  is  the  same  as  (19). 
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To  summarize,  when  t=0,  the  output  of  the  reference  loop  is  very  small 
and  the  weight  equation  is  similar  to  that  without  the  reference  loop. 

For  t  large  enough,  the  output  of  the  reference  loop  reaches  its  steady 
state  value  and  the  array  works  as  if  we  have  a  perfect  reference  loop. 

(17)  may  therefore  turn  out  to  a  good  approximation. 

It  should  be  pointed  out  that  an  information  bearing  waveform  a(t) 
will  be  comprised  of  a  random  sequence  of  positive  and  negative  steps. 

Even  after  the  initial  transient  there  will  be  signal  related  feedback  in 
resulting  from  inability  of  the  reference  loop  to  instantaneously  follow 
signal  changes.  These  will  affect  the  weight  control  loop  even  though  the 
interference  input  to  the  weight  control  loop  and  the  signal  residue  on 


The  solution  of  (17)  for  the  case  of  a  two  element  array  is  carried 
out  in  the  following.  For  a  two  element  array,  Y(t)  is  a  scalar.  There¬ 
fore  (17)  is  a  first  order  linear  differential  equation.  We  divide  both 
sides  of  (17)  by  and  rearrange  it  to  give 


+  {i_  +  £_  fA*(g2s  *s  T  +  02I)AT  +  A*a2S,*S,TATe  T1  ])y(t) 
dt  -I  -I  n  -d-d 


‘2  k2 


J-  A*[(82Si*£iT  +  on2I)  +  a2S^*S^T  e  Tl]£* 
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From  the  theory  of  ordinary  differential  equations,  the  solution  of 


is 


where 


and  C 


and 


Using 


+  f (t)y (t)  «  r (t) 


(22) 


ft 


y(t)  *  e“h^  [ 


2h(t')r(t')dt'  +  C] 


(23) 


h(t) 


f(t')dt' 


(24) 


is  a  constant.  (21)  is  exactly  in  the  same  form  as  (22)  with 


1  g  .  *  2  *  T  2  2*TTt 

f(t)  -  ~  +  -p  [A  (8  Sj  Sj  +  on  I)  +  a  ^  Sj  A  e  1] 


(25) 


G  2*T  2  2*Tt 

r (t)  -  f-  [(8ZSX  Sx  +  onZ  I)  +  a  e  Tl] 


(26) 


(24)  and  (25) ,  we  get 


h(t) 


(7-  +  f-  [A*(62S  *S  T  +  ojl) AT  +  A*aVVTAT  e  Tl]}dt' 
t  2  t  2  — 1  — 1  n  —a  —a 


,1  G  *  2  *  T  2T 

I—  +  —  A  (6  S_  S/  +  o^I)A  ]t 
1 2  *  2  a  1  n 


+  aVs^S^  AT  tx  (1  -  e  Tl) 


(27) 
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The  constant  C  in  (23)  is  determined  by  the  initial  condition.  If  we 
let  y(0)  ■  0,  then 

C  -  y(0)  -  0 

Therefore,  the  solution  of  (21)  is 

t 

eh(t’)r(t,)dt  (28) 

3 

with  h(t)  and  r(t)  given  by  (27)  and  (26),  respectively. 

2.  SIMULATION  RESULTS 

The  behavior  of  a  two  element  array  with  X/2  spacing  as  shown  in 
Figure  1,  was  simulated.  The  desired  signal  was  assumed  to  be  spread  in 
spectrum  with  an  information  sequence  of  100  bps  generated  by  a  seven  digit 
shift  register  and  a  code  sequence  of  10K  bps  generated  by  a  10  digit 
shift  register.  The  desired  signal  arrival  angle  was  assumed  to  be  at 
broadside,  i.e.  6d  *  90°  and  the  interferer  was  assumed  to  be  a  pure 
carrier  arriving  from,  0^.  =  45°,  in  phase  with  one  of  the  two  binary 
signal  alternatives. 

In  the  first  part,  we  assume  that  there  is  no  noise  and  the  SIR  is 

measured  at  eQ.  The  SIR  is  defined  as  the  ratio  of  squares  magnitudes  of 

the  samples  of  the  desired  signal  and  interference  at  each  time.  The 

signal  power  and  interference  power  are  assumed  the  same  and  are  chosen  to 

be  100.  The  steering  angle,  6,  ranges  from  90°  to  70°.  G  was  chosen  such 
2 

that  G6  >>  1.  was  chosen  such  that  the  bandwidth  of  the  reference  loop 
is  about  the  same  or  several  times  the  information  bandwidth  of  the  desired 
signal.  Tj  was  chosen  so  that  the  control  loop  bandwidth  is  much  smaller 
than  the  information  bandwidth.  Simulation  was  done  at  baseband. 
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Figure  2  is  the  result  of  simulation  with  9  *  90°,  G  -  10,  and 


■  4;  t ^  is  a  parameter.  Two  values  are  chosen  for  t^: 


200it  and  IOOOtt 


For  the  first  one,  the  bandwidth  of  the  reference  loop  is  the  same  as  the 
information  bandwidth.  For  the  second,  the  reference  loop  has  bandwidth 
five  times  the  information  bandwidth.  The  results  show  that  the  response 
with  -  200tt  fluctuates  more  than  with  but  both  have  the 

same  convergence  speed.  At  t“0;  the  initial  value  of  the  controlled 
weight  y  was  set  to  zero,  and  the  output  SIR  is  about  7  dB.  When  t 
increasing  to  half  the  information  bit  duration,  the  SIR  increases  to 
about  40  dB.  At  the  end  of  the  first  information  bit,  the  SIR  is  about  60 
dB.  However,  as  the  polarity  of  the  information  bit  changes,  the  SIR 
decreases  drastically  and  then  starts  to  increase,  again. 

An  explanation  of  Fig.  2  is  given  in  the  following.  The  weight  y 
is  controlled  by  the  correlation  of  and  U.  With  correct  pointing,  U 

contains  only  interference.  Let  S  and  I  denote  desired  signal  and  inter¬ 
ference  respectively.  U  can  be  written 


U  =  I 


The  feedback  is 


(29) 


ec  =  e  -  e 
for 

e  can  be  written 
o 

e  =*  S  +  I 
o  o  o 

e^  can  be  approximated  by 


e  *  S  +  N 
r  r  r 


(30) 


(31) 


(32) 


where  is  a  noise  term  resulting  from  interference  not  totally 
eliminated  by  the  reference  loop  filter.  Substituting  (31)  and  (32)  into 
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(30),  we  get 


(S  -  S  )  +  I  -  N 
o  r  o  r 


Using  (29)  and  (33),  we  have 


<■/  ■  <so  -  V  h  *  Vu*  -  Vu* 


*  *  * 
S,I„  +  I  I„  -  N  I„ 
f  U  o  U  r  U 


where  S,  -  S  -  S  . 
for 

* 

In  (34)  the  term  I  I„  drives  y  to  cancel  the  interference,  and  the 
o  u 

•k  * 

terms  S^I^  and  N{Iy  make  the  weight  fluctuate.  The  fluctuation  of  the 
weight  increases  as  or  G  increases. 

The  desired  signal  component  of  the  feedback  decreases  during  the 
transient  as  the  bandwidth  of  the  reference  loop  increases.  That  is 


why  the  response  for  t 


is  more  noisy  than  that  for  t 


1  IOOOtt 


However,  increasing  the  bandwidth  of  the  reference  loop  also  increases  the 
noise  component,  Nr,  in  the  feedback.  There  is  a  trade-off  between  de¬ 
creasing  and  increasing  Nf. 

The  drastic  drop  in  SIR  as  the  polarity  of  the  information  bit 

changes  is  caused  by  the  delay  of  the  reference  loop.  The  waveforms  of 

S  ,  S  and  S.  during  the  transition  time  are  shown  below, 
or  f  ° 


When  the  polarity  changes,  there  is  a  spike  in  S^.  In  the  short  interval 

* 

of  its  duration,  the  term  S^I^  drives  the  weight  away  from  its  ideal 
value  and  the  SIR  decreases  drastically. 

It  should  be  pointed  out  that  if  the  desired  signal  power  is 
much  smaller  than  the  interference,  the  weights  will  not  be  so  effected  by 
the  signal  when  its  polarity  changes. 

The  convergence  speed  is  controlled  by  G  and  t^.  Increasing  G  or 
decreasing  will  increase  the  convergence  speed.  In  Fig.  3  convergence 
rates  for  G  =  30  and  G  “  10  are  compared .  Since  the  fluctuation  increases 
as  G  increases,  was  chosen  to  be  to  reduce  the  fluctuation. 

The  results  show  that  the  response  for  G  *  30  is  about  three  times  faster 
than  that  for  G  ■  10  as  expected,  and  reaches  steady  state  in  about  one 
half  the  duration  of  an  information  bit. 

We  kept  G  -  30,  «  1QQQrr  and  t 2  -  4,  but  changed  the  pointing  angle 

to  see  the  effect  of  pointing  error.  Four  steering  angles  were  chosen, 
namely,  0  =  90°,  which  is  the  correct  angle,  89®,  85®  and  70°.  Results 
are  shown  in  Fig.  4.  With  5®  of  pointing  error,  the  array  is  able  to 
work  well  and  reaches  steady  state  in  about  half  an  information  bit, 
though  there  is  some  degration  at  steady  state.  However,  for  0  -  70®, 
the  convergence  speed  is  relatively  slow  and  can  not  reach  steady  state 
within  the  first  information  bit. 

The  slow  rise  in  SIR  when  pointing  at  70®  requires  explanation.  For 
a  large  pointing  error,  the  desired  signal  component  in  U  is  not  small.  This 
signal  component  will  correlate  with  the  signal  component  in  the  feedback 
making  the  signal  component  into  the  reference  generator  time  varying.  The 
apparent  effect  is  to  significantly  increase  the  convergence  time.  This 
indicates  that  knowledge  of  the  direction  of  the  desired  signal,  or  a  reason¬ 
able  approximation  to  it,  can  be  used  to  significantly  reduce  weight  con¬ 
vergence  time. 
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Simulation  of  a  two  element  Applebaum  array  was  also  carried  out  by 
setting  e^  to  zero.  Results  show  that  the  fluctuation  in  response  is 
much  larger  than  that  of  the  Hybrid  array.  The  reason  for  this  is  that 
the  signal  component  in  the  feedback  is  very  large  in  the  case  of 
the  Applebaum  array  while  the  signal  component  of  the  feedback  is  largely 
eliminated  by  the  reference  loop  in  the  case  of  the  Hybrid  array.  We 
have  seen  that  the  fluctuation  increases  as  increases.  In  order  to 
reduce  the  fluctuation,  we  have  to  decrease  G  or  increase  This* 

however,  will  make  the  convergence  speed  of  the  Applebaum  array  much 
slower  than  that  of  the  Hybrid  array. 

We  next  added  random  noise  to  the  simulation  described  above.  In 
addition  we  took  the  output  point  to  be  the  point  identified  as  eQut  in 
Fig.  1  and  determined  the  SINR  here.  First  the  array  element  signal 
power,  interference  power,  and  noise  power  were  taken  to  be  1,  100,  and  1, 
respectively.  Since  the  spreading  ratio  is  100  to  1,  without  adaption  the 
output  signal  and  interference  have  about  the  same  power  and  are  20  dB 
higher  than  the  noise.  With  G  ■  10,  -  JoStT  and  T2  “  ^our 

pointing  angles,  90°,  89°,  85°  and  70°  were  chosen  to  carry  out  the 
simulations.  Results  are  shown  in  Fig.  5.  It  can  be  seen  that  with  5° 
of  pointing  error  the  loss  is  not  much.  Within  half  an  information  bit, 
the  array  can  null  the  Interference  such  that  its  power  is  much  smaller 
than  noise  and  the  SINR  depends  on  the  signal  to  noise  ratio.  Therefore 
the  SINR  is  about  20  dB  which  is  the  signal  to  noise  ratio  at  the  input. 
However,  for  9  **  70°,  the  convergence  speed  is  relative  slow.  But  the 
array  catches  up  quickly  after  the  first  information  bit. 
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Ic  is  worth  noting  that  since  the  signal  power  at  the  array  elements 


is  chosen  to  be  much  less  than  the  interference,  the  weight  is  not  much 
effected  by  the  information  bit  polarity  reversals.  There  is  a  drop  in 
instantaneous  SINR  as  the  polarity  changes  since  the  signal  falls  to 
zero  at  this  time  but  the  weight  is  apparently  not  seriously  affected. 

We  then  simulated  the  case  when  signal  power,  interference  power,  and 
noise  power  are  100,  100  and  1,  respectively.  The  values  of  G,  t^,  t 
and  pointing  angles  used  are  the  same  as  those' in  Fig.  5.  Simulation  re¬ 
sults  are  shown  in  Fig.  6.  We  see  that  with  5°  of  pointing  error,  the 
array  can  achieve  40  dB  SINR,  which  is  determined  by  the  signal  to  noise 
ratio  at  the  input,  within  half  an  information  bit.  For  pointing  at 
9  =  70°,  the  response  is  slow  and  is  not  able  to  null  the  interference 
below  noise  within  one  information  bit.  As  the  polarity  of  the  informa¬ 
tion  bit  changes  the  weight  is  driven  away  from  its  original  value  be¬ 
cause  the  signal  has  the  same  power  as  the  interference  at  the  array 
input. 

To  justify  the  assumptions  made  in  analyzing  the  transient  response, 

the  controlled  weight  obtained  from  (28)  of  the  work  on  transient  response 

is  compared  with  the  weight  obtained  from  the  simulations.  Only  the 

imaginary  parts  are  compared;  the  real  parts  are  very  close  to  zero. 

Results  are  shown  in  Fig. (7)  and  (8).  The  values  of  G,  t  ,  and  t2  are  the 

same  for  both  figures,  which  are  G  =  10,  =  ~2tfQv  and  T2  =  l*‘  The 

signal  and  interference  arrival  angles  are,  0  =  90°  and  9r  =  45°  res- 

d  J- 

pectively.  Fig.  7  shows  the  results  of  the  noise  free  case  with  pointing  angle 
at  9  *  89°.  Curve  a  is  the  theoretical  results  from  (28)  and  curve  b  is 
the  simulation  results.  These  two  curves  match  very  well.  Curve  C  is 
the  theoretical  results  when  we  assume  that  * 


0,  i.e.  the  reference 


OOP  U  ideal.  Figure  8  »ho„.  t„e  t..ult.  for  the  c<8e  ^  ^  ^ 
8l8“1  -  ««  P.«,  ere  100,  100,  «„d  ! 

"■  th«  “»  *-“««»  are  ba.eO  0„  reason 

able  approximations. 
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FIGURE  1 


Two  Element  Hybrid  Array 
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FIGURE  2  Transient  Response  -  Hybrid  Array  Noiseless  Case 
Signal  Power  20dB  Interference  Power  20dB 
0  *  90°,  0T  =  45°,  0  *  90°,  G  -  10,  t-  «  4 


FIGURE  3  Transient  Response  -  Hybrid  Array  Noiseless  Case 
Signal  Power  20dB  Interference  Power  20dB 
8.  -  90°.  8.  «  45°.  6  -  90°.  t.  -  rr~~.  T  -  4 


90°,  eT 


J 

I 

I 

» 

i 


FIGURE  4  Transient  Response  -  Hybrid  Array  Noiseless  Case 


Signal  Power  20dB  Interference  Power  20dB 
0d  -  90°,  9t  -  45°,  0  -  30,  Tj  -  4 
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FIGURE  5  Transient  Response  -  Hybrid  Array 
Signal  Power  OdB  Interference  Power  20dB 


Noise  Power  OdB 
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FIGURE  6  Transient  Response  -  Hybrid  Array 

Signal  Power  20dB  Interference  Power  20dB 
Noise  Power  OdB 
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FIGURE  7  Weight  Transient  -  Hybrid  Array  Noiseless  Case 

Signal  Power  20dB  Interference  Power  20dB 

ed  -  90°,  et  -  45°,  9  -  89,  G  -  10,  .  ji-,  ^  .  4 
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FIGURE  8  Weight  Transient  -  Hybrid  Array 

Signal  Power  20dB  Interference  Power  20dB 
Noise  Power  OdB 
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transformation  matrix 
band  pass  filter 
beamformer  output 
array  output 

array  output  with  desired  signal  despread 
output  of  the  band  pass  filter 
reference  signal 

gain  of  the  weight  control  loops 
spreading  code 

time  function  as  defined  in  text 
identity  matrix 
interference  component  of  eo 
interference  component  of  U 
low  pass  filter 
covariance  matrix 
noise  component  of  e^ 

input  noise  vector 
noise  vector 
noise  vector 

time  function  as  defined  in  text 
steering  vector 

arrival  phase  vector  of  the  desired  signal 

arrival  phase  vector  of  the  interferer 
desired  signal  component  of  eQ 

desired  signal  component  of  e^ 
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(continued) 


Y 

yx(t) 

a 

a2 

a(t) 

B2 

B(t) 
ef  (t) 
e 


List  of  Symbols  -  Appendix  B 


(Continued) 


complex  conjugate  of  S_ 

input  vector  to  the  sidelobe  canceller 

vector  of  array  element  outputs 

transpose  of  V 

sidelobe  canceller  weights 

sidelobe  canceller  output 

amplitude  of  the  desired  signal 

desired  signal  power  at  array  elements 

desired  signal  waveform 

interference  power  at  array  elements 

interferer  waveform 

residue  feedback 

steering  angle 

desired  signal  arrival  angle 

interferer  arrival  angle 

noise  power  at  array  elements 

time  constant  of  the  reference  loop 

time  constant  of  the  weight  control  loops 
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Appendix  C 


EFFECTS  OF  RANDOM  AMPLITUDE  AND  STEERING  PHASE  ERRORS 
ON  THE  BEHAVIOR  OF  THE  HYBRID  ARRAY 

Chien  -  Chung  Yeh  and  Fred  Haber 

The  adaptive  array  processor  proposed  by  Applebaum  f 1]  uses  prior 
information  on  the  direction  of  signal  arrival,  through  a  steering  vector. 

A  different  system  realization  which  constrains  the  processor  to  maintain 
constant  gain  of  the  main  beam  is  later  described  by  Applebaum  and 
Chapman  [2].  Ideally,  both  systems  achieve  maximum  signal  to  interference 
plus  noise  ratio.  For  application  to  communication,  both  systems  are  very 
sensitive  to  the  accuracy  of  the  steering  vector  [3,-i],  To  overcome  this 
shortcoming,  a  Hybrid  array  processing  technique  has  been  proposed  for  point- 
to-point  spread  spectrum  communication  [5).  It  has  been  shown  that  the  hybrid 
processor  is  much  less  sensitive  to  error  of  the  expected  direction  of  signal 
arrival. 

In  addition  to  the  error  caused  by  imperfectly  information  of  the 
direction  of  signal  arrival,  the  steering  vector  might  be  subject  to  random 
phase  errors  arising  from  uncertainties  in  element  positions  or  as  a  consequence 
of  quantization.  In  addition  the  system  might  introduce  random  amplitude 
errors  on  the  array  element  outputs  prior  to  processing,  arising  from  circuit 
disimilarities.  An  analysis  of  effect  of  random  errors  in  the  steering 
vector  in  the  Applebaum  array  has  been  given  by  Compton  [6].  Here  we  study 
the  effect  of  phase  and  amplitude  perturbations  on  the  performance  of  Hybrid 
array  for  both  the  gain  constrained  and  the  non-const  rained  processors.  The 
two  versions  of  Hybrid  array  based  on  the  Applebaum  arrays  with  and  without 
gain  constraint  are  shown  in  Figures  1  and  2,  respectively.  The  analysis  assumes 


1 


narrow  band  desired  and  undesired  signals  and  a  reference  generating 
circuit  which  has  achieved  perfect  spread  spectrum  chip  synchronism  with 
respect  to  the  desired  signal.  The  analyses  and  computational  results 
are  presented  below. 

Analysis: 

Consider  an  N-element  linear  array  with  random  gain  error  c  ,  n*l,  2,  ..  N, 

n 

on  each  element.  The  c  's  are  all  assumed  to  be  random  variables  with  zero 

n 

2 

mean  and  variance  o  .  The  output  signal  of  the  array  elements  may  be 
c 

represented  by  the  complex  vector 

V'(t)  -  CV(t)  (1) 

where  C  is  an  N  x  N  diagonal  matrix  with  elements  1+c^,  l+t^*  •••  ^+cu  on 
the  diagonal,  and  V(t)  is  the  output  vector  of  the  array  elements  without 
random  gain  error.  V(t)  is  represented  by 

J 

V(t )  -  0(0^+  ^  Bj(t)S  +  N(t)  (2) 

where  S,  is  the  arrival  phase  vector  of  the  desired  signal  comprised  of  N 

unit  amplitude  components,  a(t)  is  its  complex  envelope;  is  the  arrival 

phase  vector  of  the  jth  interference  signal,  6j(t)  is  its  complex  envelope; 

and  fJ(t)  is  the  complex  noise  envelope  vector,  the  components  of  which  are 

2 

assumed  independent  with  power  o^  .  If,  for  instance,  the  elements  are  arranged 

along  the  X  axis  at  positions  x  ,  n*l,  2,  ..  N,  the  ith  component  of  the 

n 

arrival  phase  vector  of  a  desired  signal  arriving  as  a  plane  wave  at  an  angle 
6d  relative  to  the  X-axis  Is  s^  ■  exp(  jkx^ose^) . 
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A  steering  vector  intended  for  signals  with  arrival  angle  6  but  with 


random  phase  error  can  be  written 


jkx^cosB  j  <t>1 

r  ■— i 

e  *e 

si 

jkx^cosQ  j$2 

e  .  »e 

= 

S2 

jkv  cose  j<t>N 

e  *e 

N 

(3) 


where  4>  ,  n=L  2,  . .  .N  are  random  errors.  We  will  assume  the  latter  to 
n  2 

be  i.i.d.  random  variables  N(0,o  ).  From  (3),  the  mean  value  of  S  is 


S  =*  eJ 


j-  j  kx^cos6  ~| 


jkx2cos@ 


jkx„cos0 


2 

-o.  /2 


(4) 


2 

-a.  /2 


since  e 


,  and  where 


j  kx..  cos9 

_  — 1 

e 

sol 

jkx2cos@ 

s  „ 

e 

o2 

jkx^cosE 

• 

e 

soN 

(5) 


is  the  steering  vector  without  random  errors.  Equation  (4)  shows  the 

-a  2 /2 

steering  vector  £  is  not  biased;  the  constant  e  $  does  not  matter. 
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With  array  output  vector  V|  ,  and  steering  vector  £,  the  Hybrid  array 
with  constrained  mainbeam  gain,  as  shown  in  Fig.  1,  results  in  equivalent  steady 
state  weights  [  Sj 

W’  «  y  -1S*  (6) 


so  that 


e  =  V*  V  (7) 

o  —  — 

The  time  variable  is  omitted  but  implied.  In  (6),  the  random  variable  y  and 
matrix  M^'  are  given  by 


and 


N 


s1^'  _1s* 


M. '  =  CttjC 


(8) 

(9) 


"i  -  <i-fs>|«2|  s^s,T  +  Il»  +»„! 


(10) 


where 

2|  T  +  r  2, 

j 

with  Fg  being  the  complex  reference  loop  gain.  (Note  that  if  Fg  is 
complex,  is  not  Hermit ian) .  Substituting  (6)  and  (9)  into  (7)  we  have 


T  -1  -1  -1 

eQ  =  yV'  lC  C  S* 


(11) 


using  (1) ,  (11)  comes  out  to  be 


T  -1  -1 

e  =  u  V  M,  C  S* 
o—l  — 


(12) 


Define  the  equivalent  weight  vector  W  such  that 


T 

e  *  V  W 
o - 


Therefore,  W  *  uH^  ^S* 


(13) 

(14) 


Substituting  (9)  into  (8),  we  have 
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For  the  nonconstrained  Hybrid  array  processor  shown  in  Fig.  2,  the  equivalent 
weight  vector  is,  as  shown  in  Appendix  A, 

W  -  M1~1C“1S*  (16) 

The  average  SINR  can  be  defined  as  the  average  value  of  the  SINR 
treated  as  a  random  variable  which  is  a  function  of  the  random  gain  and 
phase  errors.  With  this  definition  the  average  SINR  of  these  two  system 
realizations  are  the  same.  However,  it  is  difficult  to  get,  by  analysis, 
an  explicit  expression  for  the  average  so  defined.  For  analysis  it  is 
easier  to  determine  the  separate  averages  of  the  power  of  signal,  interference, 
and  noise,  then  to  form  the  SINR  using  these  averages.  For  the  nonconstrained 
processor,  such  separate  calculations  are  especially  easy  since  p  is  one. 

We  therefore  use  this  approach  below,  setting  p  to  unity.  To  insure  that 
the  SINR  so  determined  is  meaningful  we  also  carry  out  an  average  SINR 
calculation  by  simulation  using  both  definitions,  and  we  find  the  difference 
between  the  two  quite  small.  Thus  we  conclude  that  the  SINR  obtained  by  the 
analytical  technique  we  have  chosen,  closely  matches  the  SINR  according  to 
the  first  definition,  and  apply  to  both  realizations  of  the  Hybrid  array. 

By  setting  u  to  be  one,  the  average  output  signal  power  is 


,-l  *  T  -1 


C  £  £  C  in  (17)  is  examined  in  the  following.  Since  C  is  a 
diagonal  matrix,  C  ^  is 


,-l 


1 

1+c" 


l+c„ 


‘  1 


1+c, 


Nj. 

For  cn<<l,  n«l,  2,  ...N.  (18)  can  be  approximated  by 


0 

z2  . 


-1  *  T  -1 


(18) 


(19) 


1_cnJ 

To  calculate  C  S^C  *,  take  the  expectation  with  respect  to  random  phase 
error  first. 


-I  T  -1 


1  „T  „-l 


(C  S  SC  )A  -  C  *  S  S*  C 


*  T 

We  examine  S  S  in  (20) 


*  T 
S  S 


N 

1 


S1  S2 


•••  sh] 


j  d>  jk(x2-x1)cos0 


e  e 


j  d>  jk(x1-x2)cos0 
e  •  e  1 


jk(x1-xN)cos0 

e  •  e 


(20) 


(equation  continued  on  next  page) 
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jk(x2-x1)cose 


jk(x1-x2)cos0  1 


jk^-x^cose 


+  (1-e  *  )I 


(21) 


Using  (5),  (21)  comes  out  to  be 


*  T  “04>  *  T  w  „ 

SS-e  ?S  S  +  (1-e  T  )I 
-  —o—o 


-O* 


(22) 


-0  i 


Comparing  (22)  with  (4),  we  see  that  (1-e  T  )I  in  (22)  is  contributed 
by  the  self-correlated  term,  which  is  similar  to  the  variance.  Substituting 
(22)  into  (20) ,  we  have 


(c'VsVV  -  e  *  C_1S  *S  TC_1+(l-e  ♦  )C'1C~1 

-  9  — O  — 0 

Then,  taking  the  expectation  with  respect  to  the  random  gain  error, 


(23) 


-o . 


-a. 


c'1S*STc"1  -  e  *  C-1S  *S  V1  +  (1-e  *  )C-1C_1 
-  — o  — o 


_"  V  *  t  -1 

We  examine  C  S  S  C  in  (24) . 


(24) 


(26) 


T  -1 
S  C 
— o 


(l-cl)sol  (l-c2)so2  .  .  .(1- cN)soNj 


Using  (25)  and  (26) ,  we  have 


-1  *  T  -1 

C  S  S  C 
— o  — o 


(l-cl.)  (1-cl^1“c2)sol  So2 


(l-cl)(1-c2)%iSo2*  (1-c2)< 


d-^xi-Voi* 


oN 


u-c3r 


(l-cN)- 


(27) 


The  expectation  of  (27)  is  therefore 

2 


-1  *  T  -1 

C  S  S  C 
o  o 


1+0 


Sol  3o2 


s  .s  *  l+o 
ol  o2  c 


b 


olSoN 


s  s 
ol  oN 


1+0 


*  T  2 
S  S  +  a  I 
— o  — o  c 


(28) 


C'V1  in  (24)  is 


c'V1 


(1-Ci)‘ 


0-_C2^ 


(1-CN)‘ 


(l+O  i 

c 


(29) 
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9 


Substituting  (28)  and  (29)  into  (24) ,  we  have 

_ 2  2 

-  a 


-1  *  T  -1 
C  S  SC 


-  a 


Y  (S  *S  T+a  2I)  +  (1-e  <{’)(l+o2)I 

— o  — o  c  c 

2  2 

C»S*ST+(1_e'a»)I+a  2I 
— o  — o  c 


(30) 


Substituting  (30)  into  (17),  we  have 

2 

-a . 


pd  =  !a 


T  -1  *  T  -IT*  * 

e  T  S ,  M.  S  S  M,  S, 

^-d  1  — o  — o  1  —a 

2 

T..  -1..  -1T*S  * 


+  (1-e  IS,  M,  Mi 


2  T  -1  -IT*  *1 
+  a.  S,  M,  ^  J 


'c-d*! 

milarly,  the  average  interference  output  power  is 


(31) 


J  T  r  T  * 

1  i  *  *\  t  \  A  •  •  A  A 


Pi "  ‘Jr  8i~  u- SL/j  -  -tj  *’ 


J  2,  -0<b  T  -1  *  T  -IT*  * 

Z'  r  9~  1 “  S  S  1M.  ST, 

-o-o  1  — Ij 


=  1 1 VI  I- 

j=l  J  J 


+  (l-e"°*  )SI.TM1-1M1-1T*SIj*  +  Oc2SIjTM1"1M1-1T*SIj*  ]  (32) 


The  average  noise  power  is 

pn  *  (ntw)(wtn)* 


T  -1  -1  *  T  -1  -IT*  * 
-  N  Mx  C  S  S  C  Ml  N 


(33) 


In  (33) ,  we  take  Che  expectation  with  respect  to  the  random  gain  and  phase 


errors  first,  giving 


e  a  \XWlTY  +  >»VV1TV 


2 


„  4>  T  -IT*  *  T  _i  *  -0. 

6  io«i  a  e\  \  ♦  a-e  ♦)»VV1T*a‘ 


Then,  we  take  the  expectation  with 


giving 


respect  to  the  noise  random  variables 


’■  ■ '  *  i.V1T*sY  »fV  *  a-/°*  > 

♦  °c2  ?vVl,V 


+  cr  ^Trace 
n  c 


[V1 


-IT*  • 

*1  j 


Note  that  the  quadratic  form  NTm  a  2  r  i  vr*~ 

-  M1  *1  N  reduces  to  o^Trace  Mj  V~lT  I 

because  the  noise  terms  are  independent  with  each  other.  ”  L  J 


Discussion 


“  Ch‘T‘  U  °aly  '™d°m  8al"  %  -  0,  Cl),  02,  (35) 


turn  out  to  be 


a  ■  ••VvVi.VV *  I^I-cVhi'S'1^ 


c-io 


(37) 


The  first  term  on  the  right  hand  side  of  each  of  the  equations  (39),  (40)  and 
(41)  is  similar  to  each  of  the  corresponding  terms  on  the  right  of  (36),  (37) 

~o  2 

and  (38)  except  for  the  factor  e  <t’  .  As  the  phase  variance  increases,  these 
terms  which  represent  the  mean  steering  vector  decrease.  At  the  same  time, 

the  second  term  on  the  right  of  each  of  the  equations  (39),  (40)  and  (41), 

-o  2 

containing  the  factor  (1-e  *  ),  increases  from  zero.  We  point  c  it  that  the 

assumption  of  normally  distributed  phase  errors  implies  errors  with  less  than 
say,  it / 2 . 


C-ll 


where  ReFg  is  defined  as  the  real  part  of  F  .  The  slgnal-to-nolse 
ratio  is  then 


SNR  = 


Y.e  N4Y  (1-e  *  )4Y  o  2 

_ a _ a _ a  c _ 

2  -a  2 

+(l-e  *  +ac2){Yd(N-l)-[ jl-Fg|2YdN+2-2ReFg]+l} 


(46) 


In  the  following  we  separately  discuss  the  effect  of  random  phase  error 
in  the  steering  vector  and  the  effect  of  random  gain  error. 

First  we  assume  that  there  is  only  random  phase  error.  Then  (46) 
reduces  to 


SNR  « 


2  2 

N+Y, (1-e  V  ) 
d 


2 

l+(l-e"a*  )Yd(N-l)Ul-Fg|2YdN+2-2ReFs! 


(47) 


Two  special  cases  are  of  interest.  With  Fg  *  1,  the  reference  signal 
is  ideal  and  we  have  from  (47) 

2  _o  2 

SNR  ®  Y.e  *  N+Y.  (1-e  *  )  (48) 

d  a 

Degradation  of  SNR  for  this  case  is  the  same  as  that  for  a  conventional 
beamforming  array. 

For  Fg  ■  0,  we  have  the  Applebaum  array  for  vMch  (47)  gives 


SNR 


2  2 
~o,  ~a* 

V  ^Yda-e  $  ) 

*  2 

l+(l-e  <p  )Y  (N-1)(Y  N+2) 


(49) 
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Equations  (47)  -  (49)  are  shown  plotted  in  Figures  3-5  for 


N  =  7,  30  and  100  respectively,  with  Fg  as  family  parameter.  The  improvements 
obtained  with  the  Hybrid  array  over  the  Applebaum  array  are  here  evident.  It 
should  be  noted,  though,  that  with  increasing  N  the  effect  of  a  given  phase 
error  variance  also  increases. 

Computer  simulation  was  also  carried  out  for  a  seven  element  10X  array 
without  interference  with  results  as  shown  in  Table  1.  Six  combinations 
of  F^  and  were  simulated,  each  with  100  samples  of  random  phase  errors.  The 

third  column  in  Table  1  was  obtained  by  compiling  the  SNR  for  each  random 
phase  sample  then  averaging  over  the  100  samples.  The  average  SNR  in  the 
fourth  column  is  the  ratio  of  average  signal  power  to  average  noise  power 
when  u  equals  one.  The  last  column  is  the  results  of  theoretical  computation 
based  on  equations  (47)  -  (49).  The  difference  between  the  third  and  the 
fourth  column  is  small.  For  the  Hybrid  array  with  F^  close  to  one,  they  are 
almost  the  same.  From  these  results  we  conclude  that  the  simulation  sample 
used  for  the  fourth  column  is  adequate,  and  that  the  similarity  of  results 
in  the  third  and  fourth  columns  suggests  that  the  SNR  calculated  as  a  ratio  of 
averages  and  as  an  average  of  ratios  are  not  appreciably  different. 

If  there  is  only  random  gain  error,  (46)  reduces  to 


SNR 


Y.N+ Y,o 
d  d  c 


l+o  {Y,(N-l)[|l-F  j  y,N+2-2R  F  ]+l) 
c  d  S  d  6  s 


With  F  *  1,  we  have  from  (50) 
s 


SNR 


YdN+*c  Yd 
1+0  2 


(50) 


(51) 
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Degradation  is  the  same  as  that  for  a  conventional  beamforming  array. 


i 

i 

I 


For  F  •  0,  (50)  becomes 


SNR 


Wc 


l+oc  [Yd(N-l) (YdN+2)+l] 


(52) 


Equation  (50)  -  (52)  are  shown  plotted  in  Figures  6-8.  The  Hybrid  array 
is  much  less  sensitive  to  the  random  gain  error  than  the  Applebaum  array. 

Computer  simulation  was  again  carried  out  for  the  same  array  geometry. 
Results,  as  shown  in  Table  2,  show  the  number  of  samples  is  adequate  and  SNR 
calculated  by  two  definitions  are  close. 


i 


Average  of  Each 
SNR  Sampled 


Computation  Based 
on  (47)  -  (49) 


Average  Signal 
Power  Over 


Average  Noise 
Power 

H 

18.42 

18.42 

18.42 

5° 

15.39 

15.22 

15.29 

5° 

3.73 

3.00 

3.14 

10° 

18.33 

18.33 

18.34 

0 

O 

r—i 

11.52 

11.08 

11.19 

10° 

-2.20 

-2.96 

-2.82 

TABLE  1 

Simulation  of  7  element  10X  Hybrid  Array  with  Random 
Phase  Errors. 


0.025 


0.95e~J10°  0.025 


0.025 


0.05 


0.95e"J1°°  0.05 


Average  of  Each 
SNR  Sampled 

! 

Average  Signal 
Power  Over 
Average  Noise 
Power 

Computation  Based 

on  (50)  -  (52) 

18.45 

18.45 

18.45 

18.08 

18.07 

18.09 

13.03 

12.66 

12.77 

18.44 

18.44 

18.44 

17.15 

17.10 

17.15 

8.22 

7.56 

7.73 

TABLE  2 

Simulation  of  7  element  10X  Hybrid  Array  with  Random 
Amplitude  Errors. 
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Figure  1  Constrained  Hybrid  Array  Processor 
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Figure  6  Sensitivity  to  Random  Amplitude  Errors 
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Appendix  Cl 

The  equivalent  weight  vector  of 
in  Fig.  2  is  derived  in  the  following 

equation  for  G  large  enough  is 

*  * 

V’  (e  -e  )  ■  S 
—  or  — 


a  nonconstrained  Hybrid  array  shown 
.  The  steady  state  control  loop 


(Cl-1) 


where 


and 


eQ  =  V'TW'  »  VTCW’ 


e  *  F  aS.  CW’  +  n 
r  s  ~d  —  r 


Using  (Cl-1)  and  (Cl-2)  we  have 


*  rr  ^ 

CV  V  CW'  -  CV  er 


Using (Cl-3)  V  er  in  (Cl-4)  is 


V*e  -  V*(F  qsJcw'+n  )  -  F  | ct‘ | Y ^W' 
—  r  s  — cl  —  r  s  — d  — d 


V"  t 

Substituting  (Cl-5)  into  (Cl-4),  we  have 


*  T  ,  i  2i  *  T  * 

CV  V  W'  -  CF  o  S .  S,1CW  -  S 
-  s  — d  — d  —  — 

(Cl-6)  can  be  written  as 
CM^CW'  =  S* 


(Cl-2) 

(Cl-3) 

(Cl-4) 

(Cl-5) 

(Cl-6) 

(Cl-7) 


Therefore 

W'  =  C~Hl1~1C_1S*  (Cl-8) 

The  equivalent  weight  vector  W  is  defined  such  that 

e  =  VTW  (Cl-9) 

O  ~  — 

From  (Cl-2)  and  (Cl-9),  we  have 

W  -  CW'  -  M1"1C~1S*  (Cl- 10) 
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Appendix  C2 

We  work  out  an  example  to  see  explicitly  how  the  desired  signal 
is  effected  by  the  random  gain  and  phase  errors.  Assume  that  there  are 
only  desired  signal  and  noise.  reduces  to 

»!  '  <l-Vl“2|S,V«„21 

The  inverse  of  can  be  obtained  by  using  the  matrix  inversion  lemma, 
which  gives 


(C2-1) 


-1  1  r,  (1  Fs)yd^d  ^d 

1  -  ~ 2  11 - 

o  (1-F  )y  ,N+1 

n  s  d 


M 


with  Yd  *  2 

on 


(C2-2) 


Assume  that  there  are  only  random  errors,  so  that  S  ■  S..  (31)  and  (35) 

— o  — d 

comes  out  to  be 

—  -a  2  T* 

■  2 1  d>  T  -1  *  T  -1  * 

Pd  -  |a  |e  9  ^  Mi 


-1T*  * 


and 


+  Jot2 1  (1-e  *  )SdTM1‘1M1“1  S 

+  i°2i°c2idTMrlM1"lT  *** 

2  ^  ^ 

P  =  o  2e  S,TM .-1  M  ”1S,*-KJ  2(l-e  *  )Trace[M  "V"1  ] 

nn  -dll-dn  11 


T* 

22  -1  -1 
+  o  a  TracefM,  M,  ] 
n  c  il 


(C2-3) 


(C2-4) 
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(C2-5) 


(C2-6) 


Substituting  (C2-5)  and  (C2-6)  into  (C2-3)  we  have 


la2|e"0*  N2 


|a[ (1-e  *  )N 


d  _4|/-_.„1.t|2+_4| 


°n  Ul-F8)YdN+ir  an4i(l-Fg)YdN+l|: 


| a2 | a  2N 

c _ 

5„4i<1-F»)^N+1l: 


Parts  of  (C2-4)  are  next  evaluated.  Thus 


T  -1^*  —1  *  1  T  (1~F_ ) YjS ,  S .  T* 

5d  1  «,  \  -Ai’li - “  ) 


(C2-7) 


^  1  1  -=d - 4  ^d 

o 

n 


(i-Fa)Y.S,VT  . 

[I - s--d~tl  )g* 


(1-F9>YdN+l 


(1~F  )  Y  N  +  l 
s  a 


.  _  (l-F  )*Y.NS.T  (1-F  )  Y.NS.* 

■~r  [s.T - ns.* - 3  d  ^  i 

°n  (1-Fs)  YdN+1  ^VY*1 


0  ^  (1-F  )  *Y,N+1  (1-F  )Y.N+1 

n  s  a  s  a 


(C2-8) 


7n  |(l-F8)YdN+l|* 
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and  also 


„  -L  -lT*  i  (1-F  >y.S.*S.T 

Trace[M1  ^  1  ]  -  ~  Trace[l - -g-  d~d 


(1-F  )Y.  N  +  1 
s  d 


U-V*Yd"  +  1  l(l-Fs)rd»  +  1|  2 


.  J_,N  .  “'W  .  (1-Fs>  V 

“„4  a-F5)VdS+l  (l-Fs)*YdN+l 


|d-Fs)YdN+l|: 


— r - -j  {n[|i-f  [2y  2n2+(i-f  )*y .n 

o  4  (1-F  )Y.N+ir  S  d  s  d 

n  s  d 


+  (l-Fs)YdN+l]-[|l-Fs|?Yd2N2  +  (l-Fs)YdN] 


“  Ul-Fs|2Yd2N2+(l-Fa)*YHN]  +  |1-FJ2Yj2N2} 


s'  'd 


s'  'd 


4 


Y  ,N(N-1) [ 1 1-F  |Zy  N+2-2R  F  ]+N 

^  S _ Q _  6  S 

0n''l(1-F,)''dN+1|2 


where  R  F  is  che  rnnl  r*rt  of  F  . 
e  r.  s 

Substituting  (C2-8)  and  (C2-9)  into  (C2-4) ,  we  have 


9  ■ 

a  2e  ’’  N 

n 


"  %*I»W*I* 


2  2 
an  (1-e  9  ){YdN’(M-l)[|l-FsrYdN+2-2ReFg]+N} 


a  !|  n  -v  )y.n+i|‘ 

n  s  d  1 


a  2 a  2{y..'KN-1)[  1 1-F  |  2y .N+2-2R  F  ]+N} 
n  c  d _ __  s'  d _ e  s 

an*l(1-VYdN+1|2 
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(C2-9) 


(C2-10) 


i 

3 


From  (C2-7)  and  C^-10),  the  output  SNR  is 


List  of  Symbols  -  Appendix  C 


C 


c 

n 

e 

o 

e 

r 

F 

s 

G 

I 

k 


N 


r 

N(t) 


R  F 
e  a 


h 

“IJ 


S 

— o 


on 


u 

V(t) 

V’(t) 

M 


diagonal  matrix  with  elements  1  +  c  ,  n  -  1,  ...N 

random  gain  error  on  the  n  array  element 

array  output 

reference  signal 

complex  loop  gain 

gain  of  the  weight  control  loops 

identity  matrix 

wave  number 

equivalent  covariance  matrix 
number  of  array  elements 
noise  component  of  reference  signal 
noise  vector  of  array  output 
average  output  signal  power 
average  output  interference  power 
average  output  noise  power 
real  part  of  F 

s 

steering  vector  with  random  phase  errors 
arrival  phase  vector  of  the  desired  signal 

arrival  phase  vector  of  the  jth  interferer 

tn  „  „ 

n  component  of  £ 

steering  vector  without  random  phase  errors 

th  *  „ 

n  component  of  S 
— o 

weight  factor;  a  random  variable  here 
array  output  vector  without  random  gain  error 
array  output  vector  with  random  gain  error 
equivalent  weight  vector 
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(Continued) 
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List  of  Symbols  -  Appendix  C 
(Continued) 


T 

weight  vector  defined  by  eQ  ■  V'  W' 
position  of  the  nC^  array  element 
desired  signal  waveform 
interferer  waveform 

signal  to  noise  ratio  at  array  elements 
steering  angle 

desired  signal  arrival  angle 
variance  of  random  gain  error 
noise  power  at  array  elements 
variance  of  the  random  phase  errors 
random  phase  error  on  the  nC^  element 


Appendix  D 

ESTIMATION  OF  THE  POINTING  VECTOR 


by 

Chlen  -  Chung  Yeh  and  Fred  Haber 


We  deal  here  with  the  problem  of  deriving  an  array  pointing  vector 
aimed  at  a  desired  source  assuming  that  the  environment  may  contain  an 
interferer  emitter  angularly  close  to  the  desired  source.  We  suppose 
that  a  reference  station  forms  part  of  the  system  that  station  providing 
system  management,  information  on  ground  source  location,  synchronizing 
signals,  and  also  capable  of  serving  as  a  beacon  for  organizing  the  array. 
The  scenario  assumed  is  shown  in  Figure  1.  The  array  shown  in  the  figure 
and  utilized  in  the  analysis  is  linear  though  there  is  no  Intention  to 
limit  the  approach  to  linear  arrays. 

The  procedure  to  be  examined  is  based  on  a  two-step  process.  First 
the  array  will  generate  a  pointing  vector  aimed  at  the  reference  beacon. 
This  vector  is  determined  by  correlating  array  outputs  with  the  output 
of  a  reflector  antenna  with  good  sidelobe  reduction  properties.  Then, 
with  given  information  on  angular  displacement  between  beacon  and  desired 
ground  source,  the  array  will  steer  this  vector,  aiming  it  at  the  desired 
ground  source.  The  objective  is  to  get  a  reasonably  accurate  pointing 
vector  toward  the  desired  ground  station  in  the  presence  of  strong  nearby 
interference  and  without  the  requirement  of  extremely  accurate  knowledge 
of  array  element  locations. 
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REFERENCE  ANTENNA 


S' 

« — © — o . a — © — © 


FIGURE  1 
Array  Scenario 
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1.  Estimation  of  Beacon  Pointing  Vector. 

With  the  arrangement  suggested  in  Figure  1,  assume  an  N  element 
linear  array  with  element  positions  x^,  ...  x^,  and  a  directive 

reference  antenna,  though  a  planar  array  is  likely  to  be  used  in 
applications,  the  concept  is  more  easily  pursued  using  a  linear  deploy¬ 
ment.  At  this  stage  the  desired  ground  source  is  assumed  off.  The 
beacon  signal  arrival  direction  is  denoted  6^  and  that  of  the  inter- 
ferer  is  denoted  0^.  The  inputs  of  the  array  elements,  which  are 
denoted  by  x^(t),  X2(t),  ...  x^Ct),  can  be  expressed  as 

j<|>.  jkx  cos9  j<j>  jkx  cose 

x  (t)  *  a(t)e  e  +  B(t)e  e  +  n  (t)  (1) 

n  n 

where  a(t)  is  the  beacon  signal  waveform,  B(t)  is  the  interferer  waveform, 

<f>k  is  the  electrical  angle  of  the  beacon  signal  at  the  origin,  is 

the  electrical  angle  of  the  interferer  at  the  origin  and  n  (t)  is  the 

n 

noise,  a(t),  B(t),  n^Ct),  »  and  are  assumed  to  be  independent. 

The  reference  from  the  directive  reference  antenna  is 


j0b  J*i 

r(t)  *  aa(t)e  +  bB(t)e  +  n^Ct)  (2) 

where  a  and  b  are  complex  numbers  representing  the  gains  of  the  reference 
antenna  to  signal  and  interferer.  The  vector  of  array  element  inputs  is 


X(t) 


x^t) 

Xj  (t) 


X„(t) 
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The  vector  is  pointing  vector  we  wish  to  estimate  and  we  consider 
the  estimator  S.  defined  by 

D 


s  -I 

T 


X(t)  r  (t)dt 


1 

T 


a(t)e  +  8 (t)e 


j  <J>  j  ^ 

*  jaa(t)e  b  +  b8(t)e  1  + 


%  +-(t^ 

nr(t)]  * 


(7) 


(7)  can  be  written  as 

44 


*  2  1 
a  a  (t)Sfedt  +  - 


b*82(t)STdt 


4 


r  j<f>h  *  * 

ja(t)e  ^[6  f(t)e  +  ny  (t)J 


+  8  (t)e  "Sj. 


-J4> 

a  a(t)e  “  + 


nr*(t)] 


+  N(t) 


“*  -i*b  *  -J* 

a  a(t)e  +  b  6(t)e 


1  +  nr*(t)j 


>  dt. 


(8) 


The  mean  and  variance  of  the  estimator  all  now  examined.  The 
expectation  of  the  first  term  on  the  right  hand  side  of  (8)  is 

T 


1  * 
A1  "  f  a 


<a2(t)>  dt 


*  2c 

a  a  S, 
—t> 


(9) 
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where 


<a2(t) \y  =  a2. 

The  expectation  of  the  second  term  is 

T 

a2  -  y  h\  ^2(t)?  dt 

o 

*  2 

-  b  B  Sx 

where 

<f8"(A)^>  =  82. 

The  expectation  of  the  third  term  is 

T  -  __ 

A3  *  T  ^<a(t)>  ^  bX^<e(t)Xe  3  1>+<nr*  (t)> 

'  \ 

o 

+  <^B(t)><'e  1  a*  <a(tf><fe  bN  +  <^nr  (t)'y 

‘  +  <H «=  )>  1_»*  <o )>  <e ' ^ * b>  + 

<„r*(t5>]dt. 


(10) 


(11) 


It  is  assumed  that  the  beacon  signal  is  a  zero  mean  time  continuous, 
binary  polarity  reversal  waveform  and  the  Interferer  is  monochromatic  at 
the  carrier  frequency  of  the  signal.  Then 


0 


8 
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(12) 


The  phases  and  are  uniformly  distributed  between  [0,  2rr], 
so  that 


<e  >"  0 

<e  l>-  0 

and  the  noises  are  zero  mean  so  that 
^nr(t^>-  0 
4?(t)^-  0 
with  (12),  (13),  and  (14) 


(13) 


(14) 


a3  =  0. 


(15) 


From  (9),  (10),  and  (15),  the  expectation  of  S,  is 


.  a  v  *  2  *  2 

<-b>‘“  “  4  +  b  S  *1 


(16) 


The  second  turn  on  the  right  represents  a  pointing  bias  which  should  be 
kept  small  relative  to  the  first  term.  That  is  we  will  want  to  make 


i  *  2  2  2 1 

|b  6  /a  a  ]<<  1  to  avoid  significant  bias  error.  We  now  turn  to  the 
variance  of  the  estimator  _Sfe.  For  two  vectors  A  and  B  with 


define 


a. 

b. 

1 

1 

a2 

and  B  * 

b2 

1 

• 

• 

• 

aM 

b 

N 

N 

-  — 

_  - 

AB 


a  b 
ii 

a2b2 


aNbN 


(17) 
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With  this  definition,  the  variance  of  S.  can  be  written  as 

D 

Var  [4]  -  <(4-  <4>)(4-  <4>)*> 
where  *  means  the  complex  conjugate.  From  (16),  (9)  and  (10), 
<h>  can  be  expressed  as 

T 


<%}•! 


T 

*  /  9  v  1  *  /  2 

a  <0  (t)\  4dt  +  y  b  ^(t^Sjdt. 


Since  it  is  assumed  that  the  signal  is  a  binary  code  sequence  and  the 
interferer  is  monochromatic , 

2 ,  .  2 


<a  (t)^>  -  a  (t)  ■  a' 

<B2(t))>  =■  62(t)  -  e2 


Substitute  (20)  into  (19), 

*  1 
(h}  “  T  , 

From  (8)  and  (21), 


§b 


f  T 

*  2  1  *  2 
a  a  (t)4dt  +  y  JQ  b  B^(t)SIdt. 


* .  if  ^b  r *  ~^i  *  ~i 

-  <4>  m  T  Jo  '“<t>e  4  [b  8<'>*  +  "r  <e)J 

j*.  r*  *  1 

+  B(t)e  la  a(t)e  +  n  (t)J 

ri  *  T' 

+  nr  (t)|  )  dt 


+N(t) 


1 

T 


*  JYb  * 

a  a(t)e  +  b  B(t)e 


it.  * 

a(t)e  “4b  B (t)e  dt 


(18) 


(19) 


(20) 


(21) 


it,  *  -it. 

8(t)e  Si  a  o(t)e  dt 
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(t)dt 


+  Y  1  [_a(t)e  ib  +  0(t)e  -l  “t 


1  I"  .  -j4>j  *  ' 

^  N(t)  jji  a(t)e  b  +  b  B(t)e  +  **r  (C)J  <*t 


Substitute  (22)  into  (18) 


Vat  (Sb)  -  <(4  -<4>K4  -<ib>)*> 

X.  T  JK  ,  -1*! 

.<7±  a(t)e  bSbb  B(t)e  dt 

1  rt 


i  ^i  *  ~j*b 

+  i  jo  B (t)e  a(t)e  dt 


Tp  j  4>  -i  A 

+  |  a(t)e  ^  +  6(t)e  nr  (t)dt 

Jo  L  J 


t  r  -3$,  * 

+  ^  fl(t)  3  Cx^t^e  +  b  B^t^e  +  nr  ^  dt 


[i  (T 
<  —  a 

T 

rx 


-j$K  * 

(t)e  °Sij  bB  (t)e  dt 


i  f  T  *  j*b 

+  i  6 (t)e  Sx  aa ( t ) e  dt 


i  T  r  -j*h  *  ~^±  *\  ,  . J 

+i  I a(t)e  bS^  +  6(t)e  £x  nr(t)dt 

+  4  N* (t)  Taa(t)e  b  +  b8(t)e  1  +  n  (t)|  dtW 


■<K 


(1)  -N  , 

T 


JO, 


-J0< 


a(t)e  “Sfeb  6 (t)e  Adt 


"JO.  *  JO,  V 

a(t')e  b6(t‘)e  cit’^ 


X  T 

/I  j*i  *  ~J*h  “JO,  *  J0K  v 

(2)  +\~ ?  S(t)e  S  a  a(t)e  °dt  •  8(t')e  *■«  aa(t')e  bdt\> 

T  o  Jo  1 


(3)  <i 


t  T 


j*b  * 
a(t)e  S,b  8(t)e  dt 


6(t’)e  j  1SI*aa(t’)ej*bdt^> 


r  t 


/i  i  Jo,  *  -Joh 

(4)  +C-,  S(t)e  Va  a(t)e  °df 

T  Jo  1 


—1<t>  i* 

a(t')e  bs,  bB(t')e  Sit' 


> 


/i  i  r  j^h 

(5)  +<(i  1  1  -■/‘x-  b 


JO 


j  a(t)e  uSh  +  B(t)e  ^  !  nr*(t)dt 

L  J 


T  _ 


■J0V 


a(t')e  b5b*  +  8(t')e  J0i§1*  nr(t')dt^> 


<12 


r  t 


-J0b  *  -jo 

N(t)  a  a(t)e  +  b  B(t)e 


N*(t' )  j^aa(t' 


J0b  j* 

)e  +  bB(t')  e 


1  +  nr*(t>j  dt 
1  +  dt'> 


(7) 


T- 


J4> 


“JO, 


JTb  *  JY< 

a(t)e  uSbb  B(t)e  <it  + 


JO) 


“JO, 


B(t)e  "Sja  a(t)e  “dt] 


r  “J0b  *  -Joi  * 

[a(t')e  DSfe  +  B (t')e  Sj  ]nr(t')dt' 


’  *  J0h  JO, 

N  (t')[aa(t')e  D  +  bB(t')e  1  +  n(t')]dt'}> 
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|a(t)e  + 


S(t)e  S^J  nr  (t)dt 


N(t)  £ 


‘  *  -iK  *  -i*±  *  "I  1 

a  a(t)e  +  b  B(t)e  +  t»r  (t)  dt) 


T  T 

-j<j>  j<()  -j$ 

a(t’)e  ^  b6(t’)e  6(t')e  VaoU'Je 

0 


*(t)e  +  8  (t)e  s3nr*(t)dt 


I  (t')  aa(t')e  D  +  bB(t’)e  1  +  n^t')  dt'^> 


r  *  -js  * 

N(t)  a  a(t)e  0  +  b  B(t)e  1  +  n 


nr  (t)  dt 


r  -j*b  * 

a(t')e  DSl  +  B(t')e  V 


n  (t')dt'\ 


No  Noise  Case: 


In  practice,  the  signal  is  designed  to  be  much  larger  than  the 
noise.  It  would  be  reasonable  to  assume  that  the  noise  is  zero.  Then, 
only  the  first  four  terms  on  the  right  hand  side  of  (23)  are  not  zero. 
The  first  term,  in  (23),  is 


fT 

o(t)a(t')|b!2I1B(t)3(t,)dt  dt'^> 
o 


T  r  T 


<a(t)a(t')>B2|b|2J1dt  dt' 


where 


T  J  o  Jo 

,  *[  i } 


is  an  N  x  1  all  l's  vector.  a(t) 
of  rec.tangular  chips  with  random  and  equiprobable 


(24) 

is  a  binary  code  sequence 
amplitudes  ±  a. 
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IC  can  be  shown  for  such  a  waveform 


<o(t)a(t')>  -  }  a2[l  -  1  —  ~  —  -1-1  for  |t  -  t'|  a  t 


otherwise 


(25) 


For  T  >>  t  (24)  can  be  approximated  to  be 


t'-O 


<a  (t)a(t'  )>82  |  b  ]  ^JC^dtdt ' 


t-0 


t'-O 


a2  [1-  -^^-]S2|b|2l1dtdt' 

t-t’-T 


t  T 


a2t82 1  b|  2L,dt' 

t'-O 


1  2  q2 i  .  i 2_ 

T  a  t8  |b|  1^ 


(26) 


Similarly,  the  second  term  in  (23)  is 


1 

< — 


T  r  T 


2 

T  ■’t'-O 


t-0 


8(t)8(t')IjJa|  a (t)a (t '  )dtdt ' > 


t'-O 


t-0 


S  2 | a | 2<a(t)a(t ' )>I^dtdt ' 


t'-O 


B2t|a|2a2I1dt’ 


1  2  a2 i  | 2_ 

Y  a  T8  | a  | 


(27) 
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The  third  term  is 


<1 


T 

t'-O 


1  j2t  — J  2<|> 

.  ab*S.  S.*e  be  ia(t)o(f)6(t)6(t’)dtdt’> 
t*U  o  1 


t'-O 


t-0 


*  *  j2*  -J2*.  2 

ab  SVS_.  <e  b><e  ><a(t)a(t')>8  dtdt' 
— b— I 


(28) 


because 


j2*  -J2* 

<e  D>  *  <e  >  *  0. 


Similarly,  the  fourth  term  is 


1 

T 


2 


t'-O 


t-0 


*  *  J2*i 

a  bS^  Sj.e  e 


T 

t'-O 


t-0 


*  *  -JZ»b 

a  bS^  j3j<e  ><e 


a(t)a(t')8(t)B(t')dtdt'> 

j2*.  2 

l><a(t)a(t')>8  dtdt' 


(29) 


Substituting  (26),  (27),  (28)  and  (29)  into  (23),  VarlS^]  for  the  no 
noise  case  is 


VartSfa] 


(30) 
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where 


|a|  ■  amplitude  of  signal  at  reference  antenna  output 
relative  to  element  signal  level 

| b |  •  amplitude  of  interference  at  antenna  output 
relative  to  element  interference  level 

t  -  PN  sequence  chip  width 

T  ■  integration  duration 

2 

o  /2  ■  signal  power  at  each  element 
2 

8  /2  -  interference  power  at  each  element 


Example : 

We  obtain  numerical  values  for  the  estimator  bias  in  the  no  noise 
case.  Assume  at  the  inputs  of  array  elements  the  signal  is  lOdB  lower 
than  the  interference,  and  at  the  reference  from  the  dish  antenna  the 
signal  is  20dB  higher  than  the  interference;  i.e.  the  antenna  gain  is 
30dB.  Thus 

8 2  -  10a2 

|a|2a2  -  100|b|282 

so  that 

| a | 2  -  1000 | b | 2 

From  (16) 

*  *  2  *  2 
<Sb>  *  a  a  +  b  B  S j 

-a  + 

a  a 


where 


i  *  2 1 

lb  8  1 

I  *  2I 
|  a  a  | 
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The  biasing  error  vector  is  l//lO  of  the  desired  pointing  vector 
magnitude.  This  may  not  be  good  enough.  If  at  the  inputs  of  array 
elements  the  signal  is  the  same  as  the  interference  (a  condition  which 
is  not  unreasonable  if  the  pilot  source  is  a  fixed  station  capable  of 
substantial  power)  and  the  antenna  gain  is  30dB. 


2 


a 


8 


2 


and 

|  a } 2  =■  1000 1  b  | 2 


a*«2[Sb  +  Vi  V 

a  a 


where 


i  *  2 1 

lb Vl _ L_ 

]  a*a2 1  10  /TO 


This  ratio  is  much  better  though  we  would  advise  using  a  dish  with  better 
than  30dB  mainlobe  to  sidelobe  ratio  to  get  even  greater  suppression  of 
the  bias  term. 


With  Noise: 


Let  the  power  spectral  density  of  the  noise  at  each  array  element 

be  denoted  <f>  (f)  and  assume  it  to  be  as  shown  below 
n 
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and 


r 


The  variance  of  S,  is  given  in  (23).  Here  the  first  four  terms  are  the 
~o 

same  as  those  for  the  no  noise  case,  which  are  given  in  (26),  (27),  (28) 
and  (29).  The  fifth  term  in  (23)  is 


f  T 


1 

< — ■ 


j*h  3*±  * 

[a(t)e  DS,  +  6(t)e  S_]n  (t)dt 
- d  — I  r 


-J4>k  *  -j**  * 

[a(t')e  V  +  8 (t ' )e  V  ]n  (t')dt’> 
)  — o  i  r 


T 

t'=0 


t=0 


<a(t)a(t')><nr(t')nr  (t)>I^dtdt' 


T  (1 
t’-O 


ta0<B(t)e(t,)><nr(t,)nr  (t)>_I^dtdt ' 


(33) 


substituting  (25),  (31)  and  (32)  into  (33)  and  using  T  >>  t,  (33)  can 
be  approximated  by 


1 

T 


2 


rT 

t'-O 


t'+T 


t-f-T  “  t1' 


sin 


2n(t-t') 


2TT(t— t') 
T2 


Ijdtdt' 


T 

t '  =0 


rT 

t=o 


BV  -  Vtdf 


(34) 


Examining  the  two  terms  in  (34)  separately,  the  first  term  *s 


t’-O 


t'+x 


sin 


2  ir  (t-t ' ) 


2  r,  t-t’  ,  2  ‘2 

t-t ' -f  a  2irTt-t~' ) 


I, dtdt ' 


D-17 


I'T 


t'-O 


sin- 


2irt 


.2a-  M1o  2  -J2. 


.  -  a-  -^lo  ■  -  r-i.dtdf 

t— t  t  r  2nt  — 1 


1 

T 


.  2nt 

t  sin - 

a2[l-  i^-lo  2  2-  I  dt 

t»-T  t  r  2irt  — 1 


sin: 


2irt 


|  Jo  “2[l-  <>°r2  -rkh* 


(35) 


In  the  second  term  of  (34)  let  t-t'  =  u  and  t'  -  v,  giving 


sin- 


2it  (t-t 1 ) 


t'-O 


2  2 
S  o 


TT— Ldtdt* 


t-0 


r  2tt  (t-t* )  *1 

t2 

(0  ('T  ~ 

6  o  2  I^dvdu 

u— T  v--u  r  2nu 


*K 


u-0 


T-u  sin— - 

2  2  t2  • 

n  8  o  — = -  I.  dvdu 

v-0  r  2ttu  —1 

t- 


t  T 


siir 


2ttu 


o'1-'*162^2  TiT  Iid“ 
t2 


(36) 


From  (35)  and  (36),  the  fifth  term  in  (23)  can  be  replaced  by  the 
approximate  equivalent 

2irt 


2 

T 


t  sin- 

2 . .  t,  2 
0  a  t^°r  2itt 

x2 


I,  dt 
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X 


2*TTt 


Iidt 


fT 

o 


sin- 


2irt 


(T-t)  |  b 


2  2 
6  a 


2jrt  -1 

T 


I,dt 


rT 

0 


(T-t)o 


The  7th,  8th,  9th  and  10th  terms  in  (23)  are  zero  because  of  the 
independence  of  the  random  variables  representing  signal  and  interference 
phase.  Substituting  (26),  (27),  (37)  and  (38)  into  (23),  we  have 


Var  [S^]  -  y  <*2T82(|a|2  +  [  b|  2)ix 


2irt 


+  l 


2 r ,  t,  2  l2  _  , 

0  a  [1"  t  °r  kdt 


+  k 


2 

2wt 


T  s  in-  t 

0  <T-t)B V 


*1 


rT 


sin 


2irt 


,  |2  2r,  t  ,  2  T1  . 

0  |a|  a  [1-  T  ]on  -27F  ILdt 


1 
2tt  t 


+ 


T  sin' 

n  <T-t)|b|2B2o  2 


I,dt 


*K 


a 


sin- 


2irt  ,  2ut 


2  2  T1  T2 
rt  (T-t)o  a  1 - - - 5 — —  I.dt 

0  »  r  ^  “l 

t1t2 


To  simplify  (38a)  somewhat  assume  the  rf  bandwidths  of  array  elements 
and  reference  antenna  equal  so  that  “  T2  <<  T>  we  c^en  have 


Var(Sb]  -  |  a2TB2(|a|2  +  ibj2)^ 


fT 


*1 


sit*- 


2irt 


2 .  2  .  ,  1 2  2,  tv  1  . 

0  a  (or  +  I  a |  on  )(1-  ~) 


sin 


1 
2irt 


0  82(ar2+|b|2on2)(T-t)-^-I1dt 


+  k 


T  sin2  ^r- 

o  V  -T-T--tidt 

(2zt) 

T1 


To  examine  (39),  the  second  term  on  the  right  hand  side  of  (39)  is 

2nt 


2 

T 


sin 


lalVm-f>-2^-ildt 


2irt 


2  2.  2  |  .2  2 

-  a  (of  +  I  a  I  on  )^ 


0  2*t 

T. 


dt 


(38a) 


(39) 
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2  2/  2^,  ,2^2.1  T  T1 

T  v  r  1  1  n  2it  1  t 


2wt 

0  T, 


dt 


In  (39a),  two  integral  terms  are  examined  below 

2irt 


sin- 

fT  T, 


2irt 


dt 


J. 

2ir 


2ttt 

^  3inx 

0  x 


dx 


J. 

2ir 


sinx 
0  x 


dx 


The  other  integral  term  is 


rr 


sin^1Tt  dt  ■  0 
0  T1 


(39a) 


(39b) 


(39c) 


for  t  -  mTi«  ®  is  a  integer.  The  second  term  in  (39a)  will  also  be 
negligible  for  other  values  of  t  provided  that  t  >>  Substituting 

(39b)  and  (39c)  into  (39a),  the  second  term  on  the  right  hand  side  of 
(39)  is 
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T  J0°  (°r  +  !al  °n  )(1“T)  "lit—  h 


T,  2  ,  2  I  |2  2V  T 

_1  o  (a  +  | a  I  o  )  I. 
2T  n  1 


The  third  term  on  the  right  hand  side  of  (39)  is 


i  S2(»t2  +  |b|2«„V-t>  -2^-  h 


3  82(«r2+  l"lV>Til  )  -pfdt 


T  .  2irt 
sin - 


In  (41) 


2  .2 ,  2  | . 2  2.  T1  _  J_2irt 

3  6  (or  +  |b|  °n  ^  2tT  — 1  SinT7d,: 

T  'o  1 


dt”0  as  specified  by  (39c),  (41)  then  reduces  to 


•T  sin - 

3  62<°t2  *  W 2«„2)tt-t)  -jji-  4  dt 

T  - 

;0  T, 


(-T  sit?— — 

2  „2 ,  2  .  u  .2  2.  „  T1  T1  2 nt 


?  *W V>  h  27 

'o  -  1 


Tl„2,  2  .  il|2  2,,  sinx 

—  6  (o_  +  b  o_  )I,  — — 


(42) 


2*  *V  +  l»l: 


The  fourth  term  on  the  right  hand  side  of  (39)  will  be  approximated  by 
an  upper  bound,  as  follows: 


.  2  2rt 
_  «  sin  - 

«-'>“«  °r  - li  Jl 

<^) 

T1 


2  2  2  . 

<  —  o  a  I, 
T  n  r  -1 


T  sin2^ 

_ ]l  dt 

2 

o  -27rt, 

T1 


2  2  2  _  ll 

T  °n  °r  -1  2ir 


r2jl T  2 

t,  sin  x  , 

1  — 2 —  dx 

x 


T1  2  2 
itT  n  r  — 1 


sin2x 


o  x 


dx 


Substituting  (40),  (42)  and  (43)  into  (39),  we  get  our  approximation 
for  Var[S,  ]. 

— D 


(43) 


Var(Sb]  *  ^  ctV(|a|2  +  M2)^ 
+  ^a2(or2+|a|2on2)Ii 
+  ^32(or2+|b|2on2)Ii 


(44) 
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To  examine  the  estimated  more  thoroughly,  we  go  back  to  (8). 
can  be  written  as 

rT 


(8) 


rT 


*  2  1 
a  a  (t)S^dt  +  y 


bV(t)SIdt 


+  hk 


* 

a(t)e  D[b  8 (t)e  +  nr  (t)]dt 


*  *lk 


j<K  *  -J*h  * 

8(t)e  x[a  a(t)e  +  nr  (c)Jdt 


♦I 


N(t)[a  a(t)e  +  b  B(t)e  +  n^  (t)]dt 


(45) 


The  2nd,  3rd  and  4th  terms  in  (45)  are  zero  mean  with  variance  not 
zero  when  T  is  finite.  However,  the  2nd  and  3rd  terms  can  be  written  as 

and  respectively  where  and  C ^  are  uncorrelated  random  variables. 

So  the  2nd  term  contains  signal  information  and  the  3rd  term  contains 
interferance  Information.  The  4th  term  is  random  fluctuation. 

From  (44)  the  variance  contributed  by  the  4th  term  of  (45)  is 


T_1 

2T 


2  2 
a 

n 


2 


Ib| 


2  2 
o 

n 


+ 


li 

2T 


2 


(45)  can  be  written 


(46) 
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^  -  (aV  +  C1)Sfe  +  (bV  +  C^Sj  + 


Jwi 


Ul6 


V 


jwK 


u  e 
N 


(47) 


In  (47),  u  u  .  .  .  u  are  independent  random  variables  with  zero 
1  Z  N 

.  .  22  22 

mean  and  variance  o  «  a  ■  0  =  a  ,  . 

U1  U2  “n  rf 


(48) 


wl’  w2’  ‘  '  ‘  WN  3re  in<*ePen<*ent  random  variables  uniformly  distributed 

between  [0,  2ir].  The  u  's  are  independent  of  the  w  '  s. 

n  n 

From  (44)  and  (45)  we  can  see  that  the  variance  of  is 


2„2 


B‘|b|  +2Ta° 


2  2 


(49) 


and  the  variance  of  C 2  is 


i  .2  2, 

P  01  | 


.  1  .2  2 

+  2T  6  °r 


(50) 


*  2 

Since  the  standard  deviation  of  is  much  smaller  than  |a  a  |,  can 
be  ignored  in  (47).  Thus  we  have 


4 


*  2 
a  a 


<4 


*  2 

b  B  +  C2 

*~2 
a  a 


+ 


jwl 

*  2  e 


(51) 


The 


*  2 

constant  a  a 


may  be  dropped. 


*  2 
a  a 


(51)  becomes 


D-26 


bV  +  C2 

4-4  +  — — lx  + 

a  a 


V 


V 


K 


(52) 


where  v 

N 


*  2 
a  a 


is  a  r.v.  with  zero  mean  and  variance 


I  I2  4 
a  o 


2 


,  1  2 1  1 2  2 

[2Ta  |a|  4 


2 


I  b| 


2  2 
a 

n 


2  2, 


(53) 


2.  Scanning  from  Beacon  to  Ground  Station. 

We  now  turn  to  the  problem  of  steering  the  beacon  phase  vector  to 
the  direction  of  the  ground  station.  From  (16)  &e  have 


"V 


*  2 
b  B 

4  +~^2  h 

a  a 


(54) 


2  2 

SL  and  £  are  beacon  and  interferer  phase  vector,  respectively.  a  and  6 
are  the  power  of  beacon  and  interferer.  a  and  b  are  complex  numbers  re¬ 
presenting  the  voltage  gain  of  the  dish  antenna  on  beacon  signal  and 
interferer,  respectively. 
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To  examine  (54)  more  thoroughly,  define 


*  2 

.  .  A  b  6 

Y  "  YRe  +  jYIm  "  — 
a  a 


(55) 


Then, 


.§b  +  yJLj  • 


(56) 


It  is  reasonable  to  assume  that  beacon  and  interferer  have  the  same 


order  of  power  and  that  the  dish  antenna  has  a  30  dB  mainlobe  to  side- 

2 

lobe  ratio  i.e.  3  8^  and  3  10^.  Then  |y|  is  in  the  order  of  10 


From  (5b), 


|b| 


jkx, cosQ,  jkx, cos0. 

Id,  11 

e  +  ye 

jkx2coseb  jkx2cos01 
e  +  ye 


JkXjjCosOh  +  jkxNcos@1 
i  ye 


(57) 


where  0,  and  0.  are  beacon  and  interferer  directions. 

□  i 

To  see  the  phase  of  the  nth  element  of  <S,  >  we  write 


jkx  cos0,  jkx  cos©, 

„  n  b  .  J  n  i 
e  +  ye 


jkx  cos0,  jkx  (cos0.  -  cos0, ) 

n  l)r«  *  n  i  b  « 

e  [1  +  ye  ]. 


(58) 


The  term  in  brackets  of  (58)  is 


1  +  Ye 


jka  (cosQ.  -  cos0.  ) 
n  l  b 
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■  1  +  a„  cosfkx  (cos0.  -  cos9, )]  -  y_  sin[kx  (cos9.  -  cos9,  )] 

Re  n  l  b  Im  n  l  b 

+  j{  Y„  sin[kx  (cos9  -  cos9.  )]  +  yT  cos[kx  (cos0.  -  cos0. )]}  (59) 

Re  n  i  b  Xm  n  i  b 

j*_ 

-  A  e 
n 

The  phase  term  il>  is  what  we  are  interested  in  and  is  specified  by 

y_  sinfkx  (cos9.  -  cos9, )]  +  yT  cosfkx  (cos0.  -  cos9, )] 

Re  n  l  b  -Lin  n  i  b 

tantj^  =*  -  (60) 

n  1+y „  cosfkx  (cos0.  -  cos0.  )]  -  yT  sin(kx  (cos9.  -  cos9. )] 

Re  n  l  b  xm  n  i  b 

For  | y | <<1 ,  (60)  can  be  approximated  by 


tani(/  2  y_  sin [kx  (cos0.  -  cos9,  )]  +  yT  cos[kx  (cos9,  -  cos9,  )]  (61) 

n  Re  n  i  blmni  b 


Since  taniji  is  small,  tan  i|*  *  . 

n  n  n 


>/<  =  YD„sin[kx  (cos8  -  cos9,  )]  +  y  cosfkx  (cos9,  -  cos0  ) ]  (62) 

n  n  i  d  im  n  l  b 


From  (58)  and  (59),  the  nth  element  of  <S,>  is 

— D 


jkx  cos9  jkx  cos9.  j  (kx  cos9,  +  ^  ) 

_  n  „  n  i  .  „  n  b  n 

e  +  ye  =  A  e 

n 


(63) 


The  normalized  mean  value  of  the  estimated  beacon  phase  vector  is  therefore 
j(kx1cos9b  +  <1^) 


j(kx2cos9b  +  1 1>2) 


j(kxNcos0b  +  ^N) 


(64) 
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In  (64),  i|»^,  ip^,  ...  i|/N  are  phase  errors.  (62)  shows  that  ipn  is  correlated 

with  kx  cos8.  which  is  the  correct  pointing  phase.  However  is  more  like 
n  □  n 

a  random  phase  error  as  the  following  example  shows.  The  phase  error  of  the 
nth  element  as  given  by  (62)  involves  trigonometric  functions  of  [kx  (cos 


9.  -  cos0. )].  Thus  we  need  only  consider  the  mod-2ir  version  of  this 
i  b 

quantity,  assume  both  0.  and  0.  are  near  ^  rad  (both  beacon  and  interferer 

lb  / 

are  near  broadside  and  9.-9.  -  -i^r).  Then 

1  b  ZU 


COS0.  -  COS0,  “ 
i  b 


sin(j  -  0±)  -  sin  (j  -  0b) 


Assume  an  array  of  10  elements  with  array  size  ~  500  wavelengths,  the 
average  spacing  between  elements  is  50  wavelengths.  Then 


kAxA0  =  2rr  *  2.5 


Adjacent  elements  will  therefore  differ  by  the  order  of  5tt  radiuses  and 
the  modulo  2tt  phase  differently  will  tend  to  be  random.  With  non- 
uniform  spacing  the  modulo  2n  phase  differences  will  be  further  ran¬ 
domized.  Therefore  may  approximately  be  treated  as  a  random  variable 
with  zero  mean  and  variance 

2  2 

>  (take  expectation  over  xn) 
n 

2  2  2  2 

-  <YRgsin  [kxn(co80i  -  cos0b)]  +  YIm  cos  [kxn(cos0i  -  co0b>] 

+  2YReYIm3in^kxn^cosei  "  cos6b)  ]cos[kxn(cos0jL  -  cos@b)]> 

1  2  1  2 
2  YRe  +  2  YIm 

M  \  M2  (65) 


D-30 


In  (65)  mod  (kx  (cos0.  -  cos9.))  is  assumed  to  be  uniformly  distributed 
n  i  d 

between  [0,  2ir]. 

If  the  integration  time  in  estimating  is  finite,  the  estimated 

beacon  phase  vector  is  given  by  (52) j-  ,  ~ 

JW1 

V 

be2  +  c  j», 

4*4  +  —  2 "  4  +  V 

a  a 


where  is  zero  mean  with  variance 


2  t  a2  2,  ,2  ^  1  2  2 

%  8  a  |al  +2T  8  °r 


and  where  v  ,  v  ,  ...  v  are  independent  random  variable  with  zero  mean 
1  £  N 


and  variance 


2  1  /l  2,  ,2  2  .  T1  „2i.  1 2  2 

Tv  "  7T2 T  2T  01  |al  °n  8  |bl  °n 

n  a  a 


,  T1  22. 
+  2T  °n  °r  ] 


Wl’  W2’  ***  WN  3re  in^ePen<^ent  random  variables  uniformly  distributed 

between  [o,  2tt ] .  The  v  's  are  independent  of  the  w  's. 

n  n 


Define  -  -r— „  »  X  =■  X  +  IX  T 

*  2  re  J  I® 

a  a 


X  is  zero  mean  with  variance 

2 


I  i  2  A 
a  a 
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(70) 


7 


1*1 


1  r x  a2  2,  , 2  ^  T1  .2  2. 

T4  l?  6  °  lal  +  2T  e  °r  1 


using  (55)  and  (69),  (66)  becomes 


4-4+  (y+*>4  + 


V 


V2e 


V 


N 


_l 


(71) 


The  nth  element  of  4  ia 


jkx  cos0, 
n  b 


jkx  cos0 


jv 


+  (Y  +  X)e  “  *  +  v  e 


n  i  +  v  n 
n 


jkx  cos0  jkx  (cose.  -  cose.)  j (w  -kx  cos6  ) 

e  °  b[l+(Y+X)e  n  1  b+vne  n  n  b](72) 


Let  v  *  v  „  +  jv  _  and  define 

n  nRe  J  nlm 


J*' 

A'  e  “  =  1  +  (y  +  X)e 
n 


jkx  (cos9  -cos0,  )  j (w  -kx  cose  ) 


b  ~  '  n  n 

+  v  e 
n 


(73) 


Proceeding  as  was  done  in  going  from  (7)  to  (9)  we  have 


Ip'  -  (yd  +  X_Q)  sin[kx  (cos6  .-cos6  )  ] 
n  Re  Re  nib 

+  (y  .  +  )  cos[kx  (cos0  .-cosQ,  )] 

Im  Im  nib 

+  v  sin(w  -kx  cose, )  +  v  _  cos(w  -kx  cos0. )  (74) 

nRe  nn  b  n  In  nn  b 
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The  normalized  beacon  phase  vector  estimated  is 
j_j(kx1cos0b  +  irJ 


4-  e 


j(kx2cos6b  + 


j(kVoseb  +  *V 


As  In  Che  case  of  we  can  argue  Chat  ip'a  can  be  Created  as  a  random  error 
ip'  is  zero  mean  with  variance 


V  -<(*'nr> 

n 


1  2  1  2  1,2,1,2 

—  v  +  —  v_  +  -r-  <X_  >  +  —  <X,  > 

2  Re  2  ’  Im  2  Re  2  Im 


1  2  1  2 
+  -=-<v  >  +  t-  <v  T  > 

2  n  2  n  Im 

Ke 


1.  ,2^1  21  , 

yWl  +2°\  +2°v2 


1  bV 

2  !a|V 


+  TT2-4  £  32a2|a|2+^BV) 

2  a  a 


1  rTl  2,  ,2  2  .  T1  „2 , ,  ,2  2 

+  ^TTT  2T  0  °n  +  2T  6  >bl  an 
2  a  a 


1  2  2. 

+  2T  °n  °r  1 


In  (76),  the  dominant  terms  are 


bV 

|a|V 


2  a 


t  .2  2,  , 2 

I4T6  °  |al 


In  order  to  form  a  beam  in  the  direction  of  the  beacon  station, 

the  array  processor  will  put  a  phase  shifting  vector  which  is  the  con- 

*  * 

jugate  of  the  estimated  vector,  i.e.  . 

-j(kx1cos9b  +  <1^) 

-j(kx2cos0b  +  ^'2^ 

-jCkXjjeosOh  +  <rN> 

ft 

We  now  turn  to  the  problem  of  steering  the  pointing  vector  from  the 
beacon  to  the  desired  source.  The  basis  for  this  approach  is  developed 
in  Steinberg  [  1  ]  pp.  246-249. 

Assume  we  have  a  linear  array  on  the  X  axis  with  element  positions 
imperfectly  known.  The  beacon  station  is  in  the  direction  ©b  as  shown 
in  the  figure  below. 


'  beacon 


The  steering  vector  for  forming  a  beam  on  the  beacon  is  given  by 

(77);  the  nth  element  phase  shift  is  given  by  -kx  cos0.-i|i*.  We  now 

want  to  move  the  beam  through  an  angle  <|>  from  9^.  A  phase  shift  ? 

needs  to  be  added  to  -kx  cos9, -<fi'  ,  f  being  determined  in  the  following 

n  d  n  n 

way. 

The  array  processor  imperfectly  knows  that  the  location  of  the  nth 

element  is  x  +  Sx  where  5x  is  the  location  error.  To  form  a  beam  in 
n  n  n 

9  ,  the  phase  shift  calculated  by  the  array  processor  is  -k(x  +6x  )cos0.. 
d  n  n  b 

to  form  a  beam  in  8^  +  the  phase  shift  calculated  by  processor  is 
-k(Xi+ 6Xn)cos(9b  +  <J>).  The  phase  increment  determined  by  array  processor 
is,  therefore, 

4*  “  -k(x  +  6x  )  cos(0  +  $)  +  k(x  +  Sx  )  cos0  (78) 

n  n  n  b  n  n  b 

The  total  phase  shift  when  pointing  at  angle  (0fe  +  )  is  the  sum  of  the 

phase  determined  by  focusing  on  the  beacon  (the  phase  in  the  exponent 
in  (77))  plus  the  increment  in  (78);  i.e., 

A4>  *  -kx  cos9,  -  f '  +  ¥ 
n  n  b  n  n 

“  -kx  cos0,  -  <f>'  -  k(x  +6x  )  cos (6,+  $)  +  k(x  +  6x  )  cos0. 

non  nn  b  n  n  b 

*  -kx  cos(0,  +  $)  -  i P'  ~  k5x  cos(e,+  4>)  +  k6x  cos6.  (79) 

ti  d  n  n  b  n  b 

In  (79)  -kx  cos(0,  +  4>)  is  the  correct  phase  shift  and  -<J>'  -k6x  cos 
n  b  r  n  n 

(9.  +  I)  +  kix  cos0  is  error.  L#t  i (/„  denote  the  error  of  the  nth 
b  n  b  en 

element. 
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(80) 


iji  »  ~<p'  -  kfix  cos (0t  +  $)  +  kfix  cos6. 

en  n  n  b  n  b 


-ib'  -  k<5x  cos8.cos$  +  k6x  sin6.  sin<J>  +  k6x  cos8. 
nnb  nb  nb 


If  the  scanning  angle  <j>  is  small,  cos$  =  1  and  sin<J>  =  $.  (80)  is  then 

approximated  by 


+en  “  “*'n  +  k<S>«xnsin8b 


(81) 


4»en  is  a  function  of  9^.  If  the  array  is  being  used  for  near  broadside 

reception  8,  «  90°  and  sin  8  can  be  set  equal  to  unity.  If  6x  is  zero 
b  b  n 

2 

mean  variance  cr.  ,  is  also  zero  mean  with  variance 
6x  en 


en 


2  ^  ,A2  ,  2.  2 

a ,  ,  +  k  ip  sin  8,  a 

*  n  b  Sx 


(82) 


where  a  „  is  given  by  (76). 

*  n 

Example: 

We  assume  a  reference  dish  antenna  with,  (1)  a  30dB  mainlobe  to  sidelobe 

ratio,  (2)  signal  and  interference  of  equal  level,  both  20dB  above  the  noise 

of  each  element,  and  (3)  equal  noise  generated  in  the  array  elements  and 

the  reference  antenna  amplifers.  Thus 

a2  -  1000.  b2  -  1.  a2  -  B2  -  100c  2  -  100a  2 . 

n  r 


The  variance  due  to  estimation  is  given  by  (76) 
2 


V 


bV 


_  I  |  2  4 
n  2  a  a 


♦TT-TT1*  s2'-2!.!2  ♦  S  »Vl 

2  a  a 


1  ,T1  2,  ,2  2  .  T1  „2 1 ,  ,2  2 

+  ~;2~4  [2T  a  lal  °n  +  2T  6  |b|  an 
2  a  a 


.  T1  22. 
+  2T  °n  °r  1 


P-36 


The  total  variance  a.  is 

ven 


2  2  2 

+  k  ♦  sin  eba6x 


2 


with  sin  9,  *  1, 
b 


_1_ 

20 


and  a&  in  wavelengths 


2  ,  <2”>2 

2 

(20)Z 


o 


2 

5x 


*  5x10 


-jio° 


For  F  »  0.95e  ,  N  *  7,  and  Y  =  10,  the  standard  deviation  for  1  dB 

s  a 

loss  is  about  2°  or  0.04  radians.  Let  o  «  0.04^  ■  0.0016  then 

ven 

=  0.1X;  this  is  a  fairly  tight  tolerance  requirement. 

However,  for  F  *  0.95,  N  -  7  and  Y  -  10,  the  standard  deviation  for 

S  Q 

1  dB  loss  is  about  7°  or  0.14  radians.  The  tolerance  of  element  positions 
is,  in  this  case,  -  Q.44X;the  tolerance  demands  are  more  relaxed  in 
this  case. 

We  point  out  that  in  these  examples  the  dominant  contribution  to 
the  total  pointing  phase  error  is  the  error  induced  by  steering  from 
beacon  to  ground  station.  That  is,  if  the  values  chosen  for  reference 
antenna  sidelobe  levels  are,  in  fact,  realistic,  the  method  is  shown  by 
these  examples  to  result  in  adequate  estimates  of  the  beacon  pointing 
phases.  The  adequacy  of  the  final  estimates  of  pointing  phases  to  the 
ground  station  depends  now  mainly  on  the  scanning  range,  and  the  un- 
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certainty  in  element  position.  If  the  error  induced  by  the  scanning 

process  were  totally  eliminated  the  pointing  vector  would  be  left  with 

2-4  o 

an  error  variance  of  o  *  5x10  (a  standard  deviation  of  about  1  ) 

en  -110° 

which,  even  with  Fg  -  0.95e  J  ,  causes  a  loss  in  SINR  is  negligible. 
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random  variable 
imaginary  part  of  y 
real  part  of  y 

position  error  of  the  array  element 

total  phase  shift  of  the  nC^  array  element  when  pointing  at  angle 
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beacon  arrival  angle 
inter ferer  arrival  angle 
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real  part  of  A 
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noise  power  at  array  elements 
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variance  of  the  error  of  array  element  locations 
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variance  of 

en 
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duration  of  signal  chips 


inverse  of  the  bandwidth  of  $  (f) 

n 

inverse  of  the  bandwidth  of  $r(f) 
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List  of  Symbols  -  Appendix  D 
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scanning  angle 

electrical  phase  angle  of  the  beacon  signal  at  the  origin 

electrical  phase  angle  of  the  interferer  at  the  origin 

(f)  power  spectral  density  of  noise  at  array  elements 

(f)  power  spectral  density  of  noise  at  the  dish  antenna 
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I.  INTRODUCTION 


Adaptive  Arrays  for  Communication  based  on  the  LMS  algorithm  of  Widrow  [1], 
but  utilizing  a  self-generated  reference,  have  been  subjects  of  study  in  the  past 
[2].  These  methods  are  based  on  the  exploitation  of  a  distinguishing  character¬ 
istic  of  the  desired  signal  which  enables  it  to  be  separated  from  unwanted  signals. 
Typically,  the  desired  signal  is  viewed  as  being  spectrum  spread  by  a  known  code 
so  that  with  despreading  a  moderately  clean  replica  of  the  desired  signal  is  de¬ 
rived  from  the  array  output.  With  this  as  the  reference,  array  focusing  toward 
the  desired  signal  is  improved,  in  turn  improving  the  reference.  The  scheme  pro¬ 
posed  and  analyzed  heretofore  is  shown  in  Figure  1  [3],  and  Is  here  called  the 
"LMS  adaptive  array  with  coded  reference  signal  loop."  Each  array  element  is 
s to  be  weighted,  the  weight  being  determined  via  an  LMS  loop  driven  by  the 
correlation  between  element  output  and  array  output  with  reference  subtracted. 

The  limiter  functions,  in  part,  to  prevent  the  weights  from  converging  to  zero. 

A  difficulty  observed  with  this  arrangement  Is  that  Incidental  phase  shift  In 
the  band-pass  circuitry  of  the  reference  generating  loop  will  cause  the  weights 
to  cycle  indefinitely  with  cycling  amplitude  range  and  frequency  depending  on 
the  main  feedback  loop  gain  and  on  the  magnitude  of  the  reference  loop  phase 
shift  [4,5].  Methods  of  adaptive  compensation  to  overcome  this  effect  have  been 
pursued  [6] ;  non-adaptive  correction  is  not  apt  to  work  because  of  uncertainty 
of  the  center  frequency  caused  by  doppler  shifts  and  oscillator  frequency 
errors. 

An  alternative  circuit  in  which  one  element  is  left  unweighted  and  the 
limiter  omitted  is  shown  in  Figure  2.  To  distinguish  this  from  the  previous 
circuit  we  refer  to  this  as  the  "LMS  interference  canceller  with  coded  reference 
signal  loop". 
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Constraining  the  weight  of  one  element  to  be  constant,  the  all  zero  array  weight 
condition  cannot  arise  making  the  limiter  unnecessary  on  this  account.  Mow,  how¬ 
ever,  operation  will  be  affected  by  both  phase  and  amplitude  differences  between 
reference  and  signal  component  of  the  array  output.  The  analysis  presented  below 
Is  Intended  to  reveal  the  properties  of  this  circuit  with  particular  attention  to 
the  weight  convergence  and  the  output  slgnal-to-noise  ration  (SNR).  It  is  shown 
chat  under  simple  attainable  conditions  the  weights  do  converge  to  a  constant  final 
value.  The  configuration  of  Figure  2,  however,  generates  a  variable  SNR  at  the 
summer  output,  the  SNR  possibly  being  as  low  as  that  of  a  single  array  element 
depending  on  the  reference  loop  gain  and  phase  shift.  Because  this  may  be  trouble¬ 
some  to  the  reference  generator  a  variant  of  this  circuit,  in  which  the  signal  to 
the  reference  generator  is  the  sum  of  weighted  clement  outputs  only  (see  Figure  6), 

is  also  analyzed.  The  SNR  at  the  reference  generator  input  is,  in  this  case, 
the  full  value  expected  out  of  (N-l)  coherently  combined  elements.  Finally,  the 
SNR  at  the  point  from  which  signal  output  is  to  drawn  is  obtained. 


II.  N-Element  LMS  Interference  Canceler 

An  N-element  interference  canceler  is  shown  in  Fig.  3.  The  signal  from  each 
antenna  auxiliary  element  yi(t),  i  ■  0,1... N-l,  is  passed  through  a  quadrature 
hybrid  forming  pairs  of  components  yAI(t)  -  y*(t)  and  y^qXt)  ■  yt(t),  i  ■  1,2,.. 
N-l ,  where  the  caret  stands  for  the  Hilbert  transform  (90*  phase  shift).  Each 
of  these  quadrature  signals  is  weighted  by  real  factors  Wjj(t)  and  U^q(t), 
respectively,  and  then  summed  together  with  the  signal  from  the  main  element 
yjj(t)  to  produce  the  actual  array  output  va(t).  The  difference  between  the 
array  output  and  the  reference  signal  d(t)  is  the  error  signal  e(t).  Utilizing 
eXt),  the  feedback  loops  adjust  the  weights  Wii(t)  and  Wiq(t)  so  that  the 
mean-square  error  |e(t)  |2  is  minimized. 

The  operations  performed  in  Fig.  3  are  represented  by  the  following:  The 
pairs  of  signals  yj/t)  and  y'i(t)  are  used  to  express  the  analytic  signal, 

*i(t)  -  tyA(t)  -  jy^t )]//!  (1) 

The  complex  weights  W^(t)  Vf£j(t)  +  jWiq(t)  i  ■  1,2. ..N-l  act  on  the  xj^t) 
to  generate  an  output  which,  in  analytic  signal  form,  is  given  by, 

v(t)  »  xT(t)W(t)  +  xN(t)  (2) 

where  W*(t)  -  (W^t),  W2(t),. .  .W^U) } ,  xt(t)  -  {x^t),.  ,.xN1(t)}  , 

T  stands  for  transpose,  and  Xfj(t)  ■  [yfl(t)  -  Jy^t)]/*^  is  the  analytic  signal 
representation  associated  with  the  signal  from  the  main  (unweighted)  element. 

The  actual  signal  at  the  output  of  the  array  is  given  by 

va(t)  -  Re{v(t)}  (3) 

and 

e(t)  -  d(t)  -  v(t),  (4) 

where  e(t)  and  d(t)  are  the  analytic  signal  representations  associated  with  the 
error  and  reference  signals,  respectively. 
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The  weight  are  assumed  to  be  adjusted  using  the  steepest  descent  algorithm 
[1];  that  is,  they  are  determined  by 

dW(t)/dt  -  -kAw  E[e(t)2]  (5) 

where  Ay  is  the  gradient  with  respect  to  W,  and  is  understood  to  be  a  complex 
vector  whose  components  are  the  gradients  with  respect  to  the  real  and  imaginary 
parts  of  W,  respectively,  and  k  is  the  main  feedback  loop  gain.  E(- )  denotes 
mathematical  expectation.  Following  the  derivation  in  [1],  (5)  reduces  to 

dW(t)/dt  -  2kE[x*(t)e(t)]  (6) 

where  the  asterisk  denote  complex  conjugate. 

Using  (6)  with  the  instantaneous  product  [x*(t)e(t)]  as  an  estimate  of 
the  corrcRoonding  expected  value  this  becomes  the  complex  LMS  algorithm  whose 
equivalent  network  for  one  antenna  element  is  displayed  in  Fig.  4.  Substituting 
for  e(t)  using  (2)  and  (4)  we  get  the  differential  equation  that  governs  the 
weight  vector  W(t) 

dW(t)/dt  +  2kRxM(t)  -  2^4  -  2kE[x*(t)xN(t)]  (7) 

where  Rx  ■  E[x*(t)xT(t)]  is  the  input  covariance  matrix  and  Rx<j  -  E[x*(t)d(t)] . 

In  obtaining  (7),  using  the  network  of  Fig.  4  as  a  basis,  an  assumption 
commonly  made  in  the  analysis  of  adaptive  arrays  is  used.  That  is,  that  the  signal 
and  weight  processes,  the  latter  now  being  a  random  process,  are  independent.  W(t) 
in  (7)  and  in  the  subsequent  work  should  therefore  be  viewed  as  an  expectation. 

For  the  case  of  a  single  CW  signal  of  amplitude  A  arriving  at  angle  to 
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broadside  we  have(t) 

x(t)  »  Agexp[ juictj//l  +  N(t)  (8) 

xN(t)  -  APNexpl jwct]//7  +  nN(t)  ($> 


where  PT  «  (e  ^  e  ^  *...  e  P(j  ■  q  ^°N  “  1 

ai  -  (ZHLi/XcJsln*  1  -  1...N-1  (10) 

A  is  Che  signal  amplitude, Xc  is  Che  free  space  wavelength  aC  frequency  uc,  the 
are  the  distances  between  the  ith  and  the  main  element,  NT(t)  *  {n^(t) . .nN_|(t)} 
is  the  noise  process  vector  and  nN(t)  is  noise  process  at  the  main  (unweighted) 
element.  The  noise  at  the  ith  element  is  written 

1^(0  -  0^(0  -  Jn1q(t)  (11) 

with 

ElniI2(t)]  -  E[n1Q2(t)l  -  on2/2.  (12) 

Hence 

E[ni(t)nj*(t)]  -  o^2  for  i  •  J,  and  every  i 

-  0  i  *  j  (13) 


The  covariance  matrix  of  the  element  output  is  therefore 


where  4 


5x  "  *  +  'll 

A2p*px/2. 


Therefore  (7)  can  be  written  as 
dW(t)/dt  +  2k(4  +  o„I)W(t) 


2k*xd  -  KA2P% 


(14) 


(15) 


(t)  We  assume  without  loss  of  generality  that  x(t)  contains  no  Interference.  In 
fact,  the  interference  tqrm  will  only  change  the  matrix  4. 


III.  lf-Element  Interference  Canceler  with  Reference  Signal  Loop 
a)  Solution  of  the  Equation  for  Expected  Array  Weights. 

When  the  array  output  v(t)  is  processed  to  produce  a  reference  signal  as  in 
Fig.  5,  then 

d(t)  -  {APTexp[juct]//I  W(t)  +  APNexp[juct)//T]ae”J*  (16) 

where  a  is  a  positive  constant  depending  on  loop  gain  and  filter  attenuation 
and  $  is  a  phase  shift  introduced  by  the  reference  loop  filter.  Noise  and 
interference  entering  the  reference  loop  will  be  substantially  attenuated  fol¬ 
lowing  the  filter  and  as  a  consequence  d(t)  is  represented  in  (16)  as  free  of 
these  influences. 

Using  the  definition  of  Rxcj  we  have 

5xd  "  E{x*(t)(^expIjuct]/^W(t)+A//IPNexpljuct])ae"J4' 

-  I*W(t)  +  A2P*PN/2]ae"J,>  .  (17) 

Thus  (15)  becomes 

dW(t)/dt  +  2k [ ( l-ae“ J ♦) 4+cj2  i ] w ( t )-  -kA2(l-ae-J<i>)P*PM.  (18) 

*•  -  ^  n 

The  matrix  *  is  Hermltian  having  rank  equal  to  one. 

Therefore  there  exists  a  constant  unitary  matrix  Q  such  that  Q'*Q  ■  A  (the  prime 
denotes  transpose  conjugate)  where  A  is  a  diagonal  matrix  with  only  one  non-zero 
element  X^-l* 

Premultiplying  (18)  by  (}'  we  get 

dr(t)/dt  +  2k[(l-ae"J*)A  +  ahlr(t)  -  -kA2( l-ae_J*)Q'p*PN 

-  *  (19) 

where 

r(t)  -  Q’w(t)  (19e) 

aw  «*  w 

(19)  is  obtained  also  using  the  fact  that  q'q  ■  I.  It  can  be  easily  shown 
that  §'p*  is  a  complex  vector  with  only  one  non-zero  entry  which  equals 
(/2Ajj_j/A]eJ$,  where  6  is  some  arbitrary  constant.  We  will  choose  6-0. 
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With  rT(t)  -  lYt(t) . .  Yji_i(t)]  (19)  can  be  written  [5] 


dYA(t)/dt  +  2ko|Y1(t)  -  0  i  -  1,2 . N-2 

(20) 

4yn_i  ( t  )/dt+2k[  ( l-ae“ J ♦)  XN_i+o|l  Y(j— 

-  -2kA/XN_1/2(l-ae~j4) 

(21) 

Note  that  in  (20)  we  have  used  the  fact  that  ■  0  for  1  -  l,2...N-2. 

The  solution  to  (20)  Is  given  by 

Yi(t)  -  Yi(0)  exp[-2kont]  1  -  1,2. ..N-2  (22) 

(21)  is  solved  by  decomposing  it  into  its  real  and  imaginary  parts  which 
are  given  by 

dY§-i (c )/dt+2k( ( l-acos^)X[I.1+o^]  y^-x ( c )+2kasin$XN_iYN-i ) 


»  -2kA/XN»i/2[l-acos$]  (23) 

dYN-i  ( t  )/dt-2kasin$XN_iY§_i  ( t  )+2kf  ( l-acos*)  XN_x+  °n  Wn-x  (  c  ) 

•  -2kA/Xfj_x/i Easing]  (24) 


where  yn_x(0  -  Yg_x(t)  +  jY^U) 

These  two  equations  can  be  written  in  matrix  form  as. 


dYfl^O/dt 

dYN-l(t)/dtj 


+  2kXN-l 


-2kA/X[j-l/2 


(25) 


where 


c  ■  On^N-i  (26) 

Using  a  dlagonalization  transformation  (see  Appendix  A-l)  we  find  the  solution 
to  (25)  to  be  (see  Appendix  A-l) 

YjS-i(t)  ■{(Yj5-x^°J-cr^)l  cos(2kXN_j  asin$)t 

“1yn-i (0)”Cj(^)]  sin(2kXN-iasin^)t) 
exp.  l-2kX(j-x  ( l-acos$+c)t  1+Cr(^) 
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(27) 


»  {(y^-I*0*-0!^1  cosUkX,,.!  asin$)t 

+  Iy5-1  C°>“Cr  (♦>  1  Sin(2kxfl>1  asin$)t} 
exp[-2kXu-l<l-aco8^+c)t]  +  C*($) 

2  2 

(l-acos»)(l-acos*t+c)  -  a  sin  6 
(l-acos$+  )2  +‘  a2sin2$  (28) 


where 


Cr(#)  -  ~ 


/5xN-1 


ct(*) 


aalnj(2  -  2acosit  +  c) 
(l-acos$+  )2  +  a2sin2$ 


(29) 


b)  Stability  Condition 

Stability  is  viewed  as  requiring  the  expected  value  of  the  weight  vector 
U(t)  to  converge  to  a  constant  steady  state  value.  This,  in  turn,  requires 
the  vector  r  given  by  (19a)  to  approache  a  constant.  The  components  of  r  ,  yi> 
i  “  1,2...N-1,  must  therefore  converge.  The  first  (N-2)  components  given  by  (22) 
will  converge  for  k  >  0.  For  Y§_i(t)  and  Y(j_j(t)  to  converge  to  the 
final  values  Cr($)  and  Cj(+)  respectively,  it  is  necessary  and  suf¬ 
ficient  that 

1  -  acos  4  >  -c 

2 

be  satisfied.  Furthermore,  one  can  show  that  Xfj_i  ■  (N-l)A  /2  and  there¬ 
fore  using  (26)  this  condition  becomes 

1  -  acos  *  >  -2o2/A2(N-l)  (30) 

For  this  to  be  satisfied  for  any  signal-to-noise  ratio  (A2/2o„)  it 
is  sufficient  to  require  that 

acos  $  <  1  (31) 


In  particular  if  we  choose  a<l  the  system  converges  to  the  steady  state  re¬ 
gardless  of  the  input  SNR  and/or  the  amount  of  phase  shift  in  the  reference  loop. 


e)  Array  Output  SNR 


The  complex  envelope  of  the  desired  signal  at  the  array  output  Is  given  by 
s(t)  -  (A//I)  tP^U)  +  PN] 

-  (a/ /I)  [pTgr(t)  +  PN] 

By  using  (20)  we  get  in  the  steady  state 
lim  s(t)  -  (A//2)  tPTq  Y  (-)  +  P  I 

t-M.  -  1  N-l  N 

where  qt  is  the  last  column  of  the  matrix  Q.  By  definition  (j  (Q  •  (A^/2)Q  P  ?TQ  -  A, 

and 

Ola) 


therefore  (A//7)Ptq^  *- 


11m  s(t)  "  yn-1^**^  /aSTT  + 

t-M* 

If  the  stability  condition  is  satisfied  then  lim  s(t)  exist  and  the  output  desired 

t-M* 

signal  power  S0  is  given  by 

So  -  IYn-^-)  +  A//2|2  (31b) 

With  YN-l(®)  "  Cr($)  +  jCi(^),  and  using  (28)  and  (29) ,  we  obtain 

SG  -  A2[(l  +  aJ(*  }  +  a12(^)]/2  (32) 

where 

03) 

(34) 


ar($)  *  *"Cr(^)/2XjJ_j/A 


a^(^)  **  ~C£(^)/2Xjj_j/A 
Also  the  noise  power  is  given  by 

■  on^lA2/2ljj_j(ar2(^)  +  ai2(p)  +  1] 

Since  XN_j  ■  (N-l)A2/2,  we  have 

Nq  -  on2[(l/N-l)(or2(^)  +  o^*)  +  1] 

From  (32)  and  (35)  we  finally  write  for  the  output  SNR 


(35) 


S  _  i2/m  (  \  .  (1^  2($)  +  at2(4)  +  2ar($)  +  l 


o  -  A  (N-l)  r 

k 


orz(^>  +  a^O)  +  N-l 

.2 


A^(N-l) _ (2-2acoso  +  c)*  _ 

2onz  N[ ( l-acos$+c)2  +  azsinz$)-c(2-2acos$+c) ] 


(36) 
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The  last  step  Is  obtained  In  (E.13)  of  Appendix  A-2.  Note  that  in  contrast  to 
the  LMS  adaptive  array  with  coded  reference  loop  as  treated  by  Compton  and  DICarlo, 
here  the  output  signal-to-noise  power  ratio  depends  on  the  reference  loop  phase 
shift  and  gain.  This  ratio  is  less  than  A2(N-l)/2on2  for  most  practical 
cases  (see  Appendix  A-3)  if  we  choose  a  to  be  less  than  one.  The  output  SNR  may 
even  fall  to  the  output  SNR  of  a  single  element,  A2/2on2  (see  Appendix  A-4). 

From  these  results  it  is  evident  that  in  the  network  of  Fig.  2.  The  array 
element  outputs  are  not  being  maximally  combined  to  give  the  expected  array  gain 

of  an  N  element  array.  In  particular,  the  variability  of  the  SNR  at  this  point  is 
apt  to  affect  the  performance  of  the  reference  generating  loop.  However,  we  now 
show  that  utilizing  only  the  weighted  element  outputs  as  reference  loop  input  one 
gets  an  SNR  close  to  the  maximum  and  not  dependent  on  the  parameters  of  the  refer** 
ence  generating  loop. 

The  sum  of  the  weighted  array  element  outputs  (auxiliary  output)  alone  Is 
vau(t)  -  xT(t)W(t) 
and  the  signal  component  is 

sau(t)  -  APTW(t)//I 

This  can  be  expressed  in  a  form  similar  to  (31a)  with  the  component  contributed  by 
the  unweighted  element  deleted.  Thus 
sau(c)  *  YN-l("^*N-l» 

t-M» 

and  the  power  of  the  desired  signal  at  the  auxiliary  output  is 

soau  *  Iyn-1<*>I  Vl  <V) 

The  corresponding  noise  power  is  similar  to  (34a)  with  the  contribution  by  the 
unweighted  element  deleted.  Thus 

»» *  »„2'  W>i2  <“> 

and  the  signal-to-noise  power  ratio 

<Wa  ‘  W°n2  “  A2<N-l)/2on2  <»> 
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The  array  gain  of  the  (N-l)  unweighted  elements  is  therefore  realized;  that 
is,  the  weights  adjust  themselves  to  coherently  combine  the  output  of  the  weighted 
elements  and  to  provide  maximum  SNR.  However,  the  resultant  weighted  output  is 
not  coherently  combined  with  the  output  of  the  unweighted  element;  the  SNR  of  the 
combination  depends  on  the  parameters  of  the  reference  generating  loop  and  has  a 
value  generally  well  below  that  of  the  combined  unweighted  elements. 

He  now  propose  a  modified  circuit  which  takes  advantage  of  the  foregoing  to 
provide  constant  SNR  to  the  reference  loop,  but  which  also  combines  the  weighted 
and  unweighted  element  outputs  to  give  a  good  final  output  SNR. 

IV.  Alternative  Arrangement  for  Extracting  the  Reference  Signal 
Signal 

The  system  arrangement  depicted  in  Fig.  6.  is  proposed  to  overcome  the 
problem  discussed  above.  Here  the  reference  signal  d;(t)  is  given  by 
dj(t)  -  a  xT(t)W(t)e“j* 

where  a  and  4  are  gain  and  phase  shift  in  the  reference  loop.  Consequently, 

Jgxdl  “  E  tx*(t)  xT(t)W(t)]  ae-J* 

-  «W(t)ae-j*  (40) 

Instead  of  (18)  the  weight  vector  W(t)  will  be  governed  by 
dW(t)/dt  +  2k((l-ae'j*)«  +  o„2I}W(t)  -  -kA2P*P„ 

where,  P  and  P^  and  4  were  defined  in  connection  with  (8),  (9)  and  (14). 

Following  the  same  procedure  as  used  earlier  two  equations  in  f Jj_i (t)  and 
Y^_j(t)  similar  to  (25)  are  obtained  with  the  right  hand  side  being 
-2kA/AN_j/2^j.  The  solution  to  these  equations  is  given  by  (27)  with 


E-ll 


cri<*>  -  - 


A _  (1  -  acosa  +  c) _ 

Xkj_t  (1  -  acos$  +  c)2  +  a2sin2* 


(42) 


The  stability  condition  remains  exactly  the  same  (eqs.  (30)  and  (31)).  Even 
though  the  final  weights  different  from  what  they  were  in  the  previous  system  set 
up,  the  output  desired  signal  and  noise  powers  remain  as  they  are  given  by  (37) 
and  (38),  respectively.  Therefore  Che  SNR  of  the  weighted  element  output,  which 
is  also  SNR  fed  through  the  reference  loop,  equals  a2(N-1)/2o2. 
a)  Output  SNR 

Signal  output  in  the  configuration  of  Fig.  6  is  taken  from  the  point  where 
weighted  and  unweighted  elements  are  summed,  as  shown. 

The  desired  signal  and  noise  power  output  at  this  point  are  given  by  (32)  and 


(35)  respectively,  with  arl(^)  and  oil($)  used  instead  of  q r($)  and 
a^(^),  and  with 

ari($)  *  “  Cr i ($)/2XjjJj/ A  (43) 

On ($)  m  "  Cu (♦)^l[|_j / A  (44) 

Crj($)  and  Cjjfifc)  are  given  by  (41)  and  (42).  Using  appendix  (A-5)  we  write 


for  the  output  SNR 

,so,  _  A2(N-1)  ,arl2<$>  +  °il2(*>  +  2arl<4>>  +  1 

( — » - - -  [ - - - - - 

N0  2on2  ®rl2(f)  +  Oil2  (41)  +  N-l 


-  A2 (N-l)  r 

- 9 -  1 

2an2 


2  2  2 

(2  -  acos$  +  c)  +  a  sin  a 
(N-l )( ( 1  -  acos$+c)2  +  a2sin2$) 


(45) 


As  in  the  case  of  the  configuration  of  Fig.  2,  with  the  configuration  of  Fig.  6 
the  output  SNR  depends  on  the  reference  loop  phase  shift  and  gain.  This  ratio 
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exhibits  a  maximum  Value  when  the  phase  shift  is  the  smallest,  $  ■  0  (provided 
that  N  >  3,  or  else  for  N  •  2,  A^/2on^  <  2.5  is  necessary).  If  furthermore 
a  «  1  then  for  most  practical  cases  the  output  SNR  is  approximately  equal  to 
2  2 

its  optimal  value  A  (N-l)/2on  (see  Appendix  A-6).  Therefore,  from  the 
point  of  view  of  output  SNR  the  configuration  of  Fig.  6  is  preferable  to  that 
of  Fig.  2. 

V.  Conclusion 

He  have  here  reported  an  analysis  of  the  properties  of  an  adaptive  array 
based  on  the  LMS  algorithm  with  a  self-generated  reference.  Circuits  hereto¬ 
fore  proposed,  analyzed,  and  tested  had  been  reported  to  suffer  oscillatory 
weight  Instability  when  the  reference  generating  loop  shifts  the  phase  of  the 
signal  component  of  the  array  output.  The  scheme  examined  here  is  a  variant 
of  the  earlier  scheme  in  which  one  array  element  is  left  unweighted  and  the 
reference  generating  loop  is  operated  without  a  limiter.  As  a  consequence  the 
signal  output  of  the  reference  generating  loop  is  not  fixed  in  amplitude  and 
may  differ  from  the  signal  at  the  loop  input  in  both  amplitude  and  phase. 

It  was  found  that  in  this  mode  of  operation,  stable  nonoscillatory  weights 
are  obtained  in  the  steady  state  provided  certain  mild  conditions  are  satisfied  on 
the  gain  and  phase  shift  of  the  reference  generating  loop.  Two  alternative  ways 
of  driving  the  reference  generating  loop  were  examined.  In  the  first,  the  sum  of 
all  array  element  outputs  were  used  to  drive  the  loop.  Here  it  was  found  that  the 
SNR  entering  the  loop  was  highly  variable  being,  under  some  conditions,  nearly  equal 
to  the  SNR  of  a  single  array  element.  This  result  was  not  viewed  as  a  desirable  one 
and  an  alternative  in  which  the  reference  generating  loop  was  driven  by  the  weighted 
elements  alone  was  then  studied.  Here  a  constant  high  SNR  to  the  loop  is  obtained 
and,  under  reasonable  conditions  of  loop  gain  and  phase  shift,  a  high  (though  mod¬ 
erately  variable)  signal  output  SNR. 
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Arrays  of  the  sort  here  discussed  are  used  where  adaptive  interference  can¬ 
celling  is  sought.  The  analysis  reported  here  as  well  as  that  reported  earlier  on 
other  circuits  with  a  self-generated  reference  signal  deal  only  with  desired  signal 
plus  noise.  Clearly,  all  candidate  schemes  ought  to  also  be  examined  with  inter¬ 
ference  present  and  the  output  ratio  of  signal  to  Interference  plus  noise  deter¬ 
mined.  This  will  undoubtedly  be  done.  Aa  an  interim  measure  we  have  also 
determined  the  component  of  signal  power  in  the  residue  entering  the  array 
control  processor.  The  details  are  not  reported  here  but  the  analysis  showed 
that  the  first  circuit  resulted  in  domewhac  less  signal  power  at  this  point 
than  the  second  circuit.  The  reason  for  concentrating  on  this  quantity  is 
the  expectation  chat  the  scheme  generating  a  smaller  signal  power  residue  will 
be  a  more  effective  interference  canceller;  the  desired  signal  would  then  not 
mask  the  Interference  and  it  would  be  the  latter  which  mainly  influences  the 
weights  generated. 


Appendices 

A-l  By  using  Che  following  information 


Pl<t> 

m 

l 

j 

tn-i<c> 

p2(t) 

j 

l 

»  m 

(25)  becomes 


dp^O/dt" 

l-aaos$+c-jasln$  0 

Pl(t) 

+2kXN_i 

dp2(t)/dt 

0  l-acos$+c+jasin$ 

a 

P2^c^ 

-  -2kA^ 


l-acos$+jasin$ 


(E.2) 


[j(l-acos$)+asin$  J 

(E.2)  yields  two  disjoint  differential  equation  in  pj/t)  and  p2(t).  The 
solution  of  which  is  given  by 

Pj(t)  -IpjW-CjC*)]  exp [-2kXjj_ j ( 1-acos $+c- jasin$) t ] 

+  Cj(t)  (E. 3) 


p2(t)  *Ip2(°)"c2^*^  exp{-2kAjj_|( l-acos^+c+Jasin^)t) 


where 


+  C2<^) 

A  1  -  acos<t  +  jaslna 

✓2An— i  1  -  acos$  +  c  -  jasin$ 


(E.  4) 


(E.  5) 


C2($) 


A  j(l  -  acosit)  +  jasinj 

/2Xn-1  1  "  »cos$  +  c  +  jasin$ 


(E.6) 
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Transforming  back  by  using  the  transformation  of  (E.l)  we  get 

-  {[?§_!«»  -  Cr($)J  cosUklH.jasin^t 
~  Iy^-i(O)  -  0^$)]  sin(2kXN_iasin$)t> 
expI-2kXN-l(l  -  acos$+c)t]  +  Cr($) 

Cr(ip)  is  given  by  (28). 

A-2  From  (33)  and  (34)  together  with  (28)  and  (29)  we  write 

.  b,2($)  -  b22 ( 4»)  +  cb,($) 

ar(^;  ■  _ _ _ t _ x  - 

t>l^($)  +  b2  ^  ( ^)  +  2cbj($)  +  c2 


2bi($)b2($)  +  cb2($) 

oi(4>)  ”  — = - = - 

bi  ($)  +  b2z($)  +  2cb^(<^)  +  c 


where 


bj(^) 

b2($) 

Thus 

or2(4>)  +  ai2(ij») 


l-acos<{> 

asin$ 


(bI2(»)  -  bz2(»))2  +  c2(b12(»)  +  b22(<fr)) 
bi2(^)  +  bz2($)  +  2cbj;($)  +  c2 


2cb1(^)(b12(^)  +  bz2(«)) 

bj2($)  +  bz  2  ( <(>)  +  2cbj(.$)  +  c2 


Also 
2ar($)  + 


b12(^)  +  b22(^) 

bi2(^)  +  b2Z($)  +  2cbj(^)  +  c2 

3bt2(»)  ~  b22(»)  +  4cb^($)  +  c2 
bj2(^)  +  b2^(<>)  +  2cb}(t)  + 


so  that 

ar2($)  +  a^2(^)  +  1 


(2bt(»)  +  c)2 _ 

b!2(4>)  +  b22($)  +  2c  b !  ( )  +  c2 


(E.7) 


(E.8) 

(E.  9) 

(E. 10) 

(2.11) 


(E. 12) 
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Finally 

ar^($)  +  a^($)  +  2af($)  +  1 
ar2($)  +  a^2($)  +  N-l  ”” 

(2b1(^)  +  c)2 

N(bj2(<fr)  +  b22($)  +  2cbj($)  +  c2)  -  c(2b}($)  +  c) 


_ (2  -  2acos<fe  +  c) _ 

N((l-acos$+c)2  +  a2sin$)  -  c(2-acos$+c) 


where  In  Che  lasC  step  we  used  (E.10)  and  (E.ll). 


A- 3  To  examine  Che  dependence  of  S0/No  on  Che  reference  loop  phase  shlfc,  of 
Che  clrcuic  of  Fig.  2  we  rewrice  (36) 


(E. 13) 


£o(a,<j>)  -  1  (_ 


(2  -  2acos<fe  +  c)J 


N((l-acos$+c)z  +  a4sinz$)  -  c(2-2acos$  +  c) 


(E. 14) 


where  c  ■  2cn2/A2(N-l) 

SeCCing  c3S0/N0(a,$)/3$  *0,  yields  Che  equaClon 

2(2  -  2acos<$r*-c)2asln$[N((  l-acos<|F+-c)2  +  a2sin2<j,)  -  c(2  -  2acos$+c)] 

-  2(2  -  2acos$+c)2[N(2( 1  -  acos$+c)asin$  +  2a2sinij>cos$)  -2acsin$] 

“  0  (E. 15) 

The  sCadonary  poincs,  cherefore  occur  under  one  or  more  of  Che  following 


condidons. 


1)  2-2aco8$  +  c  +  0 


(E. 16) 


or,  acos#  ■  (c+2)/2.  This  condldon  cannoC  hold  for  a  <  1  and  c  *  0 


2)  2a  sin$  -  0 

or,  ^  "  0  since  a  *  0 


(E. 17) 


3)  2N((l“acos^+c)  +  a2sin2($))  -  2c(2-2acos$+c) 

-  (2-2acos$+c)  (N+c(N-l)) 

which  Implies 

2-2acos$+c  -  2N(l-a2+c)/(N+cN-c)  (E. 18) 

The  values  of  S0/N0  at  these  points  are,  respectively 

1)  S0/N0  -  0 

2)  S0/N0  -  1  - (2-2a+c)2 - 

c  N(l-a+c)2  -  c(2-2a+c) 


-  A2(N-1)  _ (2-2a+c  )2 _  (E.19) 

2  0n2  N(l-a+c)  -  c(2-2a+c) 

It  is  possible  to  show  that  if  we  choose  a  <  1  for  this  case;  (i.e.,  $  *0), 
then  the  second  factor  of  (E. 19)  decreases  with  c.  Therefore  S0/N0  is  .upper 
bounded  by 


VNo 


<  A2(N-1) 
2°n2 


4(l-a)2 

N(l-a)2 


-  A2  (N-l  +  (4-NHN-l))  (E.20) 

The  upper  bound  is  larger  than  A2(N-l)/2on2  only  if  N  <  4.  It  equals 
A2N/2on2  when  N  ■  2,  and  approaches  A2/2o2.4  for  a  very  l^rge  N 

3)  S./N  -  A2^H-1)  [4N(l-a2+c)|  (E.21) 

2on2  (N+cN-c)2 


Again,  it  is  possible  to  show  that  if  we  choose  a  <  1  then  the  second  factor  of 
(E.21)  decreases  with  c  and  hence 

So/No  <  tA2(N-l)/2an2]4(l-a2)/N 

The  upper  bound  can  be  greater  than  A2(N-l)/2on2  only  when  N  <  4,  and 

equals  A2N/2on2  when  N  <  2.  This  upper  bound  is  attained  when  c  is 

very  small.  To  conclude  we  notice  that  the  output  signal-to-noise  ratio  depends 
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on  t  and  It  Is  less  than  A2(N-l)/2on2  for  most  cases  of  practical  Interest. 


A-4  To  examine  the  dependence  of  S0/N0  on  the  reference  loop  gain  In  the  cir¬ 
cuit  of  Fig.  2  use  (E.  14)  to  write 

a(S0/N0)/3a  ■  2[2-2acos$+c)(-ccos$(N+Nc-c-2) 

+  2a(cos2*(N+Nc-c)-N-Nc/2J 

Stationary  points  occur  under  one  of  the  following  conditions 

1)  (2-2acos$+c)  -  0  (E.22) 

or  a  cos^  -  ( c+2 ) / 2 .  Again,  this  is  not  valid  for  a  <  1  and  c  #  0 

(E.23) 


2)  a  -  -ccosi»[(N-2)+c(N-l)] _ 

l(l-cos2t)(N+cN/2)+(c-cN/2) 

cos2tl 


This  point  is  a  maximum  or  minimum  depending  on  whether  cos2t  is  less  or  greater 
than  (M+Nc/2)/(IH-Nc-c)  respectively.  But 

(E.24) 


1/2  <  N+Nc/2  <  l 
N+Nc-C 


for  N  >  2  and  any  c.  Therefore  for  cost  “  l»  the  minimum  point  occurs  at 
a  ■  2  if  we  assume  c  is  small.  From  (E.14)  for  t  *  0  and  a  **1  S0/N0  «  A^/2 on^, 
while  for  t  -  0  a  *  0  Sq/No  -  A2(N-l)/2onN. 

A-5  Using  (43)  and  (44)  together  with  (41)  and  (42)  we  write 


bj(t)  +  c 

a ri<$)  "  ■  — 

bj^(t)  +  t>2  ^  ( t )  +  2cbj(t)  + 


ail 22  2 

bj  ( t)  +  b2  (t)  +  2cbj^(t)  +  c 


where  bj(t)  and  b2(t)  were  defined  in  (E.10)  and  (E.ll).  As  in  Appendix  A-2 
we  now  have 

1 


o?i (t)  +  d?i (t)" 


bi^(t)  +  b22(t)  +  2cbi(t)  +  c2 
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Also 


20^1  ($)  +  1 


(bj($)+c)(bj($)+c+2)  +  b2($) 
bi2(^)  +  b22($)  +  2cb^  C  <J>)  +  c2 


so  Chat 


Finally 


a|i(^)  +  a|j  +  2orl(^)  +  1 

(b^*)  +  c+1)2  +  bf($) 

m 

(bj($)  +  b|($)  +  2cbj($)  +  c2 

a2^*)  +  aii  +  2orlC^)  +  1‘ 
o^(^)  +  a^(^)  +  N-l  ” 

Cbi < ♦)  +  c+1)2  +  b2($) 

(N“l) [b2($)  +  bg<$)  +  2cb1(*)+c2]  +  1 

(2-2acos$+c)2  +  a2sin2<j, 

“  . _ _ _ - 

(N-l){(l-acos$+c)2  +  a2sln2$]  +  1 


A-6  To  examine  che  dependence  of  (S0/N0)j  on  the  reference  loop  phase 
shift  $  of  the  circuit  of  Fig.  6,  we  write  (45) 

(Wit4**5  "  -  (<2“*co»*tc>2  *  a2s^"2» _ 

c  (N-l){(l-acos$+c)2  +  a2sin2$]  +  1 
Hence  c3(S0/NQ) ^(a,^)/3<(i  -  0  Implies  that 

2(c+2)  asin$[(N-l)(( l-acos$+c)2  +  a2sin2$)  +  1] 

-(c+l)(N-l)asln$[(2-acos<fr+c)2  +  a2sln2$J  «  .0 
After  some  manipulation  one  can  show  that  (E.25)is  equivalent  to 
2asln^[(N-l)a2  -  N(c+l)(c+2)  +  (c+2)2]  -  0 
A  solution  to  (E.26)  is  given  by 
asin  $  ■  0 

which  implies  $  ■  0  since  a  *  0 
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(E.25) 


(E.25) 

(E.26) 

(E.27) 


For  Chis  value  Co  be  maximum  it  la  necessary  and  sufficient  to  have 

a2  <[N(c+l)(c+2)  -  (c+2)2]/(N-1)  (E.28) 

The  right  hand  term  of  (E.28)  is  Increasing  with  c*  To  make  this  term  greater 
than  or  equal  to  one  and  consequently  (E.28)  eo  be  satisfied  for  every  a  <  1,  it 
sufficient  to  require 

3  <  N  +  (N-2)(c2+3c)  +  c2  +  2c  (E.29) 


This  is  satisfied  for  every  c  if  N  >  3  or  when  N 
(or  A2/3oJ<2.5). 


The  value  of  (S0/N0)i  at  $  ■  0  is  given  by 

iSQ/VQ)l  -  1  [ - 1 

c  (N-l)( l-a+c)2+l 


2  provided  that  c>0.4 


If  c  Is  small  enough  then 


<S0/N0>i  -  A2g=U 


2o; 


[  (2-a)2  ] 

(N-l)(l-a)^+l 


(E.30) 


If  c  Is  very  large  then 

<So/No>l  ■  k2fZol  (E*31) 

However,  c  large  is  not  a  practical  alternative.  If  a  Is  close  to  unity  then 

(So/No)l  *  (A2(N-l)/2o2)l(l+c)2/(N-l)c2+l)] 

(N-l)c2  -  (2on/A2)2/(N-l)  which  Is  much  smaller  than  unity  for  reasonably 
large  N  and  practical  value  of  nolse-to-signal  power  ratio  of  a  single  array 
element.  Therefore  (S0/N0)l  is  approximately  equal  the  optimal  slgnal-to^nolse 
ratio  A2(N-l)/2o2 

To  conclude  we  notice  that  the  output  signal-to-noise  ratio  of  this  config¬ 
uration  also  depends  on  the  reference  signal  loop  parameters,  a  and  +.  However, 
for  0  and  a  »  1  the  output  signal-to-nolse  power  ratio  approaches  the  optimum 
value  A  (N-l)/2on  for  most  of  the  practical  cases. 
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Figure  1  The  LMS  Adoptive  Array  with  Coded  Reference 
Signal  Loop  for  Sp.  Sp.  Communication 
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Figure  2  The  LMS  Interference  Canceler  with  Coded  Reference 
Signal  Loop  for  Sp.  Sp.  Communications 
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Figure  3  Quadrature  Weighted  LMS  Interference  Canceler 
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Figure  4  Main  Feedback  Loop  for  LMS  Algorithm  Realization 


▼ 


Figure  6  The  LMS  Interference  Conceler  With  Different 

Arrangement  for  Extracting  the  Reference  Signal 
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List  of  Symbols  -  Appendix  E 
desired  signal  amplitude 

constant  dependant  on  the  limited  level  and  filter 
attenuation  of  the  reference  extraction  loop 

2 

defined  to  be  a  /\  . 

n 

reference  signal 

reference  signal  for  alternative  arrangement 
denotes  expected  values 
minimized  error  signal  for  LMS  algorithm 
identity  matrix 

gain  constant  for  LMS  algorithm 
number  of  elements 
noise  process  vector  (1  x  N) 
noise  power  at  array  output 
noise  power  in  v 

aw 

noise  process  at  main  element 

input  signal's  phase  vector  (1  x  N) 

phase  at  main  element 

unitary  diagonalizing  matrix 

input  signal  autocorrelation  matrix  (II  x  N) 

input-reference  signal  cross-correlation  matrix  (N  x  N) 

input-reference  signal  cross-correlation  matrix  for 
alternative  arrangement  (N  x  N) 

output  desired  signal  power 

desired  signal  at  array  output 

desired  signal  at  array  output 

signal  component  of  v 

(Continued) 
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V(t) 

V  (t) 

SL 

V  (t) 
aw 

wC(t) 


Wil(t) 


xC(t) 

y±(t) 

yn(t) 

yiQ(t) 

yN(t) 

“i 

rc(t) 


v^t) 

Yxr(t) 

y^oo 

A 

A 


2 

a 

n 


w 

c 


List  of  Symbols  -  Appendix  E 
(Continued) 

array  output 
real  part  of  v  (t) 

sum  of  the  weighted  array  element  outputs 

array  weight  vector  (1  x  N-l) 
weighting  factor  for  y^(t) 

weighting  factor  for  y^(t) 

the  system's  input  vector  (1  x  N) 
signal  at  1C^  antenna  element 
equivalent  to  y^Ct) 

Hilbert  transform  of  y^(t) 

signal  at  unweighted  main  element 

Incoming  phase  of  input  signal  at  the  ith  element 

Q’  v(t) 

uncoupled  complex  array  weight  vector  (1  x  N) 
elements  of  rC(t) 

real  part  of  the  ith  uncoupled  weight  Y^(t) 

imaginary  part  of  the  iC^  uncoupled  weight  Y^t) 
uncoupled  signal  autocorrelation  matrix  (N  x  N) 
free  space  wavelength  at  frequency 
the  only  nonzero  element  of  A 

noise  power  present  at  each  element 
arrival  angle  of  input  signal 
narrowband  center  frequency  of  input  signal 
denotes  magnitude 

(Continued) 
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List  of  Symbols  -  Appendix  E 
(Continued) 

denotes  Hilbert  transform 
*  denotes  conjugate 

t  denotes  transpose 

’  denotes  conjugate  transpose 


E-31 


Appendix  F 

EFFECT  OF  INTERFERENCE  ON  THE  BEHAVIOR  OF 
AN  LMS  ADAPTIVE  ARRAY  WITH  REFERENCE  SIGNAL  LOOP 

by 
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and 
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Electrical  and  Computer  Engineering  Department 
Philadelphia,  Pa.  19104 


Effect  of  Interference  on  the  Behavior  of 


an  LMS  Adaptive  Array  with  Reference  Signal  Loop 


Abstract 

A  phase  shift  in  the  reference  loop  extraction  scheme  of  adaptive  arrays 
utilizing  the  Widrow  LMS  algorithm  has  been  shown  to  cause  cycling  of  the  array 
weights  in  an  Interference  free  environment.  Under  the  same  conditions  a  com¬ 
pensation  scheme  proposed  by  Bar-Ness  [3)  was  shown  to  force  the  weights  to 
converge  to  a  steady-state  value.  This  paper  investigates  the  compensated 
scheme  in  a  multiple  and  single  Interference  environment.  For  the  case  of 
single  interference,  it  is  shown  that  the  array  weights  of  the  compensation 
scheme  will  converge  to  a  steady  state  constant  value  while  those  of  the  uncom¬ 
pensated  scheme  will  continue  to  oscillate. 

I.  Introduction: 

Recently  adaptive  arrays  have  been  attracting  much  attention.  One  reason 
is  that  these  arrays  can  be  used  to  null  an  undesired  signal  coming  from  a 
different  source  location  than  a  desired  signal.  By  doing  this,  adaptive  ar¬ 
rays  can  be  used  as  automatic  beam  steerers  as  well  as  Interference  cancellers. 

One  popular  scheme  used  in  adaptive  array  systems  utilizes  the  least  mean 
square  (LMS)  algorithm  proposed  by  Widrow  et  al  [1].  The  scheme,  however,  re¬ 
quires  a  reference  signal  which  correlates  with  the  desired  signal  to  be  ex¬ 
tracted.  Obviously,  the  reference  signal  cannot  be  available  apriorl  to  the 
receiver  and  must  somehow  be  extracted  from  the  incoming  inputs  to  the  system. 

One  way  of  extracting  a  reference  signal  is  by  using  a  narrowband  signal 
while  having  broadband  interference.  Then  by  use  of  a  narrowband  filter  the 
Interference  is  attenuated  while  the  desired  signal  passes  through  and  is  used 
as  a  reference  signal.  This  condition  can  be  met  by  use  of  spread  spectrum 
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methods.  By  this  method,  the  desired  signal's  spectrum  is  spread  by  multiply¬ 
ing  the  signal  by  a  pseudo-random  code  at  the  transmitter.  Then  at  the  receiver 
the  incoming  signals  are  multiplied  by  the  same  pseudo-random  code.  This 
process  despreads  the  desired  signal  making  it  narrowband  and  at  the  same  time 
spreads  the  Interference  making  it  broadband.  Thus  the  reference  signal  can 
be  extracted  by  use  of  a  narrow  band  filter. 

One  problem  with  this  scheme  was  first  shown  by  Compton  and  DICarlo  [2]. 

They  demonstrated  using  an  Interference  free  input  that  a  phase  shift  Introduced 
by  the  narrow  band  filter  in  the  reference  signal  extraction  loop  causes  the 
array  weights  to  cycle.  This  work  used  a  single  array  element  and,  DICarlo  [4] 
later  extended  the  analysis  to  an  N-element  array  and  was  able  to  attain  the  same 
results.  Here  a  basic  problem  lies  in  the  fact  that  the  reference  loop  phase 
shift  is  a  function  of  the  incoming  signals  center  frequency  and  cannot  be  known 
apriorl.  Thus,  it  cannot  be  compensated  for  all  Incoming  signals,  and  under  a 
wide  range  of  nonstationary  environments. 

With  the  weight  cycling  problem  in  mind,  Bar-Ness  [3J  added  an  adaptively 
controlled  complex  weight  into  the  reference  loop  (Fig. 5).  This  complex  weight 
was  shown  by  Bar-Ness  to  effectively  compensate  for  the  reference  loop  phase 
shift.  Bar-Ness  also  showed  in  his  analysis  that  with  an  interference  free 
input  this  extraction  scheme  causes  the  array  weights  to  converge  on  a  steady — 
state  value. 

It  is  the  purpose  of  this  paper  to  examine  the  effect  of  interference  on 
Bar-Ness's  results.  This  first  will  be  done  for  the  general  case  of  a  multiple 
Interference  environment,  then  we  consider  the  special  case  of  a  single  inter¬ 
ference. 

For  the  multiple  interference  case  it  is  shown  that  the  reference  loop 
phase  shift  is  compensated  for,  however,  a  general  condition  for  weight  conver- 
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gence  has  not  yet  been  found.  For  the  single  Interference  case  a  sufficient 
condition  for  convergence  is  found.  The  uncompensated  case  of  Compton's  array 
will  also  be  examined  and  it  is  concluded  that  the  existence  of  Interference 


will  still  result  in  weight  cycling. 

II.  Array  Element  Weight  Controlling  Equations: 

A  basic  schematic  of  the  narrowband  adaptive  array  processor  is  shown  In 
Fig.  1.  The  incoming  signal  to  each  of  the  N  antenna  elements  is  split  into 
its  in-phase  (Xxx)  and  quadrature  phase  (XXq)  components  giving  2N  Inputs  to  the 
system.  Each  input  is  weighted  by  a  real  factor  wxx  or  wXq  and  then  summed 
to  produce  the  array  output  v(t).  An  error  signal  e(t)  is  produced  by  taking 
the  difference  between  a  reference  signal  d(t)  and  the  array  output  v(t). 

The  weights  are  controlled  such  that  the  mean  square  of  the  error  (e^(t))  is 
minimized . 

Following  the  results  first  made  by  Wldrow  [1]  the  weights  are  controlled 
according  to  the  equation: 


dw± ( t ) 
dt 


-  k  1  [TOT] 

3wi 


(l) 


where  k  is  the  main  feedback  loop  gain,  and  the  overbar  stands  for  expected  value. 
Realizing  that  we  represent  the  array  output  analytically  be 

v(t)  -  Xc(t)  w(t)  (2) 

where : 

wc(t)  -  {wjU),  w2(t),...,  wn(t)}  (3) 

Wj[( t )  ■  Wii(t)  +  jwiq(t) 

X^t)  -  (XjU),  X2(t) . X„(t)}  (4) 

Xi<t)  -  xlx(t)  +  JXiQ(t) 
t  -  indicates  transpose 
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Using  (2)  we  can  define  the  error  signal  as 


j 


1 


* 


e(t)  -  d(t)  -  v(t)  -  d(t)  -  xc(t)w(t)  (5) 

Substituting  (5)  Into  (1)  and  evaluating  the  partial  derivative  we  get  from 
Widrow  [1] 

dw(t)/dt  -  2k  (x*(t)e(t))  (6) 

-  2k  (x*(t>d(t)  -  x*(t)xt(t)w(t))  (7) 


or 


dw(t)/dt  +  2k  Rx  w(t)  -  2k  Rjyj 

where: 


(8) 


R x  -  E[x*(t)  xt(t)]  (9) 

Rxd  -  Efx*(t)  d(t)]  (10) 

and  E[  ]  denotes  the  expected  value. 

Notice  that  Rx  is  the  Input  auto  correlation  matrix  and  Rxd  Is  the  Input- 

reference  signal  cross  correlation  matrix.  Also,  using  eq.  (6)  we  see  that 

the  feedback  network  of  Fig.  (2)  can  be  Implemented  as  a  weight  controller. 

For  the  case  of  one  continuous  wave  (CW)  signal  and  m-lnterference  signals 

arriving  at  angles  i|»  and  X(^(l«l,..m)  respectively  with  reference  to  broadside 

(Fig.  (3)),  we  can  represent  the  Input  vector  of  eq.  (4)  as 

__  m 

x(t)  -  A//2  Ps  +  l  Bj/ /2  ?i  +  N(t)  (11) 

1-1  1 

where,  Nc(t)  -  [nj(  t)  .ujC t) , . .  ,nN( t) ]  Is  the  noise  process.  A  Is  the  desired 
signal  amplitude.  Bj  is  the  amplitude  of  the  1^  Inter ferer.  And  assuming 
a  nonlinear  array 


e’jp,<l;)a.  .'J“2 


(12) 

(13) 


where  aj  -  2xLj/Xc  sin  and  bjj  -  2irLj/xc  sin  X^,  and  pa 

and  are  the  incoming  phases  of  the  complex  envelope  of  the  desired  signal 
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and  the  1th  lnterferer  respectively.  Lj  is  the  distance  between  the  first 
and  Jth  input  element.  Xc  is  the  free  space  wavelength  of  the  centerband 
frequency  of  the  input. 

Using  (11)  in  the  definition  of  the  input  autocorrelation  maxtrlx  (Rjj) 
and  assuming  that  the  noise,  signal  and  the  m-interferers  are  uncorrelated  we 
can  say  that 


*x  -  *s  +  I  *l±  *  °n* 
i-1  1 

2 

where  oQ  is  the  noise  power,  I  is  the  identity  matrix,  and  $g 
defined  as 


and  tT  are 


(14) 


♦ 


s 


V 
8  8 


(15) 


8 1 


~  PM 

2  IiIi 


(16) 


Using  the  definition  of  Rx  (14)  in  (8)  we  get 


m 

dw( t)/dt  +  2k(*a  +  l  +  o„I)w(t)  - 


(17) 


Equation  (17)  describes  the  controlling  equation  for  each  of  the  element  weights. 

Using  the  scheme  of  Fig.  (4)  to  extract  a  reference  signal  (d(t))  from  the 
array  output  (v(t)),  and  assuming  wide  band  interference  with  a  narrow  band  de¬ 
sired  signal,  it  can  be  easily  shown  [2]  that  the  reference  signal  is  given  by 


d(t)  - 


a(P8w(t)) 

I P«w( t ) | 


i-H 


(18) 


where,  a  is  a  constant  depending  on  the  limiter  level  and  filter  attenuation  of 
Fig.  (4).  <(>  is  the  phase  shift  Introduced  by  the  reference  loop  and  is  depen¬ 
dent  on  the  frequency  of  the  incoming  desired  signal. 

Investigating  the  system  described  by  equations  (18)  and  (17)  with  no  in¬ 
terferences  ($1^*  [0]),  DiCarlo  [4]  demonstrated  that  the  weights  would  os- 
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dilate  at  a  radial  frequency  dependent  upon  the  reference  loop  phase  shift  ($), 
the  signal  and  noise  powers,  the  feedback  loop  gain  (k)  and  the  number  of  array 
elements. 

III.  Phase  Compensation  Scheme 

To  contend  with  the  problem  of  weight  cycling  caused  by  the  non-zero  phase 
shift  of  the  reference  loop  extraction  scheme  (Fig.  4),  Bar-Ness  [3]  proposed 
the  adaptively  controlled  phase  compensation  scheme  of  Fig.  (5).  This  scheme 
introduces  a  complex  weight  wr  to  the  reference  extraction  loop,  where  wr  is 
adaptively  controlled  so  that  it  minimizes  the  phase  difference  at  point  1  and  2 
of  Fig.  (5).  As  seen  this  controlling  is  done  by  the  use  of  an  integrator  and  a 
multiplier  which  together  form  a  correlation  loop  between  e^(t)  and  the  orig¬ 
inal  reference  signal  now  renamed  y(t).  From  Fig.  (5)  we  see  that  e}(t)  is  the 
difference  between  the  compensated  reference  signal  and  the  array  output  v(t). 
Thus  as  the  correlation  loop  minimizes  the  mean-square  of  ej(t)  the  phase 
difference  between  the  reference  signal  and  the  array  output  is  also  minimized. 


Analysis  of  Phase  Compensation  Scheme: 

Following  Bar-Ness's  analysis,  it  is  seen  from  Fig.  (5)  that  the  governing 
equation  for  the  added  complex  weight  wr(t)  is 
dw_(t)  _ 

— - -  -  -  2k  y*(t)e,(t)  (19) 

dt  r  1 

where  from  (18) 

a(P8tw(t)) 

y(t)  -  -  e”J$  (20) 

|P8cw(t)| 

and  from  Fig.  5. 

ej(t)  -  cy(t)wr(t)  -  xT(t)w(t)  (21) 

Noticing  that  the  newly  defined  reference  signal  d(t)  is  equal  to  the  sum  of 
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the  product  (yi(t)wr(t)]  and  an  amplitude  1 sited  version  of  (yi(t)wr(t))  we 


get 


y(t)  w_(t) 

d(t)  -  y(t)w  (t)  +  b,  I-—  ■ 

r  1  |y(t)wr(t)| 


(22) 


But  |y(t)|  -  a  which  is  the  first  limiter  level.  Thus  (22)  becomes 

b.  w_(t), 

d(t)  -  y(t)  fwr(t)  +  }- r  .1  (23) 

r  «|wrCt)| 

Now  that  the  reference  signal  (d(t))  is  defined  for  our  system  we  can  use  It 
to  describe  the  input-reference  cross  correlation  matrix  (R^  ■  E[x*(t)  d(t))). 
Substituting  from  (23) 


r  *  bl  , 

Rxd  -  E[*  *<*>!«,<*>  +  7[Wr('ty|']  1 

Substituting  for  y(t)  from  (20),  (24)  becomes 

R*d  «  E[x*(t)  [a  wr(t)  +  ^d^.]e-H] 

|PgW(t)|  |wr(t)| 


(24) 


(25) 


Using  the  assumption  that  the  Incoming  signals  are  uncorrelated  with  one  another 
we  see  that 


E[x*(t)p£]  -  AP>|  -  ^  *8 
✓2  A 


Using  this  in  (25)  we  get 


.  n  [  ♦.,.(t)i[.wr(e)  ♦ 

A  |PgW(t)|  | Wr( t) | 


(26) 


(27; 


Now  having  Rxd  defined  by  (27),  the  element  weight  equation  (17)  becomes 

+  2k(*a  +  I  *T  +  ajl)  w(t)  - 
dt  9  ll  n 

-2/jk  [  *flW(t)  w  (t)  +  Wjlje-J* 

A  |p|wCt)|  |wr( t) | 
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(28) 


Having  the  element  weights  defined  by  (28)  we  turn  our  attention  back  to 
the  phase  compensation  weight  wr(t).  By  substituting  equaiton  (21)  into  (19) 


we  get 

dwr(t)/dt  -  -2krE[y*(t)(cy(t)wr(t)  -  xT(t)w(t))J 
But  y*(t)  y(t)  -  |y(t)|2  -  a2 

And  by  using  the  definition  of  y(t)  from  (20),  (29)  becomes 

dwr(t)/dt  -  -2kr[ca2wr(t)  -  E[ (aPgW*(t)/ | Ppw( t) | )xc( t)w( t)eJ*l 

where:  '  stands  for  transpose  conjugate 
Since  P^w*(t)  is  a  scalar  we  can  see 
PgW*(t)  -  w'(t)P* 

and  using  (16)  we  can  rewrite  (29b)  as 

dwr(t)/dt  -  -2krfca2wr(t)  -  a/2/A  (w'( t)*8w(t)/ |PgW(t) | )e^*] 


(29a) 


(29b) 


(29c) 


,  2/2k_a  w'(t)*_w(t) 

dw  (t)/dt  +  2k  ca2w  (t)  -  — --2L.  [ — - ? - JeJ* 

r  A  |P|w(t)| 


(30a) 


Using  the  fact  that  9S  and  (i*l  to  m)  are  all  of  rank  1  and  complex 

m 


Hermetlan  matrices 


,  we  see  that  £  is  of  rank  <  m  and  using  Appendix 
i-1  1 


can  find  a  unitary  matrix  Q  such  that 

As  -  O’  *SQ  -  r*l  .  .  .  o" 
.  0 


ui 

Aj  -  Q'  l  *i  Q  •  <111  •  •  •  •  0 

i-1  1 


| . _|nrt-l  x  m+1  . 

0 . 0 


where 


qij  -  positive,  real 

+  pij 
«ij  -  qije  J 


1  -  j 
i  *  J 


Premultiplying  (28)  by  Q'  we  get 


m 

+  2k  Q'(*  +  l  *T  +  o2I)Q  O'  w(t) 
dt  8  A1  n 

-  2/2k  QjsQcrM  u  +  Wt)|,-J« 

A  |PgW(t)|  |Wr(C) | 

Defining  r(t)  -  Q'w(t)  -  Iyi.  Y2,...  Ynl 

Using  (31),  (32),  and  (35),  (36)  becomes 

dr(t)/dt  +  2k[A3  +  Ax  +  oJl]r(t) 


.  2 /2k 

A 


|p£w(t)| 


r(t)[a  w  (t)  + 

1  I  >  V  I 


r(t)l 


Evaluating  |PgW(t)|  we  get 

|PgW(t)|  -  [[PgW(t)]’[pJw(t)]]1/2 


[W' ( t)P3P3W( t) ] 1/2 


or 

|PgW(t)|  -  [w* ( t)Q  Q*P*p£q  Q’w(t)]1/2  -  [r'(t)-2  *8  r(t)l1/2 

Az 

*  ^  [r  (t)A_  r(t ) j 1/2 

a  8 

Using  the  definition  of  As  [(31)]  and  (36),  (39)  becomes 

|Psw(t)i  -  ^  i Iti i 2  ^ ! i 1/2  - 

A  A 


(35) 


(36) 


(37) 


(38) 


(39) 


(40) 


And  (37)  becomes 

dr(t)/dt  +  2k[ Ag  +  Aj  +  o2I]  T(t) 

Ai  b,w_(t) 

-  2k-_i -  r(t)[a  w  (t)  +  J-E - ]e  (41) 

^lKlI  |wr(t)| 

From  (31)  and  (32)  we  see  that  (41)  can  be  written  in  the  form 


dYi(t)/dt  +  2k  o2  Yi(t)  -  0 

m+1 

dYi(t)/dt  +  2k[  l  9ijYj  +  o2  Yi 1  *  0 
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i-M*2,  . N  (42) 

i-2,  ....  M+1  (43) 


m+l 

dy1(t)/dt  +  2k(  £  q.  Yj  +  (Xj  +  o2)yi 1 
j-1  J 

.  2„  /n  u»r(C)  *  b»“'U)u-J* 

I Y 1 ( t ) j  |wr(t)| 

Turning  our  attention  back  to  the  reference  loop  weight  wr(t),  we  see 
that  by  substituting  (40)  into  (30)  and  using  (36)  we  get 

,  r’(t)A_r(t) 

dwr(t)/dt  +  2krcaiiwr(t)  -  2a  kp  - 2— -  e^* 


(44) 


(45) 


IyiIA 


1 


or 


dwr(t)/dt  +  2krca2wr(t)  -  2a  kf  Yj  le^ ^  (46) 

At  this  point  it  is  noticed  that  eqs.  (42) ,(43) ,(44) ,  and  (46)  are  a  set 
of  N-complex  nonlinear  differential  equations  which  describe  the  response  of  the 
system  weights.  Looking  at  (42)  we  see  that  [N-(m+l)]  of  the  weights  will 
exponentially  decay  to  a  steady  state  value  of  zero.  Or  solving  (42)  we  see  that 

-2ko£t 

yt( t)  -  Yj(0)  e  n  i  -  M+2,  .  .  .  N  (47) 

To  solve  for  the  remaining  M+l  weights  plus  the  complex  weight  wr(t)  we 
define: 

-19, (t) 

Yl(t)  -  Pt(t)e  1  (48) 

wr(t)  -  Pr(t)e'J'i'(t:)  (49) 

where  P^(t)  and  Pr(t)  are  assumed  non-negative.  Using  (48),  (49),  along  with 
(33), (34)  we  can  rewrite  the  system  weights  equations  [ (43) , (44) , (46) )  as 

-jeACt) 


d[P1(t)e 


m+l  -j  (6j  (t)+pn  ) 

]/dt  +  2k(  l  q^j  Pj  ( t)e 


(50) 


*  »2p1<U.':l9*<t,|  -  0 

m+l 


i-2,  ...  M+l 


d[Pj(t)e  J9l(t>]/dt  +  2k[  £  qj"  P,(t)e  (0j(t)+plj) 

J-1  J 


+  (X„  +  a2)Pl(t)e 


-J  9 1  (  t ) 


]  -  2k/Xje 


“J  8 1  ( t ) 


(aPr(t)  +  bj]e 


-j ( ♦+♦) 


(51) 
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dtPr(t)e"J*(t)]/dt  +  21crca2Pr(t)e-J*(t)  -  2akr/x7P1(t)ej* 


Evaluating  the  derivatives  of  (50),  (51),  (52)  we  get 


dP1(t)_-je1(t) 
dt  C 


dt  j.J  J  J 

+  oJP1(t)e  -  0  1-2 . M+l 


dPi(t)  “ j 9 1 ( t )  -j0i(t)  d9t(t)  ®tl  +  “(01+p.i) 

— ~ — e  1  -  jPj(t)e  1  (■■  -L. ]  +  2k[^q|JPJ(t)e  J 

+  (Xt  +  o2)Pl(t)e  J0i(t)J 
-  2k/x7(e'j0l{t>J[aP  (t)  +  blJe‘J(*+4’> 

(! 

ffr<j^-J»(t)  -jP  (t)e-j«Kt)[di<lIj+2krca2Pr(t)e-^(c> 
dt  r  dt  r  r 

-  2akr^P1(t)eJ*  (! 

Multiplying  both  sides  of  (53), (54)  and  (55)  by  \  e^1^  \nd  e-J'f’^^ 

respectively  and  equating  the  real  and  imaginary  parts  of  each  we  get 

nrfl  +  ? 
dP1(t)/dt  +  2k[  l  qJjPj(t)cos(01  -  0j-  PlJ)  +  o‘P1(t)]  -  0 

1  -  2,  .  .  .  M+l  (! 

nrfl  +  ? 
dP1(t)/dt  +  2k l  l  qJjPj (1)008(0!  -  9j-  Plj)  +  (XL  +  o^P^t)] 

-  2k/\7[aPr(t)  +  bilcos  ($  +  *)  (! 

dPr(t)/dt  +  2krca2Pr(t)  -  2akr/\!P!(t)cos(^  +■  <»)  (' 

d0  (t)  nr+1 

Pj/t)  — i- —  -  2k{  l  q^jPj(t)sln(01  -  0j  -  pAj > J  -  0 

1  -  2,  ...  M+l  (! 
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(60) 


de,(c)  *+1  . 

Pj(t)- -  “  2k[  J  q^jPj(t)  sln(0j  •  0j  *  p j j  )  1 


-  2k/X[laPr(t)  +  bj]  sln($  +  i|i) 

Pr(t)itiil  -  -2ak  p/X^P^t)  sin($  +  i|i) 
dt 


(61) 


Looking  at  equations  (38)  and  (61)  we  see  that  the  compensation  weights  (wr(t) 
Pr(t)e“^^t^)  controlling  equation  Is  not  affected  by  the  interference.  Thus 
we  get  from  Bar-Ness  [3] 


cot  [$(  t)-hj>(  t)  ] 


Pr(0)cos(^»10+i(>Q)+2kra/Xn/QP1(t)exp(2krca2T)dT 


Pr(0)  sln(^io  +  (62) 

Also,  if  P^(t)  is  bounded  away  from  zero,  and  assuming  It  Is  positive  then 
for  sufficiently  large  t  we  have 

*(t)  =  -g»(  t)  (63) 

Thus  the  refernece  loop  phase  shift  is  compensated  for  and  the  weight  equations 

(56),  (57),  (58),  (59),  (60)  and  (61)  reduce  to 

dp. ( t)  “+1  .  , 

— i -  +  2k[  l  q|jPj(t)  cos(01  -  9j-  p^)  +  (Xj  +  onZ  ^(t) 


-  2k/Xi[aPr(t)  +  bj 

dp  (t)  m+1  2 

— - -  +  2k(  l  q^jPj(t)  cosCSj  -  0j-  Plj)  +  (on)P1(t)  -  0 


J-l 


1  -  2,...H+1 


dP_( t)  9  _ _ 

— - - +  2k  ca^P  ( t)  -  2ak  /xTPj(t) 

dt  r  r  r  i  i 

d9  (t)  m+1 

Pt(t) — i -  -  2k[  l  q^j Pj ( t)  8in(01  -  0j-  Pjj)!  *  0 


(64) 

(65) 

(66) 


1  -  2....M+1 


d0i(t)  m+1 

^(t) — i -  -  2k [  l  q|,P.(t)  sln(0j  -  0.-  pjj)) 

dt  j "1  ^  ^ 


p  rt)Mii  -  o 

r  dt 


(67) 

(68) 
(69) 
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At  this  point  we  have  shown  that  the  added  weight  wr(t)  effectively  coa pen- 
sated  for  the  reference  loop  phase  shift  $.  However ,  we  have  not  as  yet  been 
able  to  find  a  general  solution  for  the  weight  equations  (64-69).  In  the  next 
section  we  investigate  Che  special  case  of  a  single  Interference  and  give  a 
sufficient  condition  for  convergence  of  the  array  weights. 

IV.  Special  Case:  Single  Interference  In  an  N  Element  Array 

With  the  number  of  lnterferers  (a)  reduced  to  1  the  input  vector  (eq.(ll)) 
becomes 

x(t)  -  6.  p  +  I  p.  +  N(t)  (71) 

✓2  8  n  1 

where  A  and  B  are  the  amplitudes  of  the  signal  and  Interference  respectively. 

Ps  and  P j  are  their  respective  direction  vectors  and  are  defined  in  equa¬ 
tions  (12)  and  (13). 


Following  the  derivation  for  the  multiple  Interference  case  we  see  that 

A2  *  t 

we  must  diagonalize  the  desired  signal  autocorrelation  matrix  (*8*j  P8P8) 

as  In  (31).  Also,  we  must  use  the  same  unitary  diagonalizing  matrix  Q  in  order 

2  ^ 

to  transform  the  interference  autocorrelation  matrix  (*j«j  pipi^  Co  the 
form  of  (32)  where  m«l.  Using  the  definitions  of  Appendix  A  (AS),  let 


where:  V}  is  the  eigenvector  corresponding  to  the  only  non  zero  eigenvalue  of  t8 

v2  “  vn  are  th®  eigenvectors  corresponding  to  the  n-l  zero  eigenvalues  of  t8 
'  denotes  the  transpose  conjugate  of  the  vector 
By  inspection  we  can  see  that  a  proper  choice  for  V}  is: 

Vj '  -  P8C/|P8|  -  l/v'N  (1  e  J°2  ....  e  Jan)  (73) 
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(74) 


Using  this  knowledge,  we  get  from  (31) 
As  -  q'*bQ  -  A*  q'pIpJq 


A 

7 


v\ 

1 

[1  e*302  . Jan] 

f 

jct2 

v  2 
• 

e 

e 

• 

*  • 

a 

>n 

v  n 

Combining  (72),  (73),  and  (75)  we  see  that 

X  j  0  •  •  • 

is  -  .  i  0  0 

where  from  Appendix  A  (A4)  we  know 
,  _  NA2 


•  • 

•  • 

V1  v2 
•  • 


(75) 


(76) 


(77) 


Using  Appendix  A  we  see  that  we  now  can  place  a  restriction  on  Q  such  that 

i2  i  *  *  t»2  * 


o’i 


'pipt«  '  f 


z  zl 


(78) 


where : 

Zc  *  [<VjPj>  *'V3^I^  •  •  •  •  <vnPj>] 


(79) 


(80) 


From  Appendix  A  we  know  that  this  can  be  accomplished  if  we  let  vj  and  V2 
span  the  plane  defined  by  P*  and  P*  and  making  the  other  n-2  eigenvector 
orthonormal  to  the  space  spanned  by  vj  and  V2»  Thus  using  the  diagram  below 
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we  see  Chat 


/q 


11 


p*  -  <p;,v1>v1 

Ip?  -  <p?Vvil 


Using  (79)  and  (80)  we  get 


<in  -fi<p!.vi>i2  -|^i<p*i.p;>i2 

where  we  have  also  used  Che  fact  that  |P8|2  ■  P8P*  “  N. 

Furthermore,  using  (79)  and  (80)  In  conjunction  with  Appendix  B  we  get 


(81) 


(82) 


B 


2  ■  *  v  1 2 


<*22  *  J-I<PI»  v2>l‘  "  ~  '  qH 


<*21  "  <»12  “  f-  <PI»  V  <VP  Pl> 


(83) 

(84) 


From  Appendix  C  we  know  that  Is  real  and  thus  (84)  Is 

seen  as 

l<»2l!  "  l<»  12 1  ’  f~  <PI»V2>  l<vl»pI>l 
/<l21  “  /<»12  “  /«g<v1’PI> 


(85) 

(86) 


Noting  that  we  can  make  Aj  a  real  matrix  If  we  eliminate  the  angle  Intro¬ 
duced  In  (86),  by  Inspection  using  (79)  we  see  that  we  can  define  a  new  Q  such 
that 

Q„'  -  A'Q*  (87) 
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1 


where 


"l/arg(<vi  ,PT*>) 


A’  - 


0 


0 


1 


Now  using  Q'n  we  define  Ax  as 


Anl  “  Qn'*I  Qn 


9nll  9nl2 

0 

9n21  9n22 

0 

0 

(88) 


(89) 


where 

Onll  “  <111  "fjf  l<PI»P**>l2  (90) 

9n22  "  922  "  l<v2»PI*>l2  “  ^  "  <Ul  (91) 

9n21  “  9nl2  *  ^11^22  “  l<112  I  (92) 

Obviously  Anx  as  defined  In  (89)  Is  a  positive  real  symmetric  matrix. 

Also  we  realise  from  (75)  that 

Ana  “  Qn'^s^n  “  Q'*sQ  "  As  (94) 

As  we  did  In  the  multiple  Interference  case  we  let  r  ■  Qn'w(t).  And,  using 
(90)-(94)  we  see  that  one  Interference  causes  the  multiple  Interference  system's 
freight  aviations  (64)-(69)  to  become 


dPj(t) 

dt 


dP2(t) 

dt 


+  2fc((q11  +  +  oJ)Pj(t)  +  /q11q22  cos(e1-02)P2(t)] 

-  2k  /TTla  Pr(t)  +  b x J 
+  2k(/qnq22  Pj(t)  cos(02(t)  -  0j(t)) 

+  (^22  +  ®n^P2^t^  "  0 


(95) 


(96) 
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+  2kfca2Pr(C)  •  2a  lcr  Pj(t) 


(97) 


dPr(t) 

dt 


Pt(t) 

P2(t) 

Pr(t) 


d81(t) 

dC 

d02(t) 

dt 

dipC  t) 
dt 


+  2kl  /q22qu  aln(e2(t)  -  6j(t))P2(C)]  -  0 
-  2k[  /qnq22  sln(o2(t)  -  e1(t))P1(c)]  -  0 

■  0 


(98) 

(99) 

(100) 


Trying  to  analyze  the  convergence  of  flj(t)  and  fl2(t)  we  can  preaultlply 
(98)  and  (99)  by  P2(t)  and  P^(t)  respectively  and  substract  to  get 

P1(t)P2(t)^il-^i)-  2k  /qltq22  [P12(t)+P22(t)  1  sln^-Oj)  -  0  (101) 


or 


d(02-01)_  2k/quq22[P12(t)  ♦  P22(t)I 


sin(02~O| ) 


dt  Pl(t)  ?2(t)  *  *  (102) 

Assuming  that  Pj(t)  or  P2(t)  cannot  attain  a  valve  of  zero  over  a  finite  time  in¬ 
terval  then  the  phase  difference  02(t)  -  0i(t)  will  converge  on 

02(t)  -  9l(c>  “  *•  003) 


Using  (103)  in  (98)  and  (99)  we  get 

P2(t)  d02(t)/dt  -  0  or  d02/dt  •  0  (104) 

Pj(t)  d@i (t)/dt  -  0  or  d0i/dt  -  0  (105) 

It  can  be  shown  from  (95),  (96),  (97)  by  inspection  that  In  fact  Pi(t),  P2(t), 
Pr(t)  can  never  attain  a  steady  state  value  of  zero  and  thus  equations  (100), 
(104),  (105)  hold. 

At  this  point  we  have  shown  that  the  phases  of  the  weights  f(t),  0i(t), 

02(t)  converge  on  a  steady  state  value.  Mow  using  equation  (103)  in  (95), 

(96),  and  (97)  we  get  a  linear  set  of  equations  for  the  magnitude  of  the  weights. 
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(106) 


— +  2kl(qn  +  Xj  +  a^P^C)  -  /q11q22P2(t) 

at 

-  2k  /I7[aPr(t)  +  bx) 

dP2(t) 


dc 


-  2 k [  1^22^1  ^22  ®n^2^^  "  0 


dpr(t)  2 

— i +  2k_caz 


dc 


2krca‘‘Pr(t)  -  2akr  /XjPjO:) 


(107) 


(108) 


or 


PL(C) 

P2(c) 

- 

P3(C) 

-2k((qu  +  XL  +  on)} 


+2k  /qnq22 

-2akf 


+2k[  */,<Ixi<l22^ 

-2k(q22  +  o„) 
0 


-2ka/I[ 

0 

-2krca2 


Pt(t) 

P2(t) 

P3(t) 


2k/Xj  t>i 
0 
0 


(109) 


The  eigenvalues  of  Che  maCrlx  in  (109)  can  be  obCained  by  solving  Che  equaClon 
(X  +  2krca2)[(X  +2k( Xj  +  a*  +  qn>)(X  +  2k(X2  +  o*  -  qu)) 


-  4kZqn(X2  -  qu)] 
-  4kkrXja2(X  +  2k(X2  +  a*  -  q 3 1 ) >  *  0 


(110) 


un2  Ub2 

when  we  used  Che  face  ChaC  q22  ■  -  q^  and  defined  X2  •  (HO)  can  be 


rearranged  Co  geC 


(X  +  2k(X2  +  o„))I(X  +  2kfca2)(X  +  2k(x^  +  a„))  -  4kkrx1a2) 

-  4k2qnX1(X  +  2kra2(c-l)>  -  0  (111) 

(111)  can  be  looked  upon  as  an  equaClon  for  Che  loci  of  Che  ehree  eigenvalues 
when  qu  varies  from  0,  when  Che  inCerference  direcdon  is  orChogonal  Co  the 
desired  signal  direcdon,  Co  a  value  equals  X2>  when  boch  direcdon  are  Che 
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0  then  by  (82)  qjj  is  also 


same.  Notice  that  with  no  Interference  X2  »  -  “ 

*■  2 

zero  and  (111)  reduces  to 

(X  +  2kof)[(X  +  2k_ca2)( X  +2k(X,  +  of))  -  4kk rX,a2J  -  0 

(112) 

This  is  exactly  what  we  had  In  [3].  Also  If  for  the  orthogonal  case  qu  “  0  but 
X2  *  0  then  froa  (111) 

(X  +2k(X2  +  of)) l(X  +  2k  ca2)(X  +2k(X,  +  of))  -  4kk  XiS2J  *  0 

(113) 

so  that  for  this  case  the  Interference  has  the  effect  of  changing  one  of  the 

2  2 

el con values  from  -2kon  to  -  2k(X2  +  on).  ®ne  can  show  that  the 

other  two  eigenvalues  an  unconditionally  real  and  they  are  both  negative  when 
c(X,+  of) 

- - - —  >  1.  For  ■  X2  we  get  after  rearranging  the  teras  of  (111).  (See 

Appendix  D) . 

(X  +  2kof)[(X  +  2k_ca2)(X  +2k(X,  +  X2  +  of))  -  4kkrA,a2J  -  0 

(114) 

Again  it  can  be  easily  shown  that  the  two  eigenvalues  that  correspond  to  the 

second  factor  are  unconditionally  real  and  both  ate  negative  If 

c(*l  f_X?...-.qn.)  >  1  (115) 

*  1 

The  condition  Is  obviously  easier  to  satisfy  than  the  previous  condition  we 
obtained  with  qu  -  0.  Comparing  (113)  with  (114)  we  notice  the  difference 
in  the  second  factor  of  the  first  term  in  the  square  parenthesis  where  we 
have  Xj  +  X2  representing  the  sum  of  power  of  both  the  interference  and 
the  design  signal  Instead  of  Xi  representing  the  power  of  the  desired  signal 
only.  The  first  terms  in  the  parenthesis  remain  unchanged  since  the  reference 
generating  loop  is  not  affected  by  the  Interference. 

In  order  to  get  an  idea  of  what  happens  to  the  eigenvalues  as  the  Inter¬ 
ference  direction  varies  between  the  orthogonal  and  colinear  directions  of  the 
desired  signal  we  can  plot  the  root  loci  of  (111)  keeping  in  mind  that  qu 
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F 


increases  from  zero  Co  X2  as  the  direction  of  Che  interference  moves  from  an  or¬ 
thogonal  to  a  colinear  direction.  Doing  this  we  can  prove  that  the  eigenvalues  will 

c(X.+  of) 

always  have  a  negative  real  part  as  long  as  _  ,  >  >1.  The  proof  is  given  in 

>1 

Appendix  E. 

The  steady  state  values  of  the  weights  are  derived  by  setting  the  left- 
hand  side  of  (109)  to  zero  and  evaluating.  Doing  this 


lim 

t-w 


Pi  ( t> 
P2(t) 
P3U) 


_ 2 _ Zl  --  _ 2 _ 

2klq22(c(on  +  xp  -  Xj)  +  ®n(c^ll  +0n  +  *l)“  *j)J 


c[q22  +  <Jnl  +c  l/9n922 


~ *  /xl(q22  +  of] 


+c  ^2^22  +  °n+9ll)'XlJ 


/XT 


kra 


/xlqllq22 


1  y'~/’  "  k  [(q22+  <’n)(xl+  °n)+  °nq22J 


fq22  +  On1  7  /Xlqllq22  2 

kra 


-2k  /Aibx 

0 

0 


(121) 


or 


a.  P.(t)  - 

1  Cl 


(122) 


llm  P 
t  ><» 


(t)  -/V:bl/qllq22 


xlbl 


lio  P_(t)  -  —1 
t~  f  G 


lq22  +  °n^ 


(123) 


(124) 
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where 


G  "  q22 Ic<an  +  Al>  -  *1*  +  +ffJ  +  *0~  M1  025) 

To  conclude  we  notice  thet  under  the  conditions  specified  the  added  coa- 
plex  weight  wr(t)  convergeson  a  steady  state  value  which  cancels  the  reference 
loop  phase  shift  +  .  By  doing  this  we  have  shorn  that  the  result  Is  the  elimi¬ 
nation  of  the  weight  cycling. 

He  have  also  shown  that  an  Increase  In  the  Interference  power  will  one 

of  the  eigenvalues  aore  negative,  thus  giving  a  stabilising  effect  to  the 
systea.  In  general,  we  see  that  the  array  weight's  convergence  rates  depend  on 
the  power  as  well  as  the  relative  directions  of  the  Interference  and  desired 
signal. 

V.  The  Uncompensated  Scheae  With  Single  Interference 

The  uncoapensated  scheae  will  be  considered  as  a  special  case  of  the  com¬ 
pensated  scheae  with  Pr(t)  -  1  and  <»(t)  -  0.  Using  this  In  (56)  and  (57) 
respectively  we  get 

dPj(t)  2 

— — —  +  2k(q11  +  Xj  +  on)Pi(t)  +  2kq12P2(t)  cos(e1(t)  -  02(t)) 

“  2 lea  /\7  cos  4  (126) 

dP. ( t )  2 

— ~  +  2kq21Pl(t)  cos(02(t)  -  0j(t))  +  2k(q22  +  on)P2(t)  -  0  (127) 

Similarly  (60)  and  (59)  become 

d0i (t)  _ 

— +  2k9i2p2(t)  "  82(t))  -  2ka  /Xj  sin  +  (128) 

d0o(t) 

P2(t)— -  2kq2lPL(t)  sln(0l(t)  -  0 2 ( t ) )  -  0  (129) 

A  general  solution  for  the  differential  equations  of  (126) ,( 127) ,( 128),  and 
(129)  has  not  yet  been  found.  However,  for  the  case  when  the  direction  vector 

F-21 


of  the  lnterferer  end  the  desired  signal  are  othogonal  (l.e.,  <PjPs>  •  0) 
we  can  reduce  the  problem.  For  this  case  we  get  from  (82),  (83),  and  (84) 
that 


dll  -  <121  *  <112  "  0 

(130) 

_  NB2 

<22  "  5J- 

Substituting  (130)  and  (131)  into  (126),  (127),  (128), 

and  (129)  we  get 

~1——  +  2k(X,  +  o*)P,(t)  -  2 lea  /x7  cos  6 
dt  1  1 

(131) 

2k(q22  +  P2(t) 

(133) 

d6 1 (t)  . — 

P,(t) — i - •  2ka  /X,  sin  6 

1  dt  i 

(134) 

P2(t)—§ - -  0 

dt 

(135) 

From  the  solution  of  (133)  we  see  that  the  weight  y2(c)  *  P2(t)exp[-J02(t) ) 
will  exponentially  decay  to  a  steady  state’ value  of  zero.  Also  we  note  that 
(132)  and  (134)  are  the  equivalent  to  the  weight  equations  derived  by  DeCarlo 
[4]  for  the  interference  free  case.  (n.b.  equations  (30)  and  (33)  of  reference 
[4]).  In  his  solution  DeCarlo  showed  that  the  magnitude  of  the  weight 
( I "V 1 ( t ) |  *  Pi(t))  would  converge  on  a  steady  state  value  proportional  to  the 
cosine  of  the  phase  offset  However  the  phase  of  the  weight  (arg(Yj( t))« 

01 (t))  will  oscillate  at  a  radian  frequency  proportional  to  the  tangent 
of  the  phase  offset 

To  conclude  we  see  that  when  the  direction  vector  of  the  Interference  Is 
orthogonal  to  that  of  the  desired  signal,  the  Interference  will  have  no  effect  on 
the  system  weights.  And  we  see  that  the  weights  will  approach  a  limit  cycle  in 
exactly  the  same  fashion  as  DiCarlo  [ 4 J  predicted  for  the  no  interference  case. 


VI.  Conclusion 


In  this  paper  we  have  devised  a  method  by  which  we  have  been  able  to  better 
study  the  effects  of  interference  on  the  array  weights  of  an  adaptive  array  sys¬ 
tem.  Using  this  method  we  have  further  analyzed  Bar-Ness's  compensation  scheme 
[3]  and  have  shown  that  for  a  multiple  Interference  case  the  reference  loop  phase 
shift  is  effectively  compensated  for.  We  have  also  been  able  to  show  that  for 
the  single  Interference  case  the  array  weights  will  converge  on  a  steady-state 
value.  For  the  uncompensated  Compton  scheme  [2]  it  was  demonstrated  that  the 
array  weights  will  cycle  for  the  special  case  of  one  interference  with  a  direction 
vector  orthogonal  to  that  of  the  desired  signal.  Future  work  is  planned  for  the 
study  of  a  general  solution  for  the  multiple  Interference  case  with  Bar-Neas's 
Compensation  Scheme  as  well  as  a  general  solution  for  the  single  lnterfernece  case 
with  Compton's  scheme.  Furthermore,  analysis  is  planned  for  the  signal  to  noise 
plus  interference  ratio  of  Bar-Ness's  scheme  in  a  single  interference  environment. 
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eigenvectors  where  V}  is  the  eigenvector  corresponding  to  the  one  nonsero 
eigenvalue  of  *8  and 


vl-Vlpsl*  <A6> 

Also,  v^(i«2  to  N)  are  the  eigenvectors  corresponding  to  the  N-l  aero  eige¬ 
nvalues  of  4g. 

Now,  turning  our  attention  to  the  M  interferers,  we  see  that  the  autocor¬ 
relation  matrix  (i  -l  to  M)  of  each  interference  is  defined  in  (16)  as 


4 


(A7) 


Each  matrix  4j  (i  -  1  to  M)  is  Hermetian  of  rank  1.  Thus  we  see  that  a  matrix 
M  1 

defined  by  £  will  also  be  Hermetian  and  have  rank  «  H.  Now  leaving  vt 
1-1  1 

alone  so  that  (A2)  is  always  satisfied  we  can  see  that  the  vectors  vt  (i  -  2 

to  M+l )  can  always  be  made  so  that  vj(J  -  1  to  M+l )  will  span  the  M-space 
M 

spanned  by  £  .  And,  the  eigenvectors  v^fk  -  M+2  to  N)  can  be  made  ortho- 

1-1  1 

gonal  to  that  M-spac*  as  well  as  the  direction  of  V}.  By  doing  this  we  can  get 

M 

a  transformation  of  :he  matrix  £  4j  such  that 

1-1 


M 

Ai  -  o*  i  no- 

i-1 


0 


M+l  M+l 
» 


(A8) 


where  since  this  a  similarity  transformation  Ax  must  be  Hermetian  and  there¬ 
fore 


F-25 


qjj  -  positive  and  real  i  ■  J 

dij  -  qji  -  complex  i  *  J 

For  i  *  J  we  can  write  q^j  in  the  form 

qij  "  <lij  e  J 


where 


«ij  -  hij  I 

Pij  “  arglqij 1 


(A9) 

(A10) 

(All) 

(A12) 

(A13) 


1 

1 

1 


Appendix  B 


From  (83) 


**22  “  l<pI*»  v2>l‘ 


Substituting  for  V2  from  (81) 


or 


>i2 

2  iPl^P^.v^Vjl 

q  .  B2  <Pl*.Pl*>  -  <Pl*.vI>*<PlV|2 

22  2  |PT*  -  <?!*,  vj  >  Vl| 


Evaluating  the  demoninator 


|Pj*  -  <PI*,v1>v1j  -  I<Pj*  -<PI*,v1>v1,  Pr*  -<PI*,v1>v1>l7 


[<?!*,?!*>  -  <PI*,vlXv1,PI*> 

-  <Pj*,V|>*<Pj*,Vj> 


+  <PI*,v1XPI*,v1>J 


or 


I  pi  -  <PI*,vl>vl|  -  (<PI*,PI*>  -  <PI*,v1Xv1,PI*>) 
Substituting  (B6)  Into  (B3)  we  get 

q22  “  Y  I<PI*,PI*>  ~  <PI*,v1>*<PI*v1>J 


|-  t<PI*,PI*>  -  Kp^.v^!2) 


But  <Px*,Px*>  -  N 
and  from  (82) 


Y  I<pi**vi>I2  “  *111 


Thus  (B8)  becomes 


^22  "  ~2~  qil 


(Bl) 

(B2) 

(B3) 

(B4) 

(B5) 

(B6) 

(B7) 

(B8) 

(B9) 

(BIO) 

(Bll) 
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Using  the  definition  of  V2  froa  (82)  we  get 
<P 


*  v  _  -  *  PI  -  <*1  ,V1>VK 
I  ,V2>  -  <Pj  ,  A  * > 

|pl  -<Pj  .V^Vjl 


or 


_  *  v  ”1  *rI 
<?I  ,v2>  - 


<PT*,PT*>  -  <P*,v.>*<P  *,v.> 


(Cl) 


(C2) 


1*1*  ~  <pi*  »vl>vl I 
Substituting  for  the  demoninator  from  (B6)  and  using  (BIO)  we  get 

1 

<PT*,v5>  -  I<PT*,PT*>  -  <PT*,v,>*<PT*,v,>]7 


i  *vr  *vr 

l 

-  IN  -  (<?!*, Vj>|2J7 

Thus  froa  (C4)  we  see  that  <Pj*,v2>  will  always  be  a  real  nuaber. 


(C3) 

(C4) 
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Substituting  qu  -  *2  in  OH),  we  get  after  rearranging  terms 
(A  +  2krca2)[(X  +  2k(A2  +  <j*)>0  +  2^!  +  o*))  +  4k2A1A2J 

-  4kkrA1a2(A  +  2k(A2  +  oj))  +  8k2A1A2kra2  -  0  (Dl) 

or 

(A  +  2krca2)[A2  +  2k(At  +  A2  +  2o*)A  +  4k2oJ(Xi  +  A2> 

+  4k2(o„)2J  -  4kkrAta2(A  +  2k a*)  -  0  (D2) 

i 

Factoring  the  first  and  rearranging  we  finally  get 

t 

I  (A  +  2ko*)[(A  +  2k_ca2)( A  +  2k(A,  ♦  A2  +  o*))  -  4kk_  A,a2]  -  0 

|  nr  (D3) 

i 

i 

I 

| 


i 


\ 

i 


Appendix  E 


We  know  by  solving  (113)  that  the  roots  of  (111)  for  qn  ■  0  are  uncon¬ 
ditionally  real  and  negative  If  c(Xj  +  o^)/Xj  >  l.  Thus  under  this  condi¬ 
tion  the  root  loci  of  (111)  has  3  poles  and  one  zero  all  on  the  negative  real 
axis.  We  are  only  Interested  in  the  roots,  for  0  <  qjj  <  X2  an<*  we  realize  from 
(114)  that  the  roots  of  (111)  as  qn  approaches  X2  will  be  unconditionally 
real  and  negative  under  the  previous  specified  condition.  Below  are  the  4 
possible  paths  which  the  root  loci  of  (111)  can  follow  according  to  the  above 
specifications 


where  Che  arrows  point  in  Che  direction  of  Increasing  qu«  The  poles  for  q  |  ^  •  0 
are  defined  by  Che  roocs  of  (113)  and  froa  (111)  Che  zero  is  aC  -2kra  2(c  -  l). 

Obviously  since  we  see  from  Che  above  root  loci  that  if  for  Che  maximum 
value  of  qn  “  X2  Che  roots  of  (111)  are  negative  and  real  so  ChaC  for  any  value 
of  qjj  <  X2  Che  root  must  always  have  a  negative  real  part. 

As  an  example  let's  assume  Che  following  values 


c  » 

1.2 

k  -  50 

£ 

°n  “ 

.001 

a  “ 

8.33 

kr  -  25 

b  * 

1.6262 

Xi  -  1 

x2  - 

1/2 

From  (113)  Che  equation  for  Che  roots  when  qu  -  0  is: 

(X  —50.1 ) [ (X  +  4166.33)(X  +  100.1)  -  347,194.45)  (El) 

which  has  roots  at 

X  -  -50,  -16.4,  -4249.9 

Nocice  also  Chat  Che  zero  for  this  root  loci  is  found  froa  (111)  to  be  -694.43. 

As  we  increase  qn  to  X2  •  1/2  from  (114)  we  get  roots  at  X  »  -.1,  -65.4,  -4251. 
The  root  loci  is  plotted  as 


-4249.9  -694.43  -50  -16.4  qn=0 

11  11 


t  it, 

-4251  -65.4  -0.1  qn=£ 


Notice  from  (113)  and  (El)  that  an  Increase  in  Interference  power  will  only 
affect  the  eigenvalue  at  X  ■  -50.  In  fact  an  Increase  in  power  will  cause 
this  eigenvalue  to  become  more  negative,  thus  helping  to  stabilize  Che  system. 
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Figure  2  Main  Feedback  Loop  for  LI 
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Figure  4  Reference  Signal  Loop 
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List  of  Symbols  -  Appendix  F 


A  desired  signal's  amplitude 

a  constant  dependent  on  the  limiter  level  and  filter  attenuation  of  the 

reference  extraction  loop 
ic^  interference  amplitude 
b}  compensation  scheme  limiter  level 

phase  of  the  ith  Interference  at  the  element 

c  compensation  loop  weighting  factor 

d(t)  reference  signal 

E[  ]  denotes  expected  value 

e(t)  error  signal  for  array  LMS  algorithm 

ej(t)  error  signal  of  compensation  scheme 

Vt  gain  constant  for  LMS  algorithm  in  array  weights 

1cr  gain  constant  for  LMS  algorithm  in  compensation  loop 

Lj  distance  between  l8t  and  jtl1  antenna  element 

M  number  of  interferes 

N  number  of  array  elements 

P^(t)  magnitude  of  linear  transformed  array  weights 
T 

Pj!  i  interference  Incoming  phase  vector  ( 1  x  N) 

Pr(t)  magnitude  of  the  complex  compensation  weight  wr(t) 

Pg  desired  signal's  phase  vector  (1  x  N) 

Q  diagonalizing  matrix  [N  x  N] 

components  of  Aj 
q+jj  magnitude  of 

R*  input  signal  autocorrelation  matrix 

Rx<j  input-reference  signal  cross-correlation  matrix 

t  denotes  transpose 

(Continued) 
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(Continued) 


V£  -  eigenvectors  of  $s  used  as  vector  components  of  Q 

wT(t)  -  complex  antenna  element  weight  vector  (1  x  N) 

w^(t)  -  components  of  wT(t)  corresponding  to  the  I6*1  complex  array  weight 

wr(t)  -  complex  reference  extraction  loop  compensation  weight 
xT(t)  -  system's  Input  vector  (1  x  N) 

Xjj(t)  -  inphase  component  of  the  lc^  array  element  Input 
xiQ(t)  -  quadrature  component  of  the  l*"*1  array  element  Input 
y(t)  -  Input  to  the  compensation  scheme 

aj  -  phase  of  desired  signal  at  the  jth  element 

rT(t)  -  O'w(t) 

-  uncoupled  complex  array  weight  vector  (1  x  N) 

-  components  of  r1(t)  corresponding  to  the  1  complex  transformed 
array  weight 

Gi  -  phase  of  linear  transformed  array  weight 

As  *  Q*  *s  Q 

-  diagonalized  desired  signal's  autocorrelation  matrix  (N  x  N) 

M 

A i  -  [  Q'  *u  Q 

1-1 

-  sura  of  the  diagonalized  Interference  signals  autocorrelation 

matrices  (N  x  N) 

2 

Ai  -  NA  /2  ;  only  nozero  component  of  A 

1  s 

*2  "  ;  trace  of  h 

angle  of  the  complex  value  q^j 

phase  shift  Introduced  by  uncompensated  reference  signal  loop 
lc>1  Interference  autocorrelation  matrix  (N  x  N) 
desired  signal's  autocorrelation  matrix  (N  x  N) 
phase  of  the  complex  compensation  weight  wrt 


Pij 

♦ 

*11  ' 

*8 

<»(t) 


F-39 


(Continued) 


END 

DATE 

FILMED 

-83' 

DTIC 


